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