% Tund nr. 6, ülesanne 2
% defineerime vektorfunktsiooni
% z=(x,y), tähistame z=x(1),z=x(2)
% F(x,y)=(x-0,8x^2+1,5y,x-1+y-sin(x-y))=0
% lühemalt
% F(z)=0
F=@(x,y) [x-0.8*x^2+1.5*y,x-1+y-sin(x-y)];
Fuus=@(x) [x(1)-0.8*x(1)^2+1.5*x(2),x(1)-1+x(2)-sin(x(1)-x(2))];
% uurime jooniselt, millised võiksid olla võrrandi-
% süsteemi lahendid, määrame alglähendid
ezplot('x-0.8*x^2+1.5*y=0',[-3,2,-2,4])
hold on
ezplot('x-1+y-sin(x-y)=0',[-3,2,-2,4])
hold off
grid on
% alglähendid (-2,3.6), (-1,1.2),(-0.6,0.6),(1.7,0.2)
lahend1=fsolve(Fuus,[-2,3.6])
lahend2=fsolve(Fuus,[-1,1.2])
lahend3=fsolve(Fuus,[-0.6,0.6])
lahend4=fsolve(Fuus,[1.7,0.2])