pub fn weak_error(exact_mean: f64, numerical_mean: f64) -> f64
Weak error: |E[exact(T)] - E[numerical(T)]|.