%võrrandisüsteemi lahendamine Matlabi abil

%alglähendi leiame graafiliselt
ezplot('x-0.8*x^2+1.5*y=0',[-3,3,-1,4])
hold on
g=ezplot('x-1+y-sin(x-y)=0',[-3,3,-1,4])
set(g,'color','r') %teeme teise joone punaseks
grid on
hold off
title('Joonis alglähendite leidmiseks')

%Jooniselt on näha, et joontel on 4 lõikepunkti. Loeme lõikepunktide
%koordinaadid alglähenditeks. Meile sobivad (-2,3.5) ja (-1,1.1) ja
%(-0.6,0.6) ja (1.5,0.3).
%Defineerime funktsioonide vektori, kus esimene võrrand on f1 ja teine
%võrrand f2 ja argumendiks on z, kusjuures 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])
lahend2=fsolve(F,[-1,1.1])
lahend3=fsolve(F,[-0.6,0.6])
lahend4=fsolve(F,[1.5,0.3])

kontroll1=F(lahend1)