Function rgsl::elliptic::carlson::ellint_RJ_e

pub fn ellint_RJ_e(x: f64, y: f64, z: f64, p: f64, mode: Mode) -> (Value, Result)

This routine computes the incomplete elliptic integral RJ(x,y,z,p) to the accuracy specified by the mode variable mode.