Ha oui, en effet, je ne connaisait pas cette méthode...
http://www.aiaccess.net/French/Glossaires/GlosMod/f_gm_monte_carlo.htm
Bon, ba c'est pas très compliqué alors.
Il faut trouvé un générateur de nombre aléatoire (rand() ?) et c'est parti.
L'avantage du fortran c'est que les tableaux sont dimensionné automatiquement sur la dernière dimension.
REAL*8 som,V
c nombre de point max :
parameter Nmax 500
DO 1 N=1,Nmax
som=som+d()
1 continue
write(6,*)'resultat de :',som*V/Nmax
Ensuite, il faut définir une fonction d() donnant une valeur de façon aléatoire en fonction de ta fonction à intégrer f():
function d()
x1=rand()*(x1m-x1p)
x2=rand()*(x2m-x2p)
...
return f(x1,x2,...)
Voilà en gros le programme.
Salutation !
Char Snipeur