COSMOS v7.655  COSMOSv7655
(AirShowerMC)
kzeta2a.f File Reference

Go to the source code of this file.

Functions/Subroutines

real *8 function kzeta2a (a)
 

Function/Subroutine Documentation

◆ kzeta2a()

real*8 function kzeta2a ( real*8  a)

Definition at line 30 of file kzeta2a.f.

References e.

30  implicit none
31  real*8 a
32 !
33  real*8 c(7)/
34  * 29.55194, -260.5342, 955.5774, -1843.917, 1984.555, -1129.490,
35  * 265.9021/
36 
37  real*8 s, r, ak, tmp
38  integer k, j
39 !
40  s=0.
41  if(a .lt. .5) then
42  r=1.
43  if(a .eq. 0.) r=0.
44  ak=a
45  k=1
46  do while (r .gt. 1.e-3 )
47  tmp= ak/k**2
48  s=s+tmp
49  r=abs(tmp/s)
50  ak=ak*a
51  k=k+1
52  enddo
53  else
54  do j=7, 1, -1
55  s= s *a + c(j)
56  enddo
57  endif
58  kzeta2a=s
dE dx *! Nuc Int sampling table e
Definition: cblkMuInt.h:130
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz r
Definition: ZavoidUnionMap.h:1
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
Definition: cblkHeavy.h:36
real(4), save a
Definition: cNRLAtmos.f:20
real *8 function kzeta2a(a)
Definition: kzeta2a.f:30
dE dx *! Nuc Int sampling table c
Definition: cblkMuInt.h:130