function p=newton(f_str,df_str,p0)
tol=0.00001;
f = inline(f_str);
df= inline(df_str);
while 1
　p=p0-f(p0)/df(p0);
　if abs(p-p0)<tol, break; end
　p0=p;
　disp(p0);
end
disp(p)