PDF/spdf/grsf2.F

   1 #include "pdf/pilot.h"
   2 C-----------------------------------------------------------------------
   3 #if defined(CERNLIB_SINGLE)
   4       function grsf2(x,s,alp,bet,a,b,ga,gb,gc,gd,ge,gep)
   5 #endif
   6 #if defined(CERNLIB_DOUBLE)
   7       double precision function grsf2(x,s,alp,bet,a,b,ga,gb,gc,gd,
   8      +                                ge,gep)
   9 #endif
  10 C
  11 #include "pdf/impdp.inc"
  12 C
  13       grsf2=(s*x**a*(ga+gb*sqrt(x)+gc*x**b)+
  14      +      s**alp*exp(-ge+sqrt(gep*s**bet*log(1.d0/x))))*
  15      +      (1.d0-x)**gd
  16       return
  17       end
  18