systemnewton.m

来自「Nonlinear Problem Files Collection」· M 代码 · 共 29 行

function [] = systemNewton(x0,maxIte)

format compact;

disp(sprintf('  k             x_k                     y_k'));
x=x0(1);
y=x0(2);

k=0;
disp(sprintf('%3d      %1.14f        %1.14f',k,x,y));

for k=1:maxIte

  Fx=sin(x)-y;
  Fy=x-cos(y);

  cx=cos(x);
  sy=sin(y);

  dy=(Fx-Fy*cx)/(1+sy*cx);
  dx=-Fy-sy*dy;

  x=x+dx;
  y=y+dy;
  disp(sprintf('%3d      %1.14f        %1.14f',k,x,y));

end