# Butcher Tableau
const B = [[1/5],
     [3/40, 9/40],
     [44/45, -56/15, 32/9],
     [19372/6561, -25360/2187, 64448/6561, -212/729],
     [9017/3168, -355/33, 46732/5247, 49/176, -5103/18656],
     [35/384, 0, 500/1113, 125/192, -2187/6784, 11/84]
     ]

# Step size fractions
const nodes = [1/5, 3/10, 4/5, 8/9, 1, 1]

# Weights for 5th order method
const b5 = [5179/57600, 0, 7571/16695, 393/640, -92097/339200, 187/2100, 1/40]
# Weights for 4th order method
const b4 = [35/384, 0, 500/1113, 125/192, -2187/6784, 11/84, 0]
# Error estimate
const errest = b5 .- b4

#Interpolation coefficients
const interpC = hcat(
[1,       -183/64,      37/12,       -145/128],
[0,       0,             0,             0],
[0,        1500/371,     -1000/159,   1000/371],
[0,       -125/32,      125/12,      -375/64],
[0,       9477/3392,    -729/106,    25515/6784],
[0,       -11/7,        11/3,        -55/28],
[0,       3/2,          -4,           5/2],
)