2 C-----------------------------------------------------------------------
3 #if defined(CERNLIB_SINGLE)
4 function grsf1(x,s,alp,bet,a,b,ga,gb,gc,gd,ge,gep)
6 #if defined(CERNLIB_DOUBLE)
7 double precision function grsf1(x,s,alp,bet,a,b,ga,gb,gc,gd,
11 #include "pdf/impdp.inc"
13 grsf1=(x**a*(ga+gb*sqrt(x)+gc*x**b)+
14 + s**alp*exp(-ge+sqrt(gep*s**bet*log(1.d0/x))))*