[C] Runge Kutta del quarto ordine

[Dracmaleontes]

Si capisco cosa intendi, purtroppo non sono ancora in grado di realizzare programmi "portatili" e ho sempre da modificare qualche cosa ad ogni piccola variazione. Comunque ho provato a stimare l'errore sull'energia al variare di dt e mi viene una curva del tipo $f(x) = ax^5 $ o giù di li. Il prof mi ricordo disse che per l'oscillatore e per piccoli dt era un metodo dell'ordine 5 quindi credo ci posso stare. L'oscillatore smorzato e con forzante l'ho provato già con RK II. però come dici tu ho dovuto passare quei 10 minuti a sistemare delle cosette. Comunque grazie dei consigli, proverò a migliorare nei prossimi programmi.

[vict85]

Portabilità si riferisce al trasferire il codice su un ambiente di esercuzione diverso (per esempio un diverso sistema operativo). Penso che tu intedessi che non sai realizzare programmi sufficientemente generici da poter essere riutilizzato per problemi diversi. Comunque anche quella è una questione di esercizio.

[Dracmaleontes]

Si hai ragione, ho fatto confusione.

[feddy]

Aggiungerei solamente che con RK4, nel caso del pendolo semplice (cioè senza forzante), non avrai conservazione dell'energia anche a livello della soluzione discreta. Ossia, se plotti in un grafico il tempo e in ordinata l'energia del sistema $E(t) = K+V = \frac{1}{2}m v(t)^2 + V(t) $ , osserverai un "gradino"!