正态分布随机数生成器

import math,random
u=1800.
q=800.
reslist=[]
for v in xrange(0,3600):
   y=math.exp(-((v-u)/q)*((v-u)/q)/2)/(q*math.sqrt(2*3.14159265359))*36000
   for d in xrange(0,int(y)):
      reslist.append(v)
res2=[]
for v in xrange(0,36000):
   res2.append(reslist[random.randint(0,len(reslist))])
for r in res2:
   print r