Newton Raphson en Matlab

%newton raphson
limpiar
fun=input('Ingresa la f(x) entre comillas: ');
f=inline(fun);
df=input('Ingresa la derivada de la función entre comillas: ');
df=inline(df);
x=input('Ingresa el x0: ');
tol=input('Ingresa la tolerancia: ')
i=1;
x(i)=x;
corr=tol+1;
k=0;
maxit=10;


while k<maxit & abs(corr)>tol
    fx=feval(f,x(i));
    dfx=feval(df,x(i));
    corr=fx/dfx;
    x(i+1)=x(i)-corr;
    i=i+1;
end


fprintf('iter   resultado \n')
for j=1:i
fprintf('%d  \t %10.5f \n',(j-1),x(j))
end


if abs(b-a)>tol
    error('Se alcanzó el número máximo de iteraciones')
end