*
* $Id$
*
* $Log$
* Revision 1.1.1.2  1996/10/30 08:31:02  cernlib
* Version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:29:49  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
C-----------------------------------------------------------------------
      SUBROUTINE STRDO2(X,SCALE,UPV,DNV,SEA,STR,CHM,GL)
C :::::::::::: DUKE OWENS SET 2 ::::::::::::::::::::::::::::
C
#include "pdf/impdp.inc"
#include "pdf/expdp.inc"
     +       F(5),A(6,5),B2(3,6,5)
      DATA Q0,QL2/2.D0,.4D0/
      DATA B2/3.D0,0.D0,0.D0,.3743D0,.013946D0,-.00031695D0,
     &3.329D0,.75343D0,-.076125D0,6.032D0,-6.2153D0,1.5561D0,
     &6*0.D0,1.D0,0.D0,
     &0.D0,.7608D0,-.2317D0,.023232D0,3.83D0,.62746D0,-.019155D0,
     &0.D0,-.41843D0,.035972D0,6*0.D0,1.6714D0,-1.9168D0,.58175D0,
     &0.D0,-.27307D0,-.16392D0,9.145D0,.53045D0,-.76271D0,
     &0.D0,15.665D0,-2.8341D0,0.D0,-100.63D0,44.658D0,
     &0.D0,223.24D0,-116.76D0,0.D0,.067368D0,-.030574D0,
     &-.11989D0,-.23293D0,-.023273D0,3.5087D0,3.6554D0,-.45313D0,
     &0.D0,-.47369D0,.35793D0,0.D0,9.5041D0,-5.4303D0,
     &0.D0,-16.563D0,15.524D0,.8789D0,-.97093D0,.43388D0,
     &0.D0,-1.1612D0,.4759D0,4.D0,1.2271D0,-.25369D0,
     &9.D0,-5.6354D0,-.81747D0,0.D0,-7.5438D0,5.5034D0,
     &0.D0,-.59649D0,.12611D0/
      WN=1.D0
      S= LOG( LOG( MAX(Q0,SCALE)/QL2)/ LOG(Q0/QL2))
      DO 10 I=1,5
         DO 10 J=1,6
   10 A(J,I)=B2(1,J,I)+S*(B2(2,J,I)+S*B2(3,J,I))
      DO 20 I=1,5
   20 F(I)=A(1,I)*X**A(2,I)*(WN-X)**A(3,I)*(WN+X* (A(4,I)+X*(A(5,I)+X*A
     +(6,I))))
      DO 30 I=1,2
         AA=WN+A(2,I)+A(3,I)
   30 F(I)=F(I)*GAMFUN(AA)/((WN+A(2,I)*A(4,I)/AA)*GAMFUN(A(2,I))
     +*GAMFUN(WN+A(3,I)))
      UPV=F(1)-F(2)
      DNV=F(2)
      SEA=F(3)/6.D0
      STR=SEA
      CHM=F(4)
      GL =F(5)
      RETURN
      END