A geometric programming based computer-aided method to derive a reduced order (rth-order) approximant for a given (stable) SISO linear continuous-time system is presented. In this method, stability and the first r time moments/Markov parameters are preserved as well as the errors between a set of subsequent time moments/Markov parameters of the system and those of the model are minimized.