f=@(x) x(1)+2*x(2)+0.03*sin(x(1)+x(2))-4;
g=@(x) 5*x(1)+6*x(2)+0.07*cos(x(1)-x(2))-8;
F=@(x) [f(x) g(x)]';
p0=[-1.995 2.985]';format long;
for i = 1:100
  p=F(p0);
  if norm(p-p0,inf) < 10^(-5)
    break;
    endif
    p0=p;
    endfor
i
p