matlab牛顿迭代法代码

clc
clear
syms x
f=input('请输入函数f(x):');%x^3/3-x
x0=input('请输入迭代初值:');%初值大小
fai=x-f/diff(f,x);  %求牛顿迭代函数
k=0 ; %迭代次数
err=1.000;
phi=input('请输入容许误差:');%允许误差大小
disp('迭代次数      求得根           误差');
while(err>phi)
    x1=subs(fai,x0);
    err=abs(x1-x0);
    x0=x1;
    k=k+1;
 end
fprintf('%d  %1.10f   %1.10f\n',k,x1,err);