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