We propose a new variational wavefunction to describe spin-dependent systems. Spin coherent states are used to represent the spin state of particles. Since states are parameterized by a continuous variable, correlation operators can be represented as simple integrals over c-valued functions. These integrals are simulated using standard Monte Carlo techniques. We present a general N-particle wavefunction written in the coherent state basis. Explicit Monte Carlo calculations on liquid He-3 using this wavefunction give very good agreement with standard methods.