source:
PSPA/parmelaPSPA/trunk/xi0.f
@
405
Last change on this file since 405 was 12, checked in by , 12 years ago | |
---|---|
File size: 329 bytes |
Line | |
---|---|
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.