  We present a Monte Carlo simulation for the response of the Laser
Interferometer Space Antenna (LISA) to the galactic gravitational wave
background. The simulated data streams are used to estimate the number and type
of binary systems that will be individually resolved in a 1-year power
spectrum. We find that the background is highly non-Gaussian due to the
presence of individual bright sources, but once these sources are identified
and removed, the remaining signal is Gaussian. We also present a new estimate
of the confusion noise caused by unresolved sources that improves on earlier
estimates.
