from math import * from random import * N=1000000 alpha=1.0 numerator=0.0 denominator=0.0 for i in range(0,N): ## insert code for generating sample point here numerator+=x*x*exp(-alpha*x*x) ## insert code for estimating denominator here numerator/=N denominator/=N print "Monte Carlo mean ", numerator/denominator