%defineerime funktsiooni
f=@(x,y)10*(y-x)*cos(y+x);

%lahendame DV numbriliselt
[x,y]=ode45(f,[1,5],3)

%joonis
plot(x,y,'r*')