We present a scheme for systematically reducing the number of differential equations required for biophysically realistic neuron models. The techniques are general, are designed to be applicable to a large set of such models and retain in the reduced system as high a degree of fidelity to the original system as possible. As examples, we provide reductions of the Hodgkin-Huxley system and the A-current model of Connor et al. (1977).