%program na riešenie nelin. DR 2.rádu Van-der-Pol oscilátor – zadaný v stavovom priestore
x0=[0  0.25]' ;     %inicializácia poè. Podmienok
t0=0;       tf=20;  %definícia èasového intervalu
[t,x]=ode23('vander',[t0, tf],[0; 0.25]);
[t,x]=ode45(@vander,[t0, tf],[0; 0.25]);
Subplot(211),plot(t,x(:,1))
Title('riešenie y(t)'),  xlabel('t'), grid
Subplot(212), plot(t,x(:,2))
Title('prva derivacia y'), xlabel('t'), grid
plot(t,x)  %dva grafy v jednom obr
Title('Van-der-Pol rovnica – èasova história');
Pause
