Skip to Content

simulation correcte

  1. # calcul des surfaces
  2. for n in range(N):
  3. # évaluation de chaque sommet
  4. x1 = gauss(0, stv)
  5. y1 = gauss(0, stv)
  6. x2 = gauss(100, stv)
  7. y2 = gauss(0, stv)
  8. x3 = gauss(0, stv)
  9. y3 = gauss(100, stv)
  10. x4 = gauss(100, stv)
  11. y4 = gauss(100, stv)
  12. x1 = gauss(0, stv)
  13. y1 = gauss(0, stv)
  14. x2 = gauss(100, stv)
  15. y2 = gauss(0, stv)
  16. x3 = gauss(0, stv)
  17. y3 = gauss(100, stv)
  18. x4 = gauss(100, stv)
  19. y4 = gauss(100, stv)
  20. # calcul de surface adapté
  21. 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
  22. surf=surf+[surface]
  23.  
  24. mu = pylab.mean(surfaces)
  25. sigma = pylab.std(surfaces)
  26. # résultats bruts d'une des simulations
  27. print mu, sigma
  28. 10003.1814816 200.382615819
  29. # résultats en tenant compte de la rêgle des chiffres significatifs
  30. print "Resultats :"
  31. print "Moyenne = ","%.f" % mu
  32. print "déviation standard = ", "%.f" % sigma
  33.  
  34. Resultats :
  35. Moyenne = 10003
  36. déviation standard = 200