% Tund nr. 6, ülesanne 3
% defineerime vektorfunktsiooni
% x=x(1), y=x(2)
F=@(x) [x(1)-x(2)^2+1,x(1)^2-x(2)];
% jooniselt määrame alglähendid
ezplot('x=y^2-1')
hold on
ezplot('y=x^2')
grid on
% alglähendid (-0.7,0.5) ja (1.2,1.5)
x=-0.7;
y=0.5;
G=[0,1/(2*sqrt(y));1/(2*sqrt(x+1)),0];
norm(G,1) % sobib, sest <1
x=1.2;
y=1.5;
G=[0,1/(2*sqrt(y));1/(2*sqrt(x+1)),0];
norm(G,1) % sobib, <1
x=[-0.7,0.5]; % alglähend reavektorina
i=0;
while norm(F(x),1)>10^(-5)
  x(1)=-sqrt(x(2));
  x(2)=sqrt(x(1)+1);
  i=i+1;
end
lahend1=x
sammud1=i

x=[1.2,1.5]; % alglähend reavektorina
i=0;
while norm(F(x),1)>10^(-5)
  x(1)=sqrt(x(2));
  x(2)=sqrt(x(1)+1);
  i=i+1;
end
lahend2=x
sammud2=i