regulafalsi.m

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

function out=RegulaFalsi(a,b)
fprintf(' \t\t Regula Falsi Method \n')
fprintf(' \tsteps\ta\t\tb\t\tmid\t\tf(a)\t\tf(b)\t\tf(mid) \n')
for i=1:10
if f(a)*f(b)<0
    mid=b-((f(b)*(b-a))/(f(b)-f(a)));
    if (f(a)*f(mid))<0
        b=mid;
    else
        a=mid;
    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,a,b,mid,f(a),f(b),f(mid))
end %for

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