% Harjutustund 6, ülesanne 3
% VÕRRANDISÜSTEEMI LAHENDAMINE MATLABI KÄSUGA
% vektorfunktsiooni defineerimine
% F(z)=(f1(z),f2(z)), kus z=(x,y), arvestame x=z(1) ja y=z(2)
F=@(z) [z(1)-0.8*z(1)^2+1.5*z(2),z(1)-1+z(2)-sin(z(1)-z(2))];
% võrrandisüsteemile vastavad jooned joonisele
ezplot('x-0.8*x^2+1.5*y=0',[-3,3,-1,4])
hold on
ezplot('x-1+y-sin(x-y)=0',[-3,3,-1,4])
hold off
yline(0)
grid on 
% alglähend(id) jooniselt (-2,3.6), (-1,1.2), (-0.6,0.6) ja (1.6,0.3)

% lahendame võrrandisüsteemi
lahend1=fsolve(F,[-2,3.6])
kontroll1=F(lahend1)
lahend2=fsolve(F,[-1,1.2])
lahend3=fsolve(F,[-0.6,0.6])
lahend4=fsolve(F,[1.6,0.3])