数值计算方法迭代法matlab实现

clear;
x0=1;
e=0.5e-3;%误差限
x1=phi(x0);
while(abs(x0-x1)>e)
    x0=x1;
    x1=phi(x0);
end
fprintf("%f\n",x0);
function x1=phi(x0)%迭代函数
x1=x0-(x0^8-175)/(8*x0^7);
end