Integration using quad.

savecancel

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)