Ze względu na szybkość współczesnych procesorów, rozwiązywanie setek równań różniczkowych zwyczajnych w czasie rzeczywistym nie stanowi żadnego problemu. Dołączając do tego progranowanie dynamiczne (funkcje pamietające) możemy analizować procesy opisane skompilkowanymi układami r. różniczkowych zdyczajnych (ODE) z dziecinną łatwością. Dla przykładu podaję niżej przykład funkcji pamiętającej numeryczne rozwiązania r. Lotki-Volterry dla dowolnych parametrów i warunków początkowych:

In[1]:=

NDSolve_zaawansowane_1.gif

Funkcja ta generuje i zapamiętuje całe funkcje interpolujące, co jest najszybszą metodą

In[18]:=

NDSolve_zaawansowane_2.gif

Out[18]=

NDSolve_zaawansowane_3.gif

Sprawdzamy czy funkcja generuje liczby, w przeciwnym wypadku nie będzie mozna jej narysować:

In[19]:=

NDSolve_zaawansowane_4.gif

Out[19]=

NDSolve_zaawansowane_5.gif

Możemy narysować rozwiązanie dla dowolnej kombinacji sześciu parametrów:

In[20]:=

NDSolve_zaawansowane_6.gif

Out[20]=

NDSolve_zaawansowane_7.gif

A także manipulować parametrami, np. warunkami początkowymi, obserwując jak odejście od stanu równowagi (linie przerywane) powoduje oscylacje:

In[22]:=

NDSolve_zaawansowane_8.gif

Out[22]=

NDSolve_zaawansowane_9.gif

In[28]:=

NDSolve_zaawansowane_10.gif

Out[28]=

NDSolve_zaawansowane_11.gif

In[29]:=

NDSolve_zaawansowane_12.gif

Out[29]=

NDSolve_zaawansowane_13.gif

Spikey Created with Wolfram Mathematica 7.0