pub fn single_slater(
    ao_points: &mut MatrixFull<f64>,
    c_mat: &mut MatrixFull<f64>,
    coeff: f64
) -> f64