%defineerime funktsiooni f(x,y)
f=@(x,y)10*(y-x)*cos(y+x);
[x,y]=ode45(f,[1,5],3)
plot(x,y,'r*')