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.