Matlab - Saben por que me marca errores? Quiero programar la regla del trapecio

clear

input('Programa de interpolacion por el metodo del Trapecio');

f=input('Ingrese la funcion a evaluar');

a=input('Ingrese el intervalo superior');

b=input('Ingrese el intervalo inferior');

e=input('Ingrese el error minimo');

n=2;

gx=str2func(strcat('@(x) ',f));

h=(a-b)/n;

sum=0;

for i=1:h:n

    eval=a+h*i;

    sum=sum+gx(eval);

end

x=h/2*(gx(a)+2*sum+gx(b));

rTrapecio(f,a,b,e,x,n)

function y=rTrapecio(gx,a,b,e,x,n)

n=n+1;

f=str2func(strcat('@(x) ',gx));

h=(a-b)/n;

sum=0;

for i=1:h:n

    eval=a+h*i;

    sum=sum+f(eval);

end

y=(h/2)*f(a)+2*sum+f(b);

if norm((y-x),inf)>e

       rTrapecio(f,a,b,e,y,n)

end

end