source:
PSPA/parmelaPSPA/trunk/xi0.f
@
476
Last change on this file since 476 was 12, checked in by , 12 years ago | |
---|---|
File size: 329 bytes |
Rev | Line | |
---|---|---|
[12] | 1 | function xi0(xsq) |
2 | c I0(kr) | |
3 | real y,p,t,s | |
4 | y=.25*xsq | |
5 | p=1. | |
6 | t=y | |
7 | s=1. | |
8 | 10 continue | |
9 | s=s+t | |
10 | if(t.lt.0.001*s)go to 20 | |
11 | p=p+1. | |
12 | t=t*y/p**2 | |
13 | go to 10 | |
14 | 20 continue | |
15 | xi0=s | |
16 | return | |
17 | end | |
18 | c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
Note: See TracBrowser
for help on using the repository browser.