secant.m

来自「Some methods Applied in matlab」· M 代码 · 共 22 行

function out=Secant(x0,x1)
fprintf(' \t\t Secant Method  \n')
fprintf(' \tsteps\tx0\t\tx1\t\tx2\t\tf(x0)\t\tf(x1)\t\tf(x2) \n')
for i=1:5
if f(x0)*f(x1)<0
    x2=x1-((x1-x0)/(f(x1)-f(x0)));
    if (f(x0)*f(x2))<0
        x1=x2;
    else
        x0=x2;
    end
else 
    disp('Choose diffrent values for a and b');
end
%fprintf(8d,8.4f,8.4f,8.4f,8.4f,8.4f,8.4f,)
fprintf('%8d  %8.4f%8.4f%8.4f  %8.4f  %8.4f    %8.4f \n', i,x0,x1,x2,f(x0),f(x1),f(x2))
end %for

end %function
function ret=f(x)
ret=x.^2-3;
end