function z=newtonraphson(f,x0,Tol)
syms x;
z=[x0];
error=1;
while error>Tol
dfx0=subs(diff(f(x)),x0);
x1=x0-f(x0)/dfx0;
error=abs(x1-x0);
z=[z;x1];
x0=x1;
end