%Alglähendite otsimine graafiliselt
ezplot('x-0.8*x^2+1.5*y=0',[-3,3,-1,4])
hold on
grid on
h=ezplot('x-1+y-sin(x-y)=0',[-3,3,-1,4])
set(h,'color','r') %teeme teise joone punaseks
title('Joonis alglähendite leidmiseks')
legend('x-0.8*x^2+1.5*y=0','x-1+y-sin(x-y)=0')
%Jooniselt näeme, et joontel on 4 lõikepunkti => süsteemil on 4 lahendit.
%Valime alglähenditeks (-2,3.6), (-1,1.1), (-0.5,0.6) ja (1.5,0.2).

%Defineerime süsteemi vasakud pooled vektorfunktsioonina, mille argument on
%vektor z, kus z(1)=x ja z(2)=y.
F=@(z)[z(1)-0.8*z(1)^2+1.5*z(2),z(1)-1+z(2)-sin(z(1)-z(2))];

%Lahendamine
vastus1=fsolve(F,[-2,3.6]) %alglähend on jooniselt, vektorina esitatud
vastus2=fsolve(F,[-1,1.1])
vastus3=fsolve(F,[-0.5,0.6])
vastus4=fsolve(F,[1.5,0.2])

kontroll1=F(vastus1)