/* Variablendeklaration */
var ODE; // ODE-Solver RK4
...
function setup ()
{
...
ODE = new RK_secondOrder();
}
function draw() // Hauptprogramm
{
...
//****************************************** Berechnung ***********************
if (START)
{
...
y1 = l0; // Startort = Ruhefederlänge
vy1 = 0; // Startgeschwindigkeit
}
else
{
...
result = ODE.RKsO(f, y1, vy1, dt); // Runge-Kutta
y1 = result[0];
vy1 = result[1];
t = t + dt;
}
//************************************** Gradientenfunktion ***********************
function f(y, vy)
{
return(g - 2*d*vy - sq(omega0)*(y-l0)); // g = 9,81; d: Dämpfung; l0: Ruhefederlänge
}