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]));
}