Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I am working with a set of nonlinear dynamic equations that Mathematica has problems with solving. It is of the form $$ f_1(x_{t+1},y_{t+1},x_{t},y_t) = g_1(x_{t+1},y_{t+1},x_{t},y_t),\\ f_2(x_{t+1},y_{t+1},x_{t},y_t) = g_2(x_{t+1},y_{t+1},x_{t},y_t), $$ where all the function are nonlinear.

I would like to analyze steady states. I may easily solve that for $$x_{t+1}=x_t, \\ y_{t+1}=y_t,$$ but I have a problem with determining its local stability. Do you have any idea how to do that?

The standard procedure is to put all the future terms on the LHS and then evaluating the Jacobian of the RHS, but the thing is that the system is complicated enough even for Mathematica to find the analitical expression for the future variables.

share|cite|improve this question
up vote 1 down vote accepted

In this reply I'm assuming that $f_1,f_2,g_1,g_2$ all map from $\mathbb{R}^4$ to $\mathbb{R}$. Let $w_{t}:=[x_t,y_t]^T$ and

$h(w_{t+1},w_t):= \begin{bmatrix}f_1(x_{t+1},y_{t+1},x_t,y_t)-g_1(x_{t+1},y_{t+1},x_t,y_t)\\f_2(x_{t+1},y_{t+1},x_t,y_t)-g_2(x_{t+1},y_{t+1},x_t,y_t)\end{bmatrix}$.

So, you have the equation $h(w_{t+1},w_t)=0$ which defines pairs of $(w_{t+1},w_t)$. First issue that you need to address is whether this relation even defines a system, that is a function that maps $w_t$ to $w_{t+1}$.

For example, if $h$ has a unique root $w_{t+1}=\phi(w_t)$ (that is a unique function $\phi:\mathbb{R}^2\rightarrow \mathbb{R}^2$ such that $h(\phi(x),x)=0$ for all $x\in\mathbb{R}^2$), then you are looking at the system $w_{t+1}=\phi(w_t)$ and want to evaluate the Jacobian $\frac{\partial \phi}{\partial w_t}$ at the steady state(s).

Suppose there is a unique root $\phi$ and that $h$ is continuously differentiable at the steady state of interest $\bar{w}_t$, then we have that

$h'(\phi(w_t),w_t)=\frac{\partial h}{\partial \phi}(\phi(\bar{w}_t),\bar{w}_t)\frac{\partial \phi}{\partial w_t}(\phi(\bar{w}_t),\bar{w}_t)+\frac{\partial h}{\partial w_t}(\phi(\bar{w}_t),\bar{w}_t)=0$.

Suppose that $\frac{\partial h}{\partial \phi}(\phi(\bar{w}_t),\bar{w}_t)$ is not singular, then

$\frac{\partial \phi}{\partial w_t}(\phi(\bar{w}_t),\bar{w}_t)=-\left(\frac{\partial h}{\partial \phi}(\phi(\bar{w}_t),\bar{w}_t)\right)^{-1}\frac{\partial h}{\partial w_t}(\phi(\bar{w}_t),\bar{w}_t).$

Then to asses the local stability you just need to check the eigenvalues of the above Jacobian (note that you don't need anywhere an explicit formula for $\phi$ you just need to verify it's existence and uniqueness).

If there is more than one root, then one has to decide which root they're interested in (it'll be the one whose domain contains the initial conditions) and make sure that the domain on which the root is defined is forward invariant with respect to the system it defines.

Hope this helps.

share|cite|improve this answer
Great, thank you! I am trying it now – wlq Mar 12 '13 at 13:24
Well, there is still the same problemthat I may not represent the future state as a function of current state, in your example w[t+1] as a function of w[t]. The rest is clear and straightforward but since I don't have this representation, the problem remains... – wlq Mar 12 '13 at 16:02
You don't need to - you just need to compute the partial derivatives of $h$ with respect to $\phi$ (or $w_{t+1}$ as you had in the original equation), $\frac{\partial h}{\partial \phi}$ and the partial derivatives of $h$ with respect to $w_t$, $\frac{\partial h}{\partial w_t}$ evaluated at the steady state and use these to generate the desired Jacobian, $\frac{\partial \phi}{\partial w_t}$. – jkn Mar 12 '13 at 16:23
For example, suppose that $h:=\begin{bmatrix}x_{t+1}^3+y_{t+1}^2x_t+3\\ x_{t+1}^3y_t^2+2y_{t+1}+y_t\end{bmatrix}$. Then $\frac{\partial h}{\partial \phi} = \begin{bmatrix} 3x_{t+1}^2&2y_{t+1}x_t\\ 3x_{t+1}^2y_t^2& 2\end{bmatrix}$ and $\frac{\partial h}{\partial w_t} =\begin{bmatrix} y_{t+1}^2&0\\0&2y_tx_{t+1}^3+1\end{bmatrix}$. – jkn Mar 12 '13 at 16:32
That is correct, but it is difficult to prove that the relationship is unique. In your derivation you assume it, right? – wlq Mar 12 '13 at 16:48

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.