adding = 0
n = 10000
for j in range(n):
    x = random()*3
    y = random()*2-1
    adding += x^2*y^2
avg = adding/n
integralv = avg*6
show('n=',n,' gives integral approx = ',integralv)