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