A novel order reduction approach for continuous time systems which utilizes recently developed Cuckoo Search Optimization and Routh Approximation (RA) technique is proposed to simplify transfer functions/transfer matrices. The method may be applied to continuous time systems and assures stability of lower order model when the given higher dimension system is stable. The efficient approach is provided for reduction of MIMO and SISO Linear Time Invariant (LTI) systems. To describe the proposed method two numerical examples are solved. The results obtained are judged with other popular order reduction techniques in respect of Integral Square Error (ISE). Finally, the steady state value of reduced system is matched with the given original high dimensional system.