func7p1.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 28 行

%func7p1.m二分法解方程
%(1)
s=1;f=3;tol=0.01;i=0
difference=1
while(abs(difference)>tol)
    if(sub7p11((s+f)/2)>0)
        f=(s+f)/2;
    else
        s=(s+f)/2;
    end
    difference=sub7p11((s+f)/2);
    i=i+1
end
%用四次循环得到解 1.9375，误差0.004
%(2)
s=0;f=1.5;tol=0.01;i=0
difference=1
while(abs(difference)>tol)
    if(sub7p12((s+f)/2)>0)
        f=(s+f)/2;
    else
        s=(s+f)/2;
    end
    difference=sub7p12((s+f)/2);
    i=i+1
end
%  用6次循环得到解 1.1133，误差0.0012