tom634

I changed math.sqrt to np.sqrt, but now I have another problem - quad divided to real and imaginary parts and quadpy calculating both parts give different results - why? Here...

Thank you for the comment. I have another issue with integrating integral. My code is: ```python from scipy.integrate import quad import numpy as np import scipy from scipy.special import jv...

When I change: `return quad(f_D, -0.5, 0.5, limit=50, args=(z))[0]` to `return quadpy.quad(f_D, -0.5, 0.5, limit=50, args=(z))[0]` I have error: ``` Exception has occurred: TypeError f_D() takes 2 positional arguments but...

I tried looping with an error: `Exception has occurred: TypeError f_D() argument after * must be an iterable, not numpy.float64` ``` def psi_D(z): print("psi_D",z) for e in z: print(e) return...

Full code: ``` from scipy.integrate import quad import numpy as np import scipy from scipy.special import jv import math import quadpy B=8 m=4 M_x = 0.3 M_T = 1 theta=0.3...

This is loop over z values in order to calculate internal function `f_D(x, z)` in order to calculate integral of `def psi_D(z):` . But currently I have no idea how...

I solved above problem by changing (z) to [z]: ``` from scipy.integrate import quad import numpy as np import scipy from scipy.special import jv import math import quadpy B=8 m=4...

In order to simplify the above question, I show the equations (simplified for implementation purposes) that I want to calculate numerically.  and the code: ``` from scipy.integrate import quad...

Hi @nschloe do you know maybe how to resolve this issue?

bump, I have the same problem