# Integration using quad.

In my understanding, I is the integration of normal distribution from -infinity to 2^50. The value should be 1. But it is giving 0. I need very good precision say up to 50 decimal places.

```
p=0.578
q=1-p
a=p*2^(24)
b=p*q*2^(24)
b=sqrt(b)
T=2^50
F = lambda x: 1/(sqrt(2*pi)*b)*e^(-(x-a)^2/(2*b*b))
I=quad(F, [-inf,T])
I=round(I,1000)
print(I)
```

How is

`quad`

defined ? What is`inf`

? Are you using Sage ?I used " from mpmath import *" and "import time". It is running. My Sage version is 8.2

Suggestions: