Problems with the calculate_distance_to_harmonic_approx method
I have problems with the method calculate_distance_to_harmonic_approx of the BFieldAnalyzer
class. This method produces differing output compared to Mathematica when using NIntegrate to calculate same integrand (The Mathematica file used for this calculation is attached to this issue: magnetic_fields_distance_examples.nb). I already checked the integrand, it produces the same output in the Mathematica and the Python file. Additionally, if I compare the distance of the Z-trap to its harmonic approximation in Python and Mathematica, I get the same results. I suspect that it has to do something with the quadrature which is used by quadpy.