clear %vahemälu

%joonis alglähendite jaoks
ezplot('x-0.8*x^2+1.5*y',[-3,3,-1,4])
hold on
h=ezplot('x-1+y-sin(x-y)',[-3,3,-1,4])
set(h,'color','r')
grid on
title('Joonis alglähendite leidmiseks')
xlabel('x')
ylabel('y')

%Jooniselt on näha, et süsteemil on 4 reaalset lahendit. Loeme
%lõikepunktide koordinaadid alglähenditeks.
%Seega (-2,3.5) ja (-1,1.2) ja (-0.7,0.6) ja (1.5,0.3)

%Teeme võrrandisüsteemi võrrandite vasakutest pooltest vektori, mille 
%argumendiks on ka vektor 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))];

lahend1=fsolve(F,[-2,3.5]) %alglähendiks on lõikepunkti koordinaadid
lahend2=fsolve(F,[-1,1.2])
lahend3=fsolve(F,[-0.7,0.6])
lahend4=fsolve(F,[1.5,0.3])