source: PSPA/parmelaPSPA/trunk/xi0.f @ 315

Last change on this file since 315 was 12, checked in by lemeur, 12 years ago

parmela pspa initial

File size: 329 bytes
Line 
1      function xi0(xsq)
2c    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
18c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
Note: See TracBrowser for help on using the repository browser.