MODEL Procedure

Limitations

There are limitations to the types of differential equations that can be solved or estimated. One type is an explosive differential equation (finite escape velocity) for which the following differential equation is an example:

y Superscript Super Superscript prime Superscript Baseline equals a times y comma a greater-than 0

If this differential equation is integrated too far in time, y exceeds the maximum value allowed on the computer, and the integration terminates.

Likewise, differential systems that are singular cannot be solved or estimated in general. For example, consider the following differential system:

StartLayout 1st Row 1st Column x Superscript prime 2nd Column equals 3rd Column minus y Superscript Super Superscript prime Superscript Baseline plus 2 x plus 4 y plus exp left-parenthesis t right-parenthesis 2nd Row 1st Column y Superscript prime 2nd Column equals 3rd Column minus x Superscript Super Superscript prime Superscript Baseline plus y plus exp left-parenthesis 4 asterisk t right-parenthesis EndLayout

This system has an analytical solution, but an accurate numerical solution is very difficult to obtain. The reason is that y Superscript prime and x Superscript prime cannot be isolated on the left-hand side of the equation. If the equation is modified slightly to

StartLayout 1st Row 1st Column x Superscript prime 2nd Column equals 3rd Column minus y Superscript Super Superscript prime Superscript Baseline plus 2 x plus 4 y plus exp left-parenthesis t right-parenthesis 2nd Row 1st Column y Superscript prime 2nd Column equals 3rd Column x Superscript Super Superscript prime Superscript Baseline plus y plus exp left-parenthesis 4 t right-parenthesis EndLayout

then the system is nonsingular, but the integration process could still fail or be extremely slow. If the MODEL procedure encounters either system, a warning message is issued.

This system can be rewritten as the following recursive system, which can be estimated and simulated successfully with the MODEL procedure:

StartLayout 1st Row 1st Column x Superscript prime 2nd Column equals 3rd Column 0.5 y plus 0.5 exp left-parenthesis 4 t right-parenthesis plus x plus 1.5 y minus 0.5 exp left-parenthesis t right-parenthesis 2nd Row 1st Column y Superscript prime 2nd Column equals 3rd Column x Superscript Super Superscript prime Superscript Baseline plus y plus exp left-parenthesis 4 t right-parenthesis EndLayout

Petzold (1982) mentions a class of differential algebraic equations that, when integrated numerically, could produce incorrect or misleading results. An example of such a system is

StartLayout 1st Row 1st Column y 2 Superscript prime Baseline left-parenthesis t right-parenthesis 2nd Column equals 3rd Column y 1 left-parenthesis t right-parenthesis plus g 1 left-parenthesis t right-parenthesis 2nd Row 1st Column 0 2nd Column equals 3rd Column y 2 left-parenthesis t right-parenthesis plus g 2 left-parenthesis t right-parenthesis EndLayout

The analytical solution to this system depends on g and its derivatives at the current time only and not on its initial value or past history. You should avoid systems of this and other similar forms mentioned in Petzold (1982).

Last updated: June 19, 2025