A simple shuffle

      SUBROUTINE sub1(npts,seed,output)

      INTEGER npts,seed,output(npts)
      INTEGER magic1,magic2,j

      magic1 = 727
      magic2 = 524287
      output(1) = seed
      DO j = 2,npts
         output(j) = MOD(magic1*output(j - 1),magic2)
      END DO
      seed = MOD(magic1*output(npts),magic2)

      RETURN
      END