subroutine fits
      print *, 'the following functions are provided:'
      print *, '1) poly(n,x): fit function for nth order poly.'
      print *, '   polyf(n,x): plot function, uses vector par()'
      print *, '2) pole2(x): 2-pole fit function'
      print *, '   pole2f(x): 2-pole plot function, uses par()'
      end

c polynomial of order n

      real function poly(n,x)
      integer n
      real x
      common /pawpar/p(99)
      integer i
      poly=p(n+1)
      do i=n,1,-1
        poly=poly*x+p(i)
      enddo
      end

      real function polyf(n,x)
      integer n
      real x
      common /pawpar/p(99)
      vector par(99)
      do i=1,n+1
        p(i)=par(i)
      enddo
      polyf=poly(n,x)
      end

c two-pole fit function

      real function pole2(x)
      real x
      common /pawpar/p(99)
      integer order
      parameter (order=1)
      integer i
      real xx
      xx=exp(x)
      pole2=p(order+1)
      do i=order,1,-1
        pole2=pole2*xx**2+p(i)
      enddo
      pole2=pole2*(p(3)**p(4)/(p(3)**p(4)+xx**p(4)))
      end

      real function pole2f(x)
      real x
      common /pawpar/p(99)
      vector par(99)
      integer i
      do i=1,9
        p(i)=par(i)
      enddo
      pole2f=pole2(log(x))
      end