Revision history [back]

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=pq2^(24) b=sqrt(b) T1=2^50 print(T1) F = lambda x: 1/(sqrt(2pi)b)e^(-(x-a)^2/(2b*b)) I=quad(F, [-inf,T1]) I=round(I,1000)

print(I<2^(-15))

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=pq2^(24) b=sqrt(b) T1=2^50 print(T1) T=2^50

F = lambda x: 1/(sqrt(2pi)b)e^(-(x-a)^2/(2b*b)) I=quad(F, [-inf,T1]) [-inf,T]) I=round(I,1000)

print(I<2^(-15))print(I)

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)
a=p*2^(24)
 b=p*q*2^(24)
b=sqrt(b)
T=2^50
 
b=pq2^(24)
b=sqrt(b)
T=2^50
 
F = lambda x: 1/(sqrt(2pi)b)e^(-(x-a)^2/(2b*b))
1/(sqrt(2*pi)*b)*e^(-(x-a)^2/(2*b*b))
I=quad(F, [-inf,T])
I=round(I,1000)I=round(I,1000)
  
print(I)print(I)