In this research proposal, diminution of higher order (HO) discrete interval system (DIS) is accomplished by utilizing Kharitonov polynomials. The DIS is firstly, transformed into continuous interval system (CIS). The Markov-parameters (MPs) and time-moments (TMs) are exploited for determination of approximated models. The ascertainment of model order diminution (MOD) of DISs is done by Routh-Padé approximation. The Routh table is utilized to obtain the denominator of approximated model. The unknown numerator coefficients of desired approximated model are determined by matching MPs and TMs of DISs and desired model. This whole procedure of MOD is elucidated with the help of one test illustration in which third order system is reduced to first order model as well as second order model. To prove applicability of the proposed method, impulse, step and Bode responses are plotted for both system and model. For relative comparison, time-domain specifications of proposed model are tabulated for both upper and lower limits. Further, performance indices are specified for dissimilarities between responses of system and model. The obtained results depict the effectiveness and efficacy for the proposed method.