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

}