A matrix method is used to solve the difference equations which describe the zone refining of a bar. The matrix method is designed for direct numerical calculation of the solute distribution, after any number of passes, with any initial solute distribution. The method is illustrated with a simple example.