A method is described and illustrated for explicit incorporation of and computation with ranges of initial conditions, functions, and parameter values in dynamic models using interval analysis. This approach is neither a statistical nor fuzzy set analysis but instead utilizes interval arithmetic which is particularly well suited for computerization. When a dynamic model is couched in interval analytic terms, ranges of all possible solutions are generated allowing not only an analysis of ranges of behavior modes but for sensitivity and stability analysis to be performed as a natural part of the model. Moreover, uncertainties such as specification, numerical method (e.g., numerical integration), and roundoff errors can also be analyzed in conjunction with or separate from the interval dynamic model.