An algorithm to compute a low order single input-single output linear dynamic model from a higher-order model is presented. The algorithm is based on minimizing the integral of the squared error between the Step responses of the original system and the low-order system. The method utilizes an analytical form for the original model but the technique could ho extended to treat experimental data and to handle time delays. A numerical example shows that the algorithm is quite efficient