**************************************************************
*
*     example.f
*
*     Example program of the use of the subroutines 
*     in nnqns.f 
*     (See nnqns.f for detailed description of subroutines)
*
***************************************************************


      program example
      implicit none

      integer IPT,IALPHA,K

      double precision X,QNSAV,QNSER,XMIN,XMAX

      write(6,*) "-------------------------------------- "
      write(6,*) "    The NonSinglet PDF                 "
      write(6,*) "       ---------------                 "
      write(6,*) "   The NNPDF Collaboration             "
      write(6,*) " L. Del Debbio, S. Forte, J.I. Latorre,"
      write(6,*) "   A. Piccione and J. Rojo          "
      write(6,*) "                                       "
      write(6,*) "     See hep-ph/0701xxx for details    "
      write(6,*) "---------------------------------------"

      write(6,*) " "
      write(6,*) " Enter perturbative order "
      write(6,*) " 0 -> LO "
      write(6,*) " 1 -> NLO "
      write(6,*) " 2 -> NNLO "
      read(5,*) IPT

      write(6,*) " "
      if(IPT.eq.1) then
         write(6,*) "Select value of alphas(Mz)"
         write(6,*) " 0 -> alphas=0.116 "
         write(6,*) " 1 -> alphas=0.118 "
         write(6,*) " 2 -> alphas=0.120 "
         read(5,*) IALPHA
      else
         if(IPT.eq.0) write(6,*) "alphas = 0.130"
         if(IPT.eq.2) write(6,*) "alphas = 0.115"
      endif

      call INITPDF(IPT,IALPHA)

      write(6,*) " Enter range of x (xmin,xmax) "
      read(5,*) XMIN,XMAX
      write(6,*) "        "
      write(6,51) " x          qns(x)        Dqns(x)"
      do K=0,10
         X=XMIN*(XMAX/XMIN)**(dfloat(K)/10)
         call NNQNS(X,QNSAV,QNSER)
         write(6,50) X,QNSAV,QNSER
      enddo
      write(6,*) " "

 51   format(7x,a)
 50   format(1x,3(E12.4,2x))
      stop
      end

************************************************************