por medio del metodo del trapecio calcular el valor solucion a una integral de un intervalo numero dado

program metodo_trapecio
!carmona rodriguez andrea elizabeth
!3IV5
!6/10/10
IMPLICIT NONE
INTEGER::N,K
REAL::A,B,H
REAL::SUMA,INT,F,X0
REAL,ALLOCATABLE::X(:)
f(X0)=X0**2
WRITE(*,*)"DAME EL No Total de intervalos"
read(*,*)N
WRITE(*,*)"DAME LOS VALORES DE A Y B"
READ(*,*)A,B
H=(B-A)/N
ALLOCATE (X(0:N))
X(0)=A
X(N)=B
DO K=1,N-1,1
X(K)=X(K-1)+H
END DO
SUMA=0
DO K=1,N-1,1
SUMA=SUMA+F(X(K))
END DO
INT=H/2*(F(X(0))+2*SUMA+F(X(N)))
WRITE(*,*)"LA INTEGRAL ES=",INT
DEALLOCATE(X)
END PROGRAM