Yes. You are right. I now use scipy.integrate.quad and have validated the integral values.