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`.
989
990
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
-
1006
991
```{code-cell} ipython3
1007
992
---
1008
993
mystnb:
1009
994
figure:
1010
995
caption: Stable Manifold and Phase Plane
1011
996
name: stable_manifold
1012
997
---
998
+
1013
999
fig, ax = plt.subplots(figsize=(7, 5))
1014
1000
1015
1001
K_range = np.arange(1e-1, 15, 0.1)
@@ -1058,8 +1044,21 @@ ax.set_ylabel('$C$')
1058
1044
plt.show()
1059
1045
```
1060
1046
1061
-
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.
1062
1058
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