% joonis, alglähendid
ezplot('x^2+3*y-1=0')
hold on
ezplot('2*x-5*y^2+9=0')
hold off
grid on
% alglähendid (-2,-1),(2.5,-1.7)
% uurime koondumise tingimust alglähendite korral
% a) alglähendi (-2,-1) korral
x=-2;
y=-1;
Gtuletis=[0 -3/(2*sqrt(1-3*y))
1/(5*sqrt(2*x+9)) 0];
norm(Gtuletis,1)
% b) alglähendi (2.5,-1.7) korral
x=2.5;
y=-1.7;
Gtuletis=[0 -3/(2*sqrt(1-3*y))
1/(5*sqrt(2*x+9)) 0];
norm(Gtuletis,1)
% lahendame hariliku iteratsioonimeetodiga
% defineerime F, x=z(1), y=z(2)
F=@(z) [z(1).^2+3*z(2)-1,2*z(1)+5*z(2).^+9];
% a) alglähend (-2,-1)
z=[-2,-1];
% vektorfunktsioon G
G=@(z) [-sqrt(2*z(1)+9),-sqrt(1-3*z(2))];
i=0;
while (norm(F(z),1))>10^-3
z=G(z);
i=i+1;
end
lahend1=z
sammud1=i
% b) alglähend (2.5,-1.7)
z=[2.5,-1.7];
G=@(z) [sqrt(2*z(1)+9),-sqrt(1-3*z(2))];
i=0;
while (norm(F(z),1))>10^-3
z=G(z);
i=i+1;
end
lahend2=z
sammud2=i