1. A method for execution by a processor of converting a continuous-time model (CTM) to a discrete time model (DTM), the method comprising: determining a sampling time and CTM state-space matrices [A, B, C, D] defined for a dynamical system ordinary differential equation mathematical formulation; determining an order of approximation [i,j,k] parameters of the coefficients of an Obreshkov integration formula, to be used in determining the DTM [A dm , B dm , C dm , D dm ] to provide an approximation of accuracy and to meet A- or L-stability requirements defined for the dynamical system; and discretizing the CTM into a DTM using the selected order of approximation [i,j,k], and the CTM state-space matrices [A, B, C, D] and sampling time (T) by calculating Obreshkob-coefficient and matrices [A dm , B dm , C dm , D dm ] using Obreshkov-based polynomial matrices.