import numpy as np
import time

st = time.time()
x = np.zeros(7)
x[0] = -1
x[1] = -np.sqrt((5+2*np.sqrt(5/3))/11)
x[2] = -np.sqrt((5-2*np.sqrt(5/3))/11)
x[3] =  0
x[4] =  np.sqrt((5-2*np.sqrt(5/3))/11)
x[5] =  np.sqrt((5+2*np.sqrt(5/3))/11)
x[6] =  1 

w = np.zeros(7)
w[0] = 1/21
w[1] = (124-7*np.sqrt(15))/350
w[2] = (124+7*np.sqrt(15))/350
w[3] = 256/525
w[4] = (124+7*np.sqrt(15))/350
w[5] = (124-7*np.sqrt(15))/350
w[6] = 1/21

p = np.dot(w,np.exp(-x**2/2))/np.sqrt(2*np.pi)

et = time.time()

print(.5*p-0.3413447460685429485852)
print(st-st)