Runge-Kutta 4 steps
You will find this code in the EJS console at page Initialization
lines marked // include remarks that do not influence the calculation.
xRK[0]=0;
yRK[0]=Anfangswert;
AbleitungRK[0]=yRK[0];//for exponential
for (int i=1; i<n; i++) { //loop definition)
y_a[i]=yRK[i-1]+delta/2*AbleitungRK[i-1];
Ableitung_a[i]=y_a[i];//for exponential
y_b[i]=yRK[i-1]+delta/2*Ableitung_a[i];
Ableitung_b[i]=y_b[i];//for exponential
y_c[i]=yRK[i-1]+delta*Ableitung_b[i];
Ableitung_c[i]=y_c[i];//for exponential
yRK[i]=yRK[i-1]+delta/6*(AbleitungRK[i-1]+2*Ableitung_a[i]+2*Ableitung_b[i]+Ableitung_c[i]);
xRK[i]=xRK[i-1]+delta;
AbleitungRK[i]=yRK[i];//for exponential
FRK[i]=(Anfangswert*Math.exp(xRK[i])-yRK[i])/(Anfangswert*Math.exp(xRK[i]));
//back to start of loop
}