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 |
|