Método de los trapecios

Otro método de aproximación de la integral mediante la suma de pequeños trapecios.

Programa Principal:

program integral1

use funcion

use numerico

implicit none

real::Xi,Xd,n,I

write(*,*)'elige un intervalo de integracion'

read(*,*)Xi

read(*,*)Xd

write(*,*)'el intervalo es:',Xi,Xd

write(*,*)'elige un numero de subintervalos'

read(*,*)n

call integral(f1,Xi,Xd,n,I)

write(*,*)'el area de la integral es:',I

end program

Module Numerico:

Module numerico

contains

subroutine integral(f,a,b,n,I)

interface

function f(x)

real,intent(in)::x

real::f

end function

end interface

real,intent(inout)::a,b,n

real,intent(out)::I

real::x1,x2,h,A1,A2,j

I=0

do j=1,n

  

x1=(a+((j-1)*ABS(a-b)/n))

x2=(a+(j*(ABS(a-b)/n)))

h=((f(x2))-(f(x1)))

A1=((x2-x1)*(f(x1)))

A2=((A1))+(((x2-x1))*((h)/(2)))

I=I+A2

end do

end subroutine integral

end module numerico

Module funciones:

module funcion

contains

function f1(x)

real,intent(in)::x

real::f1

f1=((x**3)+(4*(x**2))-10)

end function

end module funcion