pub fn heat_kernel_sphere( t: f64, theta1: f64, phi1: f64, theta2: f64, phi2: f64, r: f64, max_l: usize, ) -> f64