import numpy as np


#antall ledd

for n in range(20):
    p = 1
    for k in range(1,n+1):
        p += ((-1)**k)/(2**k*np.math.factorial(k)*(2*k+1))
    p = p/np.sqrt(2*np.pi)
    print(p-0.3413447460685429485852)

