secant.asv

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

function s=Secant(fx,x0,x1,n)
%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
for i=1:n
    x2=x1-f1*(x0-x1)/(f0-f1);
    x=x2;
    f2=eval(fx);
    x=[i,x0,x1,x2,f2];
    disp(x);
    x0=x1;
    
    x1=x2;
    i=i+1;
end