func7p3.m

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

%func7p3.m求解非线性方程
%(1)迭代公式，请见代码
%(2)/(3)  x0=-(8000)^(1/2),求解。运用几何知识，可知在x(i)点的坐标为x(i)/2
x0=-(8000)^(1/2)+30;x00=-(8000)^(1/2);tol=0.001;x1=0;i=1
xx=linspace(x00,0,5)
while(x00<0)
    x0=x00+1
    while(abs(x1-x0)>tol)
      x1=x0
      x0=x0-(quad('sub7p3',x00,x0,0.0001)-200)/sub7p3(x0)
      i=i+1
    end
    x00=x0
    disp(x00/2)
end
%迭代5次，得解x(1)=-88.3176，误差-1.4102e-007,切线与轴所成交点x坐标为-44.1588，以下各点类同
%在n=114时，x(n)>0，故停止，x(114)=-1.0509,切线与轴所成交点x坐标为-0.5254，其余各点请执行本文件，在command
%window将有显示