subroutine dsinqf (n,x,wsave)
      double precision x(1), wsave(1), xhold
c
      if (n .eq. 1) return
c
      ns2 = n/2
      do 101 k=1,ns2
         kc = n-k
         xhold = x(k)
         x(k) = x(kc+1)
         x(kc+1) = xhold
  101 continue
c
      call dcosqf (n,x,wsave)
c
      do 102 k=2,n,2
         x(k) = -x(k)
  102 continue
c
      return
      end