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