func7p2.asv

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

s=10;f=8;tol=0.01;
difference=1
while(abs(difference)>tol)
    if(sub7p2((s+f)/2)>0)
        f=(s+f)/2;
    else
        s=(s+f)/2;
    end
    difference=sub7p2((s+f)/2);
end
disp(difference);disp(s/2+f/2)
x=linspace(-6,0.2,100);y=linspace(-6,0.2,100);y=linspace(-6,0.2,100);xx=linspace(-6,0.2,100)
for k=1:100
    yy=1/y
plot(sin(x(k)),x(k),'-',1/y(k),y(k),'.-')
end
%从图中易看出其解区间,以二分法算得解为:1.1094,2.7813,6.4375,9.3125,如解区间出错,则无法算得结果
%作图[-6,0.2]有解区间