func7p3-fail.m

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

%(1)迭代公式，请见代码
%(2)x0=-(8000)^(1/2),求解
x0=-(8000)^(1/2)+30;x00=-(8000)^(1/2);tol=0.001;x1=0;i=1
while(abs(x1-x0)>tol)|(i<10)
    x1=x0
    x0=x0-(quad('sub7p3',x00,x0,0.0001)-200)/sub7p3(x0)
    disp(abs(x1-x0))
    disp(quad('sub7p3',x0,x00)-200)
    i=i+1
end