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.
    D.h.

    $\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)$

    Und dann kannst du RKF auf dieses System loslassen. Ich hoffe, ich habe das Problem verstanden ..


Anmelden zum Antworten