A theory for linear multistep schemes applied to the initial value problem for a nonlinear first order system of differential equations with a singularity of the first kind is developed. Predictor-corrector schemes are also considered. The specific examples given are systems derived from partial differential equations in the presence of symmetry.