  We simulate TIP3P water using a constrained Monte Carlo algorithm to generate
electrostatic interactions eliminating the need to sum over long ranged Coulomb
interactions. We study discretization errors when interpolating charges using
splines and Gaussians. We compare our implementation to molecular dynamics and
Brownian dynamics codes.
