Function quantile 

pub fn quantile<T, D>(
    a: &Array<T, D>,
    q: T,
    axis: Option<usize>,
) -> Result<Array<T, IxDyn>, FerrayError>
where
    T: Element + Float,
    D: Dimension,

Compute the q-th quantile of array data along a given axis.

q must be in [0, 1]. Uses linear interpolation (NumPy default method). Equivalent to numpy.quantile.