Méthode de Monte Carlo

La méthode de monte carlo ne sert pas à calculer directement une intégrale.
Je te conseil de te renseigner là dessus. Cette méthode sert à retrouver certains phénomènes physique par le biais de tirages aléatoire. Je n'en dit pas plus, j'ai peur de dire une connerie.

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.

program Monte
implicit none
double precision x,Rf,IMf,RI,IMI,Rand
integer N,I
complex II,JJ
logical Domaine
open(10,file='x.dat')
JJ=(0d0,1d0)
RI=0d0
IMI=0d0
N=0

do 1 I=1,100
x=20d0
Domaine=(x.le.10d0).and.(x.ge.2d0)
do while(.not.Domaine)
x=8d0*Rand()+2
write(10,*)x
write(*,*)x
Domaine=(x.le.10d0).and.(x.ge.2d0)
enddo
RI=RI+Rf(x)
IMI=IMI+IMf(x)
N=N+1
1 Continue
RI=8d0*RI/N
IMI=8d0*IMI/N
write(*,*)N
II=RI+JJ*IMI

write(*,*)'partie réelle de l''integrale est',RI
write(*,*)'partie imaginaire de l''integrale est',IMI
write(*,*)' l''integrale complexe est: II=',II
end

double precision function Rf(x)
double precision x
Rf=(x*x-3)*dcos(x)/dsqrt(x-1)
return
end
double precision function IMf(x)
double precision x
IMf=(x*x-3)*dsin(x)/dsqrt(x-1)
return
end

Bonjour,

Un point de départ par ici : http://www.physics.unlv.edu/~pang/comp6/monte.f90

je voudrais s'il vous plait que quelqu'un me fasse un programme de simulation monte carlo avec le longage de programmation fortran qui calcul la probabilité de transmission a travrs un canal droit et plein en aluminium et entouré du fer en utilisant l'albédo.
les particules incidentes sont des neutrons monocinétique.
c'est urgent voila mon e-mail: fleila1980.hotmail.fr
Merci beaucoup.

cher itralover,
mon problème ne conserne pas le calcul d'integral mais de la probabilité de transmission à travers un canal droit et plein.
si tu veux me faire ce calcul ou tous personne pouvant m'aider ou capable de le faire qu'il me contacte sur mon e-mail pour lui donner les énnoncer du problème.
je serai très reconnaissante à toute personne pouvant m'aider à résoudre ce problème.
Meri chèrs amis