Hi all,I'm trying to integrate stats.norm.pdf in the range [-np.inf, x] using integrate.quad. I'm having a few issues where the integral dives to zero when it should be one - https://gist.github.com/andyfaff/c8a94ab87333c42791e58e7d6e12a05eIt does this when norm.pdf(x) reaches the smallest number expressible by a float, and for some unknown reason when x ~ 20.Is there anyway I can mitigate against this behaviour using `integrate.quad` options? I'm hoping this is expected behaviour.
cheers,Andrew.p.s. I'm using integration instead of norm.cdf for a specific reason.-- _____________________________________Dr. Andrew Nelson_____________________________________
SciPy-User mailing list