func7p1.asv

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

%(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