format long

x = [0; 1];

J = [1+0.03*cos(x(1)+x(2)), 2+0.03*cos(x(1)+x(2));5+0.07*sin(x(1)-x(2)), 6+0.07*sin(x(1)-x(2))];
F = [x(1)+2*x(2)+0.03*sin(x(1)+x(2))-4; 5*x(1)+6*x(2)+0.07*cos(x(1)-x(2))-8];

xx = x - inv(J)*F;

for i = 1:10^4
  x = xx;
  J = [1+0.03*cos(x(1)+x(2)), 2+0.03*cos(x(1)+x(2));5+0.07*sin(x(1)-x(2)), 6+0.07*sin(x(1)-x(2))];
  F = [x(1)+2*x(2)+0.03*sin(x(1)+x(2))-4; 5*x(1)+6*x(2)+0.07*cos(x(1)-x(2))-8];

  xx = x - inv(J)*F;
end
x