In this paper, the modified matrix exponential (MME) method under zero-order hold (ZOH) assumption, is applied to solve systems of stiff ordinary differential equations. Some examples are given to illustrate the accuracy and effectiveness of the method. We compare our results with results obtained by matrix exponential (ME) method and by the Matlab ode23 solver.
