Matlab - Problema con codigo de Newton Rapson

function newton(f,x0,tol)

sym x;

df=diff(f,'x');

f=inline(f);

df=inline(char(df));

fprintf('\n it. x f(x) \n')

i=0;

fprintf('%3.0f %10.10f %10.10f \n',i,x0,f(x0))

x1=x0-(f(x0)/df(x0));

while (abs(x0-x1))/x0 >tol;

i=i+1;

fprintf('%3.0f %10.10f %10.10f \n',i,x1,f(x1))

x0=x1;

x1=x0-(f(x0)/df(x0));

end

fprintf('\n La aproximación de la raíz es: %3.10f \n\n',x1)

df=diff(f,'x');

>>syms x;   newton(1+sin((3.*x)/5)-((5.*x)./3),0,0.00000001)

 it. x f(x)

  0 0.0000000000 1.0000000000

  1 0.9375000000 -0.0291973265

  2 0.9123106074 -0.0000604211

  3 0.9122582645 -0.0000000003

 La aproximación de la raíz es: 0.9122582643