Matlab 数值计算----迭代法计算非线性方程组在指定区间的根

 iterate.m函数

function [x_star,index,it]=iterate(phi,x,ep,it_max)
if nargin<4 it_max=100;end
if nargin<3 ep=1e-5;end
index=0;k=1;
while k<it_max
x1=x;x=feval(phi,x);
if abs(x-x1)<ep
index=1;break;
end
k=k+1;
end
x_star=x;it=k;

test.m函数（测试二分法求非线性方程在指定区间解的函数）
phi=inline('(x+1)^(1/3)');
[x_star,index,it]=iterate(phi,1.5)

结果：
>> test

x_star =

    1.3247

index =

     1

it =

     7