A numerical technique is introduced for the solution of the time-dependent equation for the solar modulation of galactic cosmic rays, in which transport in heliographic latitude is ignored. We illustrate this method with a model for the solar cycle variation in cosmic ray intensity. The variation is assumed to result from changes in the number of enhanced comsic ray scattering regions, which are produced by flare-generated shock waves in the solar wind. The model can account for the observed solar cycle variation in intensity at different energies as well as the observed spatial gradients. This model also provides a natural explanation for the cosmic ray hysteresis effect.