Runge-Kutta-Fehlberg
-
Kann mir jemand mal kurz erläutern, wie man mit RKF an folgendes DGL-System rangeht?
rvek.(t) = vvek(t)
vvek(t) = -w * (mp(t) / m(t)) * evek(t) + gvek(rvek(t)) / m(t)
m(t) = m0 - INT_0_t { mp(theta) dtheta <=> m.(t) = -mp(t)
vek = vektorielle Größe
. = Ableitung
INT_0_t = Integral 0 bis t
-
Du kannst die DGLen schreiben als (n+1)-dim. System von DGLen 1. Ordnung.
$\left( \begin{array}{c}\dot{r\_1} \\ \vdots \\ \dot{r\_n} \\ \dot{m} \end{array} \right) = \left( \begin{array}{c} f\_1 \\ \vdots \\ f\_n \\ -mp \end{array} \right)$
D.h.Und dann kannst du RKF auf dieses System loslassen. Ich hoffe, ich habe das Problem verstanden ..