ALGORITMO MATLAB

CODIGO MATLAB

clear all; 
clc; 
fprintf(' IIngresar la funcion\n\n'); 
   f=input('Metodo Trapezoide:','s'); 
   fprintf(' ');
     a=input('Tiempo minimo (h):'); 
     b=input('Tiempo maximo (h):'); 
     c=input('numero de segmentos a dividir:'); 
    h=(b-a)/c; 
  z=0; 
for x=a:h:b
k=eval(f);
if x==a,d=k;
end
if x==b,e=k;
end
z=z+k;
end 
fprintf('Resultado ');
  z=z-d-e; 
    z=z*2; 
      z=z+d+e; 
     z=z/(2*c); 
   z=z*(b-a)