Integración por el método del trapecio en MATLAB

Para implementar el método del trapecio en MATLAB se deben crear los siguientes scripts. El siguiente es solo un base, es posible complementarlo para mejorar la experiencia. 

A continuación se muestra el código que será responsable de desarrollar el proceso

function integral=Metodointegracion(a,b,nintervalos,fun)

ancho=((b-a)/nintervalos) 

integral=0;

for i=(1:nintervalos)

    A=(ancho*fun(a+(i*ancho)));

    B=(ancho*fun(a+(i-1)*ancho));

    integral=(integral+((A+B)/2))

end 

end 

En un script diferente se debe implementar el siguiente código, será el encargado de recibir los datos del usuario. (El codigo tiene los datos de la integral que resolvimos anteriormente)

nintervalos=4;

a=1;

b=5;

fun=@(x)(1./x)

integral=Metodointegracion (a,b,nintervalos,fun)

Nota: Asegúrate de guardar ambos scripts en la misma carpeta, de lo contrario MATLAB lanzará un error. 

Si quieres saber más sobre le procedimiento en papel, revisa nuestra entra de la sección "Métodos numéricos en papel"