private double v_prim(double t, double x1, double x2, double v)
{
double wartosc = 0.0;
double b = wspolczynniki.b;
double m = wspolczynniki.mass;
double k1 = wspolczynniki.k1;
double k2 = wspolczynniki.k2;
double l1 = wspolczynniki.l1;
double l2 = wspolczynniki.l2;
wartosc = (1 / m) * wspolczynniki.input_func(t) - (k2 / m) * x2 - (l2 / m) * Math.Pow(x2, 3) + (k2 / m) * x1
+ (l2 / m) * Math.Pow(x1, 3) - (k1 / m) * x2 - (l1 / m) * Math.Pow(x2, 3);
return(wartosc);
}