subroutine scosqb (n,x,wsave)
      real x(1), wsave(1), tsqrt2, x1
      data tsqrt2 /  2.828427124 7461900976 0337744841 94 e0 /
c
      if (n-2) 101,102,103
  101 x(1) = 4.e0*x(1)
      return
c
  102 x1 = 4.e0*(x(1)+x(2))
      x(2) = tsqrt2*(x(1)-x(2))
      x(1) = x1
      return
c
  103 call scsqb1 (n,x,wsave,wsave(n+1))
c
      return
      end