Trisurf Monte Carlo simulator
mihaf
2014-05-13 0f9d229fb830435d92cfc73c85aa7771a420286c
commit | author | age
5bb11d 1         program cod
SP 2         double precision  co(10,21),pi
3         integer lmax,l,al, am, m
4
5         pi=4.*datan(1.d0)
6
7         lmax=10
8         do 10 l=1,lmax
9             al=l
10             co(l,l+1)=dsqrt((2.*al+1.)/2./pi)
11             do 20 m=1,l-1
12                 am=m
13                 co(l,l+1+m)=co(l,l+m)*dsqrt(1.d0/(al-am+1)/(al+am))
14                 co(l,l+1-m)=co(l,l+1+m)
15 20          continue
16             co(l,2*l+1)=co(l,2*l)*dsqrt(1.d0/(2.*al))
17             co(l,1)=co(l,2*l+1)
18             co(l,l+1)=dsqrt((2.*al+1.)/4./pi)
19 10      continue
20
21         do 30 l=1,lmax
22             do 40 m=1,2*l+1
23                 print *, 'co(',l,',',m,')=',co(l,m)
24 40          continue
25 30      continue
26
27         end program cod
28