Function integrate_kernel_sphere 

pub fn integrate_kernel_sphere(
    h: f64,
    kernel: SphKernel,
    n_samples: usize,
) -> f64

Numerically integrate the kernel over a sphere of radius h using Monte Carlo.

Useful for verifying kernel normalization (should integrate to ~1 in 3D).