main.m

极有用的数据处理程序

clear all
format long

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

max1=100;

eps1=1.0e-8;
eps2=1.0e-8;
x0=[1.0 1.0  1.0]';          %第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