COSMOS v7.655  COSMOSv7655
(AirShowerMC)
cmBremI2.f
Go to the documentation of this file.
1 ! implicit none
2 ! real*8 u, s
3 ! real*8 cmBremI2
4 ! u = 1.e-5
5 ! do while (u .lt. 300.)
6 ! s = cmBremI2(u)
7 ! write(*,*) u, s
8 ! u = u * 10.**.1
9 ! enddo
10 ! end
11  real*8 function cmbremi2(u)
12 ! This program give approximate integral value of
13 ! v^2 2zeta K(2zeta) dv from v = 0 to 1.
14 ! Original integaral is done by inteSpec1.f and
15 ! result was approximated by polynomials.
16 !
17  implicit none
18  real*8 u ! input. Upsilon
19 ! function value ! output. integral value
20 
21  real*8 lu, ans
22  if(u .lt. 3.16d-5) then
23 ! ans = ans * u**2 * 3.9257
24  ans = u**2 * 3.9257 ! 2012.Feb.4
25  elseif(u .lt. 1.584d-2) then
26  lu = log(u)
27  ans = ((-0.5219002e-02 *lu - 0.1067265)* lu
28  * -0.7421499 )*lu -0.7647935
29  ans = ans * u**2 * 3.9257
30  elseif(u .lt. 1.) then
31  lu = log(u)
32  ans = ((0.1696074e-01*lu + 0.1197285)*lu
33  * + 0.1439367e-01)*lu + 0.5595450e-01
34  ans = ans * u**2 * 3.9257
35  elseif(u .lt. 300.) then
36  lu = log(u)
37 
38  ans = (( 0.2566103e-02*lu -0.3215703e-01 )*lu
39  * + 0.9286942e-01)*lu + 0.1500733
40 
41  else
42  stop 55555
43  endif
44  cmbremi2 = ans
45  end
dE dx *! Nuc Int sampling table e
Definition: cblkMuInt.h:130
dE dx *! Nuc Int sampling table d
Definition: cblkMuInt.h:130
real *8 function cmbremi2(u)
Definition: cmBremI2.f:12