A1 Matlab

Van der Pol - ova jednaĉina:
�x
��
( x
x
x�
x�
1
1
2
�
�
�
2
x�
,
�x
x
1
�1)
x�
�
1
x
2
( x
�
2
2
x
x
�1)
�
Primer upotrebe ode23
�
x
0
2
function xprim = vdpol( t, x )
% Opis Van der Pol-ove dif. jednacine
xprim = [ x(1) * (1-x(2)^2) - x(2)
x(1) ];
» tp = 0; tk = 20; x0 = [ 0 0.25 ];
» [ t, x ] = ode23( @vdpol, [tp tk], x0 );
» plot( t, x )
Napomena: Brojevi “23”, odnosno “45” u imenu funkcija ode23 i ode45 se odnose
na stepene Taylor-ovih redova upotrebljanih u odgovarajućim algoritmima, a
ne odnose se na red diferencijalne jednaĉine koja se rešava. U principu, metod
ode45 je namenjen rešavanju diferencijalnih jednaĉina visokog reda.