# calcul des surfaces
for n in range(N):
# évaluation de chaque sommet
x1 = gauss(0, stv)
y1 = gauss(0, stv)
x2 = gauss(100, stv)
y2 = gauss(0, stv)
x3 = gauss(0, stv)
y3 = gauss(100, stv)
x4 = gauss(100, stv)
y4 = gauss(100, stv)
x1 = gauss(0, stv)
y1 = gauss(0, stv)
x2 = gauss(100, stv)
y2 = gauss(0, stv)
x3 = gauss(0, stv)
y3 = gauss(100, stv)
x4 = gauss(100, stv)
y4 = gauss(100, stv)
# calcul de surface adapté
surface = sqrt((x2-x1)**2+(y2-y1)**2)*sqrt((x3-x1)**2+(y3-y1)**2)/2+sqrt((x4-x2)**2+(y4-y2)**2)*sqrt((x4-x3)**2+(y4-y3)**2)/2
surf=surf+[surface]
mu = pylab.mean(surfaces)
sigma = pylab.std(surfaces)
# résultats bruts d'une des simulations
print mu, sigma
10003.1814816 200.382615819
# 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 = 10003
déviation standard = 200