´SIMPSON´
Comments
Content
C
C
C
PROGRAMA PARA INTEGRAR CON EL METODO DE SIMPSON
DIMENSION AREA(2)
100 FORMAT (I5,2F10.4)
101 FORMAT (1H1)
102 FORMAT ('INTERPOLACION DE SIMPSON PARA UN 1/3 ENTRE ',F10.4,'
Y
',F10.4/
1
'ANALIZANDO',I5,' SUBINTERVALOS',F18.4)
103 FORMAT ('INTERPOLACION DE SIMPSON PARA 3/8 ENTRE ',F10.4,' Y',F10.4/
1
'ANALIZANDO',I5,' SUBINTERVALOS',F18.4)
Y(X)=(X**3-12*X**2-4*X+336)/48
READ (5,*) N,A,B
IF
(N.LT.1)
CALL EXIT
H=(B-A)/N
X=A
SUMA1=0.
SUMA2=0.
M=0.5*(N-2)
P=0.5*N
DO 2 I=1,M
DO 3 J=1,P
3 SUMA2=SUMA2+0.5*Y(A-H+2*J*H)
2 SUMA1=SUMA1+Y(A+2*I*H)
AREA1=H*(Y(A)+Y(B)+2*SUMA1+4*SUMA2)/3
WRITE (6,102) A,B,N,AREA1
WRITE
(6,*) Y(A)
WRITE
(6,*) Y(B)
WRITE (6,*) SUMA1
WRITE (6,*) SUMA2
SUMA3=0.
SUMA4=0.
R=(N-1)/3
Q=N-1
DO 4 L=1,R
DO 5 T=1,Q
5 SUMA4=SUMA4+Y(A+T*H)
4 SUMA3=SUMA3+Y(A+3*L*H)
AREA2=0.375*H*(Y(A)+Y(B)-SUMA3+3*SUMA4)
WRITE (6,103) A,B,N,AREA2
WRITE
(6,*) Y(A)
WRITE
(6,*) Y(B)
WRITE (6,*) SUMA3
WRITE (6,*) SUMA4
END
SOLUCION:
6
2
14
INTERPOLACION DE SIMPSON PARA UN 1/3 ENTRE
2.0000 Y
ANALIZANDO
6 SUBINTERVALOS
48.0000
6.000000
14.00000
4.000000
11.00000
14.0000
INTERPOLACION DE SIMPSON PARA UN 8/3 ENTRE
2.0000 Y
ANALIZANDO
6 SUBINTERVALOS
48.0000
6.000000
14.0000
14.00000
1.000000
15.00000
Sponsor Documents