Euler

der aktive Code befindet sich im Modell auf den Seite Initialization

mit // gekennzeichnete Eintragungen sind Bemerkungen, die den Rechenablauf nicht beeinflussen

xE[0]=0;

yE[0]=Anfangswert;

AbleitungE[0]=yE[0];//für e- Funktion

for (int i=1; i<n; i++) {//Schleifeninstruktionen)

xE[i]=xE[i-1]+delta;

yE[i]=yE[i-1]+delta*AbleitungE[i-1];

AbleitungE[i]=yE[i];//für e- Funktion

//zurück an Schleifenanfang

FE[i]=(Anfangswert*Math.exp(xE[i])-yE[i])/(Anfangswert*Math.exp(xE[i]));

}

Heun

xH[0]=0;

yH[0]=Anfangswert;

AbleitungH_a[0]=yH[0];//für e- Funktion

for (int i=1; i<n; i++) {//Schleifeninstruktionen)

xH[i]=xH[i-1]+delta;

yH_a[i]=yH[i-1]+delta*AbleitungH_a[i-1];

AbleitungH_b[i]=yH_a[i];//für e- Funktion

yH[i]=yH[i-1]+delta/2*(AbleitungH_a[i-1]+AbleitungH_b[i]);

AbleitungH_a[i]=yH[i];//für e- Funktion

//zurück an Schleifenanfang

FH[i]=(Anfangswert*Math.exp(xH[i])-yH[i])/(Anfangswert*Math.exp(xH[i]));

}