clear %vahemälu puhtaks

%graafikud alglähendite leidmiseks
ezplot('x-0.8*x^2+1.5*y',[-3,3,-1,4])
hold on
grid on
g=ezplot('x-1+y-sin(x-y)',[-3,3,-1,4])
set(g,'color','k') %teeme teise joone mustaks
legend('x-0.8*x^2+1.5*y','x-1+y-sin(x-y)')
title('Joonis alglähendite leidmiseks')
hold off

%jooniselt on näha 4 lõikepunkti. Loeme lõikepunktide
%koordinaadid alglähenditeks. Meile sobib
%(-2;3.5) ja (-1;1.1) ja (-0.6;0.6) ja (1.5;0.4)

%Teeme võrrandisüsteemi vasakutest pooltest vektori, mille 
%argumendiks on z=(z(1),z(2)), 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))];

%Lahendame süsteemi
lahend1=fsolve(F,[-2,3.5]) %alglähend on jooniselt
lahend2=fsolve(F,[-1,1.1])
lahend3=fsolve(F,[-0.6,0.6])
lahend4=fsolve(F,[1.5,0.4])