simulation de Monte Carlo avec le module Numpy

import numpy as np
 
x1 = np.random.normal(0, stv, N)
y1 = np.random.normal(0, stv, N)
x2 = np.random.normal(100,stv, N)
y2 = np.random.normal(0, stv, N)
x3 = np.random.normal(0, stv, N)
y3 = np.random.normal(100, stv, N)
x4 = np.random.normal(100, stv, N)
y4 = np.random.normal(100, stv, N)
 
# calcul des surfaces 
surfaces = np.sqrt((x2-x1)**2+(y2-y1)**2)*np.sqrt((x3-x1)**2+(y3-y1)**2)/2+np.sqrt((x4-x2)**2+(y4-y2)**2)*np.sqrt((x4-x3)**2+(y4-y3)**2)/2
 
# résultats d'une des simulations
print mu, sigma
9999.1573234 199.631782423
# résultats en tenant compte de la rêgle des chiffres significatifs
print "Resultats :"
print "Moyenne = ","%.f" % mu
print "déviation standard  = ", "%.f" % sigma
 
Resultats :
Moyenne =  10000
déviation standard =  200