5 * Revision 1.1.1.2 1996/10/30 08:31:01 cernlib
8 * Revision 1.1.1.1 1996/04/12 15:29:49 plothow
12 #include "pdf/pilot.h"
14 SUBROUTINE STRDO1(X,SCALE,UPV,DNV,SEA,STR,CHM,GL)
15 C :::::::::::: DUKE OWENS SET 1 ::::::::::::::::::::::::::::
16 #include "pdf/impdp.inc"
17 #include "pdf/expdp.inc"
18 + F(5),A(6,5),B1(3,6,5)
19 DATA Q0,QL1/2.D0,.2D0/
20 DATA B1/3.D0,0.D0,0.D0,.419D0,.004383D0,-.007412D0,
21 &3.46D0,.72432D0,-.065998D0,4.4D0,-4.8644D0,1.3274D0,
23 &0.D0,0.D0,.763D0,-.23696D0,.025836D0,4.D0,.62664D0,-.019163D0,
24 &0.D0,-.42068D0,.032809D0,6*0.D0,1.265D0,-1.1323D0,.29268D0,
25 &0.D0,-.37162D0,-.028977D0,8.05D0,1.5877D0,-.15291D0,
26 &0.D0,6.3059D0,-.27342D0,0.D0,-10.543D0,-3.1674D0,
27 &0.D0,14.698D0,9.798D0,0.D0,.13479D0,-.074693D0,
28 &-.0355D0,-.22237D0,-.057685D0,6.3494D0,3.2649D0,-.90945D0,
29 &0.D0,-3.0331D0,1.5042D0,0.D0,17.431D0,-11.255D0,
30 &0.D0,-17.861D0,15.571D0,1.564D0,-1.7112D0,.63751D0,
31 &0.D0,-.94892D0,.32505D0,6.D0,1.4345D0,-1.0485D0,
32 &9.D0,-7.1858D0,.25494D0,0.D0,-16.457D0,10.947D0,
33 &0.D0,15.261D0,-10.085D0/
35 S= LOG( LOG( MAX(Q0,SCALE)/QL1)/ LOG(Q0/QL1))
38 10 A(J,I)=B1(1,J,I)+S*(B1(2,J,I)+S*B1(3,J,I))
40 40 F(I)=A(1,I)*X**A(2,I)*(WN-X)**A(3,I)*(WN+X*
41 & (A(4,I)+X*(A(5,I)+X*A(6,I))))
44 50 F(I)=F(I)*GAMFUN(AA)/((WN+A(2,I)*A(4,I)/AA)*GAMFUN(A(2,I))