*************************************************************
      program example
*************************************************************
*
*     calculation of BBG LO/NLO/NNLO/N3LO uv,dv & errors
*
*     example program
*
*************************************************************
*
      implicit none
*
      integer nq2,nxbj
      parameter (nq2=1 , nxbj=10)
*
      integer iini,iset
      integer iq2,ixbj
*
      double precision q2ref(nq2)
      double precision xbj(nxbj)
      double precision x,q2
      double precision uv,dv
      double precision duv,ddv
*    
      common /intini/ iini
*
*     set data:
*
      data q2ref  / 4.0d0 /
*
      data xbj / 1.0D-1, 2.0D-1, 3.0D-1, 4.0D-1, 5.0D-1, 
     +           6.0D-1, 7.0D-1, 8.0D-1, 9.0D-1, 9.5D-1 /
*
*     start reading data
*
      iini = 0                  ! reading data first time
*
*     choose order of parton sets and fit mode
*
c      iset = 1                  ! LO parton set
      iset = 2                  ! NLO parton set
c      iset = 3                  ! NNLO parton set
c      iset = 4                  ! N3LO parton set
*
*     write file header
*
      write(6,2000)
*
*     g1 and pdf vs x at Q2ref
*
      do iq2 = 1,nq2
*
         q2 = q2ref(iq2)
*
         do ixbj = 1,nxbj
*
            x = xbj(ixbj)
*
            call updf(iset, x, q2, uv, duv, dv, ddv)
*
*     output to file
*
            write(6,2001) x,q2,uv,duv,dv,ddv
*
         enddo
*
      enddo
*
 2000 format('* x,q2,uv,duv,dv,ddv')
 2001 format(6(1x,f9.6))
*
      end
*
      include 'updf.f'
*
*=======================================================================