与えられた関数
\[\text{(消費関数)} \quad C = 100 + 0.8(Y - T) \tag{1}\]
\[\text{(投資関数)} \quad I = 200 - 10r \tag{2}\]
\[\text{(政府支出)} \quad G = 200 \tag{3}\]
\[\text{(税)} \quad T = 100 \tag{4}\]
\[\text{(貨幣需要関数)} \quad L = 0.25Y - 20r \tag{5}\]
\[\text{(貨幣供給)} \quad M = 400 \tag{6}\]
\[\text{(価格水準)} \quad P = 1 \tag{7}\]
IS曲線の導出
(1),(2),(3),(4)より,財市場の均衡は以下のように書ける.
\[
Y = C + I + G = 100 + 0.8(Y - 100) + (200 - 10r) + 200\tag{8}
\]
展開すると,
\[
Y = 100 + 0.8Y - 80 + 200 - 10r + 200 = 0.8Y + 420 - 10r
\]
\[
0.2Y = 420 - 10r \Rightarrow Y = 2100 - 50r
\]
となる.したがって,IS曲線は以下のように表される:
\[
Y = 2100 - 50r\tag{9}
\]
LM曲線の導出
(5),(6),(7)より,貨幣市場の均衡は以下のように書ける.
\[
\frac{M}{P} = 0.25Y - 20r \Rightarrow 400= 0.25Y - 20r\tag{10}
\]
展開すると,
\[
Y = 1600 + 80r\tag{11}
\]
となる.
IS曲線とLM曲線の交点(均衡点)
均衡は IS = LM の交点として求められる.(8),(9)より,
\[
2100 - 50r = 1600 + 80r
\]
\[
500 = 130r \Rightarrow r = \frac{500}{130} \approx 3.85
\]
と計算できる.
この値をIS曲線((8))に代入して \( Y \) を求めると,
\[
Y = 2100 - 50 \cdot 3.85 \approx 1907.5
\]
となる.
均衡点の結論
したがって,IS-LMモデルにおける均衡は,
\[
r^* \approx 3.85,\quad Y^* \approx 1907.5
\]
と計算できる.
Mathematicaのプログラム
📋 コピー
(*パラメータ設定*)
T = 100;
G = 200;
M = 400;
P = 1;
(*関数定義*)
Cons[Y_] := 100 + 0.8 (Y - T);
Inv[r_] := 200 - 10 r;
L[Y_, r_] := 0.25 Y - 20 r;
(*IS曲線:Y=C+I+G を rについて解く\[RightArrow]r=f(Y)*)
isEq = Y == Cons[Y] + Inv[r] + G;
isSol = Solve[isEq, r][[1]]; (*r->... の形*)
(*LM曲線:L(Y,r)=M/P を rについて解く\[RightArrow]r=f(Y)*)
lmEq = L[Y, r] == M/P;
lmSol = Solve[lmEq, r][[1]]; (*r->... の形*)
(*交点を数値で求める*)
intersection =
Solve[{Y == Cons[Y] + Inv[r] + G, L[Y, r] == M/P}, {Y, r}][[1]];
yStar = Y /. intersection;
rStar = r /. intersection;
(*均衡点を表示*)
Print["均衡国民所得 Y* = ", N[yStar]];
Print["均衡利子率 r* = ", N[rStar]];
(*グラフ描画範囲を指定*)
yMin = 1000;
yMax = 3000;
graph = Plot[Evaluate[{r /. isSol, r /. lmSol}], {Y, yMin, yMax},
PlotRange -> {{yMin, yMax}, {0, 10}}, PlotStyle -> {Blue, Green},
PlotLegends -> {"IS曲線", "LM曲線"}, AxesLabel -> {"Y(国民所得)",
"r(利子率)"},
PlotLabel -> "IS-LMモデルにおける均衡(縦軸:利子率)", ImageSize -> Large,
Epilog -> {Red, PointSize[Large], Point[{yStar, rStar}], Black,
Text[Style["均衡点", 12], {yStar + 50, rStar + 0.3}]}]
(*PNG形式で保存*)
path = "/Users/macmini2025/Documents/\
Mathematica/";
filename = "ISLM_Equilibrium2.png";
(*ディレクトリが存在するかチェックして、存在すれば保存*)
If[
DirectoryQ[path],
Export[path <> filename, graph],
Print["保存に失敗:ディレクトリが存在しません \[RightArrow] ", path]]
IS-LMモデルを表示