## Exponentielles Wachstum

Aufstellen der Differentialgleichung dN/dt=k*N

dgl:'diff(N,t)=k*N;

Lösung der Differentialgleichung über den Befehl ode2

ode2(dgl,N,t);

Runge-Kutta-Verfahren - Numerische Näherung

Zuerst wird das Zusatzpaket "dynamics" geladen, über das der 
Befehl rk (Runge-Kutta 4.Ordnung) zur Verfügung gestellt wird.
Das Ergebnis ist eine Liste von x-y-Wertepaaren.
Dieses wird über plot2d ausgegeben.

load(dynamics);

loes:rk(0.2*N,N,100,[t,0,15,0.1])$

wxplot2d([discrete,loes]);

Darstellung im Richtungsfeld

Das Zusatzpaket "plotdf" (Plot Direction Field) wird geladen.
Anschließend wird für die rechte Seite der Differentialgleichung
das Richtungsfeld ausgegeben (Achtung: xMaxima muss installiert sein!

x und y-Bereich werden festgelegt, ein Schieberegler für den Parameter k eingeführt

load(plotdf);

plotdf(k*y,[x,-2,16],[y,-100,800],
  [parameters, "k=0.1"],[sliders,"k=0.01:1"],[tstep,0.1]);