main.m

来自「极有用的数据处理程序」· M 代码 · 共 29 行

clear all
format long

%%%%%%%%%%%%% Newton-Raphson method %%%%%%%%%%%

max1=100;

eps1=1.0e-8;
eps2=1.0e-8;
x0=[-1 2]';          %第1题的初始迭代值

x=x0;

for i=1:max1,
    A=Fd(x);
    b=-F(x);
    dx=A\b;
    x=x+dx;
    sprintf('第%d次迭代结果：\n',i)
    [x  dx  F(x)]      %输出每次的x(i),dx(i),F(x(i))
    if (abs(dx)<eps1)&(abs(F(x))<eps2)
        break
     end
  end