Ode Problem

[rafal__]

Jak otrzymać prawidłową trajektorię dla poniższego ukłądu równań ? Od pewnego momentu ode

daje nieprawidłowe wyniki. Zmiana metody rozwiązania i kroku czasowego nie pomaga.

clear

t0=0;

tk=5.5;

dt=0.1;

t=t0:dt:tk;

a=7;

b=3;

c=5;

d=3;

 

function w=f(t,Y) // równania ruchu

w=zeros(2,1);

w(1)=(Y(1)+a)*(b-Y(2));

w(2)=(c-Y(2))*(Y(1)+d*Y(2));

endfunction

 

x=ode([-8.5;2],t0,t,f);

plot2d(x(1,: ),x(2,: ));