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,: ));