Programa 7.1 (Regla Compuest del trapecio - pagina 393).

%Programa 7.1 (Regla Compuest del trapecio - pagina 393).
% Construccion de la aproximacion a la integral
%
% b M-1
% I f(x)dx = h/2 (f(a)+f(b)) + h Sum f(x)
% a k=1
%
% evaluando f(x) en los M + 1 nodos equiespaciados x=a+kh,
% para k=0,1,2...,M. Notese que X<0>=a y que x=b.
%
%Sintaxis
% traprl(f,a,b,M)
%
%Datos
% - f es el integrando, dado como una cadena de caracteres 'f'
% - a y b son los extremos inferior y superior del
% intervalo de integracion
% - M es el numero de suintervalos
%Resultados
% - s es la aproximacion abtenida con la regla compuesta del trapecio


h=(b-a)/M;
s=0;

for k=1:(M-1)
x=a+h*k;
s=s+eval(f);
end
x=a;
fa=eval(f);
x=b;
fb=eval(f);
s=h*(fa+fb)/2+h*s;