Indice generale

Runge Kutta di J.R.Cash, A.H.Karp codificato in Javascript

Una delle migliori e pi¨ moderne (si fa per dire) versioni di algoritmo di Runge Kutta con possibilitÓ di controllo della precisione e dunque con possibilitÓ di variazione automatica del passo di integrazione nella soluzione di un sistema di equazioni alle derivate ordinarie ( ODE ), Ŕ quello di Cash e Karp.
Qui ne propongo una versione a coefficienti interi ( il che elimina qualsiasi possibilitÓ di errore di troncamento nella specifica dell'algoritmo stesso ).
Il linguaggio di programmazione Ŕ Javascript ma Ŕ facile trasformare il sorgente in C o in Fortran.



L'algoritmo Ŕ codificato nel file: lib-cki.txt che Ŕ direttamente riutilizzabile...

Vedere inoltre: La struttura dell'algoritmo

Test elementare

Provo questo sistema:
y``+ π2·y = 0
che si splitta nel seguente sistema doppio :
v`= −π2·y
y`= v
ovvero
   y[1]` = y[2];
   y[2]` = -Math.PI*Math.PI*y[1];
   y[3]` = 1;
Con la condizione iniziale y[1]=1; y[2]=y[3]=0. Dato che il periodo vale 2 il passo di integrazione sia di 1/20 in modo che ogni 40 passi si debba riottenere la condizione iniziale.
Il valore di h deve essere dunque 1/(20*979292160.0).

Risultati...

...