Dear Sir or Madam, I don“t know to which address I should send this information, but I believe there is an error in GSL's documentation.
In the Ordinary Differential Equations chapter, it is written that one of the methods for numerical integration is the explicit embedded Runge-Kutta Prince-Dormand (8, 9) method (link <https://www.gnu.org/software/gsl/doc/html/ode-initval.html#c.gsl_odeiv2_step_type.gsl_odeiv2_step_rk8pd>). However, the highest-order method by these authors is the (7, 8), which is presented at "Prince, Peter J., and John R. Dormand. High order embedded Runge-Kutta formulae. Journal of computational and applied mathematics 7.1 (1981): 67-75.". I took a look at the rk8pd.c file and the method is indeed the (7, 8). Moreover, the above reference is also given as a comment in said file. I tried to contact the author, Gerard Jungman, but I got no reply. Kind regards, Vitor