%Joonis alglähendite leidmiseks
ezplot('x-0.8*x^2+1.5*y=0',[-3,2,-1,4])
hold on
grid on
h=ezplot('x-1+y-sin(x-y)=0',[-3,2,-1,4])
set(h,'color','r')
legend('x-0.8*x^2+1.5y=0','x-1+y-sin(x-y)=0')
%Jooniselt näeme, et neil kahel joonel on 4 lõikepunkti. Lõikepunktide
%ligikaudsed koordinaadid sobivad meile alglähenditeks, seega
%valime (-2;3.5) ja (-1;1.2) ja (-0.7;0.7) ja (1.6;0.3)

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

%Lahendame süsteemi
lahend1=fsolve(F,[-2,3.5]) %alglähendiks vektor
lahend2=fsolve(F,[-1,1.2])
lahend3=fsolve(F,[-0.7,0.7])
lahend4=fsolve(F,[1.6,0.3])

kontroll1=F(lahend1)