Function rgsl::elliptic::carlson::ellint_RD

pub fn ellint_RD(x: f64, y: f64, z: f64, mode: Mode) -> f64

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