cgUpsilon.f File Reference

Go to the source code of this file.

Functions/Subroutines
real *8 function	cgupsilon (u)

Function/Subroutine Documentation

cgupsilon()

real*8 function cgupsilon

(

real*8

u

)

Definition at line 15 of file cgUpsilon.f.

References e.

    implicit none
    real*8 u  ! input. upsilon

    real*8 z

    if(u .lt. 0.005) then
        cgupsilon = u ** 2 * (1. -5.953*u)
    elseif(u .lt. 100.) then
!                 polynomial approx. within 5 % error. roughly
        z = log(u)
                cgupsilon = 
     *          exp(((((-0.6305280e-04 * z -0.1159820e-03)*z 
     *           +0.9467500e-02 )*z-0.6861282e-01)*z  +1.080512)*z
     *           -1.770289)       
    else
          cgupsilon = 0.5563*u**0.66666666
    endif