A numerical procedure for adding the effects of binary collisions to previously collisionless plasma simulation programs has been developed. It is shown that the model faithfully reproduces a numerical solution to the Fokker-Planck equation. As an illustration of the effect of collisions, a mirror mode problem is solved in θ-pinch geometry. Solutions are displayed for varying collision frequencies and qualitative as well as quantitative differences are readily apparent.