In this contribution, two general expressions for Markov parameters (MPs) and time moments (TMs) of continuous interval systems (CISs) are proposed. The TMs and MPs are derived in such a way that they do not require the inversion of the denominator of the transfer function of a given system. Additionally, there is no requirement for the solution of interval equations for obtaining TMs and MPs. The applicability of derived TMs and MPs is shown with the help of CIS model order reduction (MOR). The method utilized for the CIS MOR is the Routh–Padé approximation. The whole MOR procedure is explained with the help of one test case where the third-order system is reduced to first-order and second-order models. To establish the efficacy of the proposed method, the step and impulse responses of the system and model are plotted. For a fair comparative study, time-domain specifications such as the peak value, rise time, settling time, peak time and steady-state value for the lower and upper limits of the model are presented. Furthermore, the integral square error, integral absolute error, integral of time multiplied absolute error and integral of time multiplied square error for the difference between the responses of the system and model are tabulated. The results presented prove the applicability and effectiveness of the proposed method.