You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
The following code generates Figure {numref}`stable_manifold`, which is patterned on a graph that appears on page 411 of {cite}`intriligator2002mathematical`.
990
989
990
+
Figure {numref}`stable_manifold` is a classic "phase plane" with "state" variable $K$ on the ordinate axis and "co-state" variable $C$ on the coordinate axis.
991
+
992
+
Figure {numref}`stable_manifold` plots three curves:
993
+
994
+
* the blue line graphs $C = \tilde C (K)$ of fixed points described by equation {eq}`eq:tildeC`.
995
+
* the red line graphs $K = \tilde K(C)$ of fixed points described by equation {eq}`eq:tildeK`
996
+
* the green line graphs the stable traced out by paths that converge to the steady state starting from an arbitrary $K_0$ at time $0$.
997
+
* for a given $K_0$, the shooting algorithm sets $C_0$ to the coordinate on the green line in order to initiate a path that converges to the optimal steady state
998
+
* the arrows on the green line show the direction in which dynamics {eq}`eq:systemdynamics` push successive $(K_{t+1}, C_t)$ pairs.
999
+
1000
+
In addition to the three curves, Figure {numref}`stable_manifold` plots arrows that point where the dynamics {eq}`eq:systemdynamics` drive the system when, for a given $K_0$, $C_0$ is not on the stable manifold depicted in the green line.
1001
+
1002
+
* If $C_0$ is set below the green line for a given $K_0$, too much capital is accumulated
1003
+
1004
+
* If $C_0$ is set above the green line for a given $K_0$, too little capital is accumulated
1005
+
991
1006
```{code-cell} ipython3
992
1007
---
993
1008
mystnb:
994
1009
figure:
995
1010
caption: Stable Manifold and Phase Plane
996
1011
name: stable_manifold
997
1012
---
998
-
999
1013
fig, ax = plt.subplots(figsize=(7, 5))
1000
1014
1001
1015
K_range = np.arange(1e-1, 15, 0.1)
@@ -1044,21 +1058,8 @@ ax.set_ylabel('$C$')
1044
1058
plt.show()
1045
1059
```
1046
1060
1047
-
Figure {numref}`stable_manifold` is a classic "phase plane" with "state" variable $K$ on the ordinate axis and "co-state" variable $C$ on the coordinate axis.
1048
-
1049
-
Figure {numref}`stable_manifold` plots three curves:
1050
-
1051
-
* the blue line graphs $C = \tilde C (K)$ of fixed points described by equation {eq}`eq:tildeC`.
1052
-
* the red line graphs $K = \tilde K(C)$ of fixed points described by equation {eq}`eq:tildeK`
1053
-
* the green line graphs the stable traced out by paths that converge to the steady state starting from an arbitrary $K_0$ at time $0$.
1054
-
* for a given $K_0$, the shooting algorithm sets $C_0$ to the coordinate on the green line in order to initiate a path that converges to the optimal steady state
1055
-
* the arrows on the green line show the direction in which dynamics {eq}`eq:systemdynamics` push successive $(K_{t+1}, C_t)$ pairs.
1056
-
1057
-
In addition to the three curves, Figure {numref}`stable_manifold` plots arrows that point where the dynamics {eq}`eq:systemdynamics` drive the system when, for a given $K_0$, $C_0$ is not on the stable manifold depicted in the green line.
1061
+
1058
1062
1059
-
* If $C_0$ is set below the green line for a given $K_0$, too much capital is accumulated
1060
-
1061
-
* If $C_0$ is set above the green line for a given $K_0$, too little capital is accumulated
0 commit comments