CODEWITHSUNDEEP
pythoncontinuous-integrationintegration
Felipe
Felipe

Reputation: 99

Numerical integration with big errors

I'm trying to compute this integral, but the errors seems too high.

I tried changing the precision with mp.mp.dps = 100 and maxdegree= 3000 but without any good results.

Is it possible to compute such an integral ?

import mpmath as mp
A = 106.829
B = 2546.40
DD = 1000.0
mp.mp.dps = 100

nom = lambda oo, mm: (2 * mm - 1) * mm**( DD * mm) * (
              mp.exp( (2 * mm - 1) * A -  (mm/oo) * B)) /  (mp.gamma(mm) ** DD * oo ** (DD * mm) )

nomRes = mp.quad(nom,[0.5,3.0],[0.5,3.0], error = True, maxdegree= 3000)

I got the following error :

nomRes = (5.2385053e-792,5.1842e-792)

This error suggests that I can't trust the value obtained.

Upvotes: 0

Views: 46

Answers (1)

Igor Rivin
Igor Rivin

Reputation: 4864

You are raising things to power DD*mm which is huge. What do you expect to get?

Upvotes: 1

Related Questions