Para el metodo de RK-Clasico (pagina 6) del Numerical AlgoritmsWhitC
se tienen los coeficientes

1   1/6,   0,
2   1/3,   1/2,     1/2
3   1/3,   1/2,      0         1/2
4   1/6,   1,        0          0         1


Se traduce como

k1 = h * f(t, w)
k2 = h * f(t + h/2, w + k1/2)
k3 = h * f(t + h/2, w + k2/2)
k4 = h * f(t + h, w + k3)

w = w + k1/6 + k2/3 + k3/3 + k4/6

La segunda columna sirve para multiplicar a las Ks y formar la nueva w
La tercera columna sirve multiplicar a h y sumarce al tiempo
La cuarta y las demas columnas sirven para multiplicar a las Ks y sumarlas a la w de cada Ks



Para RK-Fehlberg de 5 orden pagina 7, se tiene los coeficientes

1   16/135,         0,
2   0,              1/4,          1/4
3   6656/12825,     3/8,          3/32              9/32
4   28561/56430,    12/13,     1932/2197        -7200/2197     7296/2197
5   -9/50,          1,          439/216            -8          3680/519      -845/4104
6    2/55,          1/2,         -8/27              2         -3544/2565     1859/4104    -11/40


Se traduce como

k1 = h * f(t, w)
k2 = h * f(t + h/4, w + k1/4)
k3 = h * f(t + 3*h/8, w + 3*k1/32 + 9*k2/32)
k4 = h * f(t + 12*h/13, w + 1932*k1/2197 -7200*k2/2197 + 7296*k3/2197)
k5 ...
k6 ...

w = w + 16*k1/135 + 6656*k3/12825 + 28561*k4/56430 -9*k5/50 + 2*k6/55


De esta forma se implementa cualquiera de los metodos del libro