double precision function andar(yy,xx,vvunif,vvnorm,n)
      implicit real*8(a-h,o-z)
      real*8 yy(n,1),xx(n),vvunif(n),vvnorm(n)
      s1 = zmeanc(yy,n,1,1)
      s2 = stdc(yy,n,1,1)
      do 100 i=1,n
 100     xx(i) = yy(i,1)
      call dqksort(xx,n)
      do 200 i=1,n
 200     vvunif(i) = qpnorm(xx(i),s1,s2)
      do 300 i=1,n
 300     xx(i) = -yy(i,1)
      call dqksort(xx,n)
      do 400 i=1,n
       xx(i) = -xx(i)
 400     vvnorm(i) = qpnorm(xx(i),s1,s2)
      a1 = 0.
      a2 = 0.
      do 500 i=1,n
      ii = (2 * i) - 1
      a1 = a1 + (ii * dlog(vvunif(i)))
      a2 = a2 + (ii * dlog(1 - vvnorm(i)))
 500  continue
      aa   = -((a1+a2) / n) - n
      andar = aa * (1 + (.75 / n) + (2.25 / (n * n)))
      return
      end



      double precision function qpnorm(w,s1,s2)
      implicit real*8(a-h,o-z)
      bb = cdfn((w - s1) / s2)
      qpnorm = bb
      return
      end