pub fn ellint_RJ(x: f64, y: f64, z: f64, p: f64, mode: Mode) -> f64
This routine computes the incomplete elliptic integral RJ(x,y,z,p) to the accuracy specified by the mode variable mode.