aynur

syms x;
eps = 0.000001;
x0 = input('enter x ');
f = input('enter your function ');
df = diff(f, x);
xn = x0-double(subs(f, x0))/double(subs(df, x0));
while(abs(xn - x0)>eps)
    x0 = xn;
    xn = x0-double(subs(f, x0))/double(subs(df, x0));
end