secant.m

来自「Numerical analysis （Gauss、Bisection、jaco」· M 代码 · 共 32 行

function s=Secant(fx,x0,x1,eps)
if nargin==3
    eps=1.0e-6;
elseif nargin<3
    return
end
x=x0;
f0=eval(fx);
x=x1;
f1=eval(fx);
if abs(f0)<abs(f1)
    temp1=x0;
    temp2=f0;
    x0=x1;
    f0=f1;
    x1=temp1;
    f1=temp2;
end
i=1;
f2=1;
while abs(f2)>eps
    x2=x1-f1*(x0-x1)/(f0-f1);
    x=x2;
    f2=eval(fx);
    x=[i,x0,x1,x2,f2];
    disp(x);
    x0=x1;
    f0=f1;
    x1=x2;
    f1=f2;
    i=i+1;
end