gx.m

来自「有关matlab r2007的实用教程」· M 代码 · 共 10 行

function x2=GX(x0,x1)
x2=x1-ff(x1)*(x1-x0)/(ff(x1)-ff(x0));   %第一次割线选代计算
n=1;
while (abs(x1-x0)>=1e-4)&(n<=100000)
    x0=x1;
    x1=x2;
    x2=x1-ff(x1)*(x1-x0)/(ff(x1)-ff(x0));    %割线迭代求解
    n=n+1;
end