We simulate two competing firms in a Cournot duopoly market structure where each firm decides its production quantity ($q_1$ and $q_2$).

The firms interact under the following assumptions:

• Inverse Demand Function:
  $( P = a - b(Q) )$, where $(Q = q_1 + q_2)$ is the total market quantity.
• Cost Function:
  $(C(q_i) = c \cdot q_i)$ (linear cost).

Given:

• $(a = 100)$: Intercept of demand curve.
• $(b = 2)$: Slope of demand curve.
• $(c = 20)$: Marginal cost for each firm.

Objective:
Find the Nash equilibrium where both firms choose output quantities that maximize their profits.

Python Implementation

We solve this problem using an optimization approach to find the equilibrium quantities for both firms.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67

import numpy as np
from scipy.optimize import minimize
import matplotlib.pyplot as plt

# Parameters
a = 100  # Demand intercept
b = 2    # Demand slope
c = 20   # Marginal cost

# Profit function for Firm 1
def profit_firm1(q):
    q1, q2 = q  # Quantities produced by Firm 1 and Firm 2
    P = a - b * (q1 + q2)  # Price
    profit = (P * q1) - (c * q1)
    return -profit  # Negative for maximization

# Profit function for Firm 2 (symmetric)
def profit_firm2(q):
    q1, q2 = q
    P = a - b * (q1 + q2)
    profit = (P * q2) - (c * q2)
    return -profit

# Nash equilibrium condition: both firms maximize profit simultaneously
def nash_equilibrium(q):
    q1, q2 = q
    grad1 = (a - b * (q1 + q2)) - b * q1 - c  # First-order condition for Firm 1
    grad2 = (a - b * (q1 + q2)) - b * q2 - c  # First-order condition for Firm 2
    return [grad1, grad2]

# Initial guess for quantities
initial_guess = [10, 10]

# Solve for equilibrium quantities
result = minimize(lambda q: np.sum(np.square(nash_equilibrium(q))), initial_guess, bounds=[(0, None), (0, None)])
q1_opt, q2_opt = result.x

# Compute price and profits
Q_total = q1_opt + q2_opt
P_eq = a - b * Q_total
profit1 = (P_eq * q1_opt) - (c * q1_opt)
profit2 = (P_eq * q2_opt) - (c * q2_opt)

# Display results
print(f"Optimal Quantity (Firm 1): {q1_opt:.2f}")
print(f"Optimal Quantity (Firm 2): {q2_opt:.2f}")
print(f"Market Price: {P_eq:.2f}")
print(f"Profit (Firm 1): {profit1:.2f}")
print(f"Profit (Firm 2): {profit2:.2f}")

# Visualization
q1_range = np.linspace(0, 30, 100)
q2_best_responses = [(a - c - b * q) / (2 * b) for q in q1_range]
q1_best_responses = [(a - c - b * q) / (2 * b) for q in q1_range]

plt.figure(figsize=(8, 6))
plt.plot(q1_range, q2_best_responses, label="Firm 2's Best Response", color="blue")
plt.plot(q1_best_responses, q1_range, label="Firm 1's Best Response", color="orange")
plt.scatter(q1_opt, q2_opt, color="red", zorder=5, label=f"Nash Equilibrium ({q1_opt:.2f}, {q2_opt:.2f})")

# Labels and legend
plt.xlabel("Firm 1 Quantity (q1)")
plt.ylabel("Firm 2 Quantity (q2)")
plt.title("Nash Equilibrium in Cournot Duopoly")
plt.legend()
plt.grid()
plt.show()

Code Explanation

1. Parameters:

   • a, b: Define the linear demand curve $(P = a - bQ)$.
   • c: Marginal cost per unit.

2. Profit Maximization:

   • The profit for each firm depends on the price and their quantity.
     First-order conditions (gradients) are derived and solved numerically.

3. Optimization:

   • We use the scipy.optimize.minimize function to find Nash equilibrium quantities $(q_1, q_2)$.

4. Visualization:

   • Best-response curves for Firm 1 and Firm 2 are plotted.
   • The Nash equilibrium is where the two curves intersect.

Explanation of Results

Optimal Quantity (Firm 1): 13.33
Optimal Quantity (Firm 2): 13.33
Market Price: 46.67
Profit (Firm 1): 355.56
Profit (Firm 2): 355.56

1. Optimal Quantities:
   Both firms (Firm 1 and Firm 2) produce 13.33 units of output.
   This is the symmetric Nash Equilibrium in a Cournot Duopoly, where neither firm has an incentive to unilaterally change its production.

2. Market Price:
   Given the total output of $( Q = 26.67 )$ ($13.33$ from each firm), the market price is determined to be 46.67.

3. Profits:
   Both firms earn a profit of 355.56.
   The symmetry reflects the identical cost structures and competitive conditions.

Graph Interpretation

• Best Response Curves:

  • The blue line represents Firm 2’s best response curve, showing Firm 2’s optimal output for each quantity produced by Firm 1.
  • The orange line represents Firm 1’s best response curve.

• Nash Equilibrium:

  • The intersection of the two best response curves, marked with a red dot, represents the equilibrium where both firms’ strategies are mutually optimal.

• At the equilibrium point $(13.33, 13.33)$, both firms’ outputs are stable, and no firm can improve its profit by unilaterally changing its output.

Problem Statement

A firm uses two inputs, labor $(L)$ and capital $(K)$, to produce output $(Q)$ using a $Cobb$-$Douglas$ production function:
$$
Q(L, K) = A \cdot L^{\alpha} \cdot K^{\beta}
$$
where $(A)$ is total factor productivity, and $(\alpha)$ and $(\beta)$ are output elasticities of labor and capital, respectively.

Given:

• $(A = 1)$, $(\alpha = 0.4)$, $(\beta = 0.6)$
• Budget for inputs: $(B = 1000)$
• Prices: $(p_L = 10)$ (per unit of labor), $(p_K = 20)$ (per unit of capital)

Objective:
Maximize output $Q(L, K)$ subject to the budget constraint:
$$
10L + 20K \leq 1000
$$

Python Implementation

We solve this optimization problem using SciPy.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62

import numpy as np
from scipy.optimize import minimize
import matplotlib.pyplot as plt

# Define the production function (to be maximized)
def production(vars):
    L, K = vars
    A, alpha, beta = 1, 0.4, 0.6
    return -(A * (L**alpha) * (K**beta))  # Negative for minimization

# Budget constraint
def budget_constraint(vars):
    L, K = vars
    return 1000 - (10 * L + 20 * K)  # Budget must not be exceeded

# Bounds and initial guess
bounds = [(0, None), (0, None)]  # Non-negative inputs
initial_guess = [10, 10]

# Constraint dictionary
constraints = {'type': 'ineq', 'fun': budget_constraint}

# Solve the optimization problem
result = minimize(production, initial_guess, bounds=bounds, constraints=constraints)

# Extract optimal values
L_opt, K_opt = result.x

# Display results
print("Optimal labor (L):", L_opt)
print("Optimal capital (K):", K_opt)
print("Maximum output (Q):", -result.fun)

# Visualization
L = np.linspace(0, 100, 100)
K = (1000 - 10 * L) / 20  # Budget constraint line
def production_curve(L, production_value):
    return (production_value / (1 * (L**0.4)))**(1/0.6)

plt.figure(figsize=(8, 6))

# Budget constraint line
plt.plot(L, K, label="Budget Constraint", color="blue")

# Isoquants (levels of Q)
production_values = [-result.fun * 0.6, -result.fun * 0.8, -result.fun]  # Example production levels
for q in production_values:
    plt.plot(L, production_curve(L, q), label=f"Output = {q:.2f}", linestyle="--")

# Optimal point
plt.scatter(L_opt, K_opt, color="red", label="Optimal Point", zorder=5)
plt.text(L_opt + 1, K_opt, f"({L_opt:.2f}, {K_opt:.2f})")

# Labels and legend
plt.xlabel("Labor (L)")
plt.ylabel("Capital (K)")
plt.title("Production Maximization with Budget Constraint")
plt.axhline(0, color="black", linewidth=0.5, linestyle="--")
plt.axvline(0, color="black", linewidth=0.5, linestyle="--")
plt.legend()
plt.grid()
plt.show()

Explanation

1. Production Function:
   The $Cobb$-$Douglas$ function is modeled as negative for minimization since scipy.optimize.minimize minimizes objectives by default.

2. Budget Constraint:
   The constraint ensures the total cost of labor $(10L)$ and capital $(20K)$ does not exceed the budget.

3. Optimization:
   SciPy‘s minimize finds the optimal levels of $(L)$ and $(K)$ that maximize production $Q(L, K)$ within the budget.

4. Visualization:

   • The budget constraint is plotted as a line.
   • Isoquants represent levels of production.
   • The optimal point is shown at the intersection of the budget line and the highest achievable isoquant.

Explanation of Results

Optimal labor (L): 39.99936839032418
Optimal capital (K): 30.000315804863785
Maximum output (Q): 33.658654360606

1. Optimal Values:

   • Labor (L): 40.00
   • Capital (K): 30.00
   • Maximum Output (Q): 33.66

2. Budget Constraint:
   The blue line represents the budget constraint $(10L + 20K = 1000)$, showing the feasible combinations of labor and capital.

3. Isoquants:
   The dashed curves (orange, green, blue) are isoquants representing different output levels.
   The highest isoquant achievable ($33.66$) touches the budget constraint at the optimal point.

4. Optimal Point:
   The red dot at $(L = 40, K = 30)$ indicates the optimal allocation of inputs that maximizes output given the budget constraint.

Here’s an example of a Consumer Utility Maximization problem and its $Python$ implementation.

Problem Statement

A consumer has a utility function:
$$
U(x, y) = x^{0.5} \cdot y^{0.5}
$$
where $(x)$ and $(y)$ are the quantities of two goods consumed.
The consumer has a budget $(B = 100)$, and the prices of the goods are $(p_x = 2)$ and $(p_y = 5)$.

Objective:
Maximize $U(x, y)$ subject to the budget constraint:
$$
2x + 5y \leq 100
$$

Python Code Implementation

Here, we will solve the problem using SciPy‘s optimization tools and visualize the result.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61

import numpy as np
from scipy.optimize import minimize
import matplotlib.pyplot as plt

# Define the utility function (to be maximized)
def utility(vars):
    x, y = vars
    return -(x**0.5 * y**0.5)  # Negative for minimization

# Budget constraint
def budget_constraint(vars):
    x, y = vars
    return 100 - (2 * x + 5 * y)  # Budget must not be exceeded

# Bounds and initial guess
bounds = [(0, None), (0, None)]  # Non-negative consumption
initial_guess = [1, 1]

# Constraint dictionary
constraints = {'type': 'ineq', 'fun': budget_constraint}

# Solve the optimization problem
result = minimize(utility, initial_guess, bounds=bounds, constraints=constraints)

# Extract optimal values
x_opt, y_opt = result.x

# Display results
print("Optimal consumption of good x:", x_opt)
print("Optimal consumption of good y:", y_opt)
print("Maximum utility:", -result.fun)

# Visualization
x = np.linspace(0, 50, 100)
y = (100 - 2 * x) / 5  # Budget constraint line
def utility_curve(x, utility_value):
    return (utility_value / np.sqrt(x))**2

plt.figure(figsize=(8, 6))

# Budget constraint line
plt.plot(x, y, label="Budget Constraint", color="blue")

# Indifference curves
utility_values = [-result.fun * 0.6, -result.fun * 0.8, -result.fun]  # Example utility levels
for u in utility_values:
    plt.plot(x, utility_curve(x, u), label=f"Utility = {u:.2f}", linestyle="--")

# Optimal point
plt.scatter(x_opt, y_opt, color="red", label="Optimal Point", zorder=5)
plt.text(x_opt + 1, y_opt, f"({x_opt:.2f}, {y_opt:.2f})")

# Labels and legend
plt.xlabel("Good x")
plt.ylabel("Good y")
plt.title("Utility Maximization with Budget Constraint")
plt.axhline(0, color="black", linewidth=0.5, linestyle="--")
plt.axvline(0, color="black", linewidth=0.5, linestyle="--")
plt.legend()
plt.grid()
plt.show()

Explanation

1. Utility Function:
   The function $U(x, y) = x^{0.5} \cdot y^{0.5}$ is modeled as negative because optimization functions typically minimize objectives.

2. Budget Constraint:
   $(2x + 5y \leq 100)$ is converted into a function for scipy.optimize.minimize.

3. Optimization:
   The minimize function from SciPy solves for $(x)$ and $(y)$ that maximize utility within the budget.

4. Visualization:
   The budget constraint and indifference curves (utility levels) are plotted.
   The optimal point is marked for clarity.

Result Graph

Optimal consumption of good x: 24.99849398099728
Optimal consumption of good y: 10.000602407600788
Maximum utility: 15.8113882721523

The graph shows:

• The budget constraint as a line.
• Indifference curves for different utility levels.
• The optimal consumption bundle at the intersection of the budget constraint and the highest achievable indifference curve.

A Numerical Example with Python Implementation

The Principal-Agent problem involves designing a contract such that the agent performs actions aligning with the principal’s interests, despite having different objectives.

This example will focus on a simple contract design where a principal hires an agent to perform a task, and the agent’s effort influences the outcome.

Problem Setup

1. Principal:
   • Wants to maximize profit: $ \pi = R - w $, where $ R $ is revenue (affected by agent’s effort $ e $) and $ w $ is the wage paid.
2. Agent:
   • Maximizes utility: $ U = w - C(e) $, where $ C(e) $ is the cost of effort $ e $.
   • Cost function: $ C(e) = \frac{1}{2} e^2 $.
3. Revenue Function:
   • $ R = \theta e $, where $ \theta > 0 $ is the productivity factor.
4. Participation Constraint:
   • Agent’s utility $ U \geq \bar{U} $, where $ \bar{U} $ is the reservation utility.
5. Incentive Compatibility Constraint:
   • The agent chooses effort $ e $ to maximize $ U $.

The principal’s problem is to choose $ w(e) $ and $ e $ to maximize $ \pi $ while satisfying the constraints.

Python Implementation

We will solve this using numerical optimization with SciPy.

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80

import numpy as np
from scipy.optimize import minimize
import matplotlib.pyplot as plt

# Parameters
theta = 2.0  # Productivity factor
U_bar = 1.0  # Reservation utility

# Cost function and its derivative
def cost(e):
    return 0.5 * e**2

def cost_prime(e):
    return e

# Revenue function
def revenue(e):
    return theta * e

# Agent's utility function
def utility(w, e):
    return w - cost(e)

# Principal's profit function
def profit(e, w):
    return revenue(e) - w

# Objective function for the principal (negative profit for minimization)
def objective(vars):
    e, w = vars
    return -profit(e, w)

# Constraints
def participation_constraint(vars):
    e, w = vars
    return utility(w, e) - U_bar

def incentive_compatibility_constraint(vars):
    e, w = vars
    return theta - cost_prime(e)  # Agent maximizes utility when this equals zero

# Initial guess
e0 = 1.0  # Initial effort guess
w0 = 2.0  # Initial wage guess
initial_guess = [e0, w0]

# Bounds (effort and wage should be non-negative)
bounds = [(0, None), (0, None)]

# Constraints as a dictionary
constraints = [
    {"type": "ineq", "fun": participation_constraint},
    {"type": "eq", "fun": incentive_compatibility_constraint}
]

# Solve the optimization problem
result = minimize(objective, initial_guess, bounds=bounds, constraints=constraints)

# Extract results
e_opt, w_opt = result.x

# Display results
print(f"Optimal effort (e*): {e_opt:.2f}")
print(f"Optimal wage (w*): {w_opt:.2f}")
print(f"Principal's profit: {-result.fun:.2f}")

# Visualization
# Effort vs. Profit
effort_values = np.linspace(0, 2 * e_opt, 100)
profit_values = [profit(e, revenue(e) - cost(e)) for e in effort_values]

plt.figure(figsize=(10, 6))
plt.plot(effort_values, profit_values, label="Profit vs. Effort", color="blue")
plt.axvline(x=e_opt, color="red", linestyle="--", label=f"Optimal Effort (e* = {e_opt:.2f})")
plt.xlabel("Effort (e)")
plt.ylabel("Profit")
plt.title("Principal's Profit as a Function of Effort")
plt.legend()
plt.grid()
plt.show()

Explanation

1. Optimization Objective:

   • The principal seeks to maximize profit while ensuring the agent’s participation and incentivizing optimal effort.

2. Constraints:

   • Participation ensures the agent’s utility meets a minimum threshold.
   • Incentive Compatibility ensures the agent’s choice of effort maximizes their utility.

3. Results:

   • Optimal effort $(e^*)$ and wage $(w^*)$ are calculated.
   • The graph shows the relationship between effort and profit, with the optimal effort marked.

Optimal effort (e*): 2.00
Optimal wage (w*): 3.00
Principal's profit: 1.00

Phillips Curve

The Phillips Curve illustrates the inverse relationship between the rate of unemployment and the rate of inflation within an economy.

Let’s create a simple example where we model this relationship using synthetic data.

We’ll use $Python$ to perform the calculations and visualize the results.

Step 1: Create Synthetic Data

First, we’ll generate some synthetic data for unemployment rates and corresponding inflation rates.

Step 2: Analyze and Plot the Data

We’ll analyze this data to see the relationship and plot the Phillips Curve.

Here is the $Python$ code for generating the data, analyzing it, and plotting the curve:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import linregress

# Step 1: Create synthetic data
np.random.seed(42)  # for reproducibility
unemployment_rate = np.random.uniform(3, 10, 50)  # Unemployment rate between 3% and 10%
inflation_rate = np.random.uniform(0, 5, 50) - 0.5 * unemployment_rate  # Simple inverse relationship

# Add some noise to make the data more realistic
inflation_rate += np.random.normal(0, 0.5, inflation_rate.shape)

# Step 2: Analyze and plot the data
slope, intercept, r_value, p_value, std_err = linregress(unemployment_rate, inflation_rate)

# Calculate regression line
regression_line = intercept + slope * unemployment_rate

# Plotting
plt.scatter(unemployment_rate, inflation_rate, label='Data points')
plt.plot(unemployment_rate, regression_line, color='red', label=f'Regression line: y={slope:.2f}x+{intercept:.2f}')
plt.xlabel('Unemployment Rate (%)')
plt.ylabel('Inflation Rate (%)')
plt.title('Phillips Curve')
plt.legend()
plt.grid(True)
plt.show()

Explanation

1. Data Generation:
   We create synthetic data for the unemployment and inflation rates.
   The relationship is modeled with some random noise to simulate real-world data.
2. Regression Analysis:
   We use linear regression to find the best-fit line that describes the relationship between unemployment and inflation rates.
3. Plotting:
   We plot the data points and the regression line to visualize the Phillips Curve.

Graph Interpretation

• The scatter plot shows individual data points representing different unemployment and inflation rate pairs.
• The red line is the regression line, indicating the trend or the average relationship between unemployment and inflation rates.
  According to the Phillips Curve, as the unemployment rate decreases, the inflation rate tends to increase, and this inverse relationship is depicted by the slope of the line.

Example Problem

We aim to simulate a simplified BoP analysis to understand the interactions between the Current Account (CA), Capital Account (KA), and the Financial Account (FA).

Problem

1. The Current Account (CA) records trade in goods and services, income, and transfers.
2. The Capital Account (KA) reflects capital transfers and acquisition/disposal of non-produced assets.
3. The Financial Account (FA) captures investments in financial assets.

By the BoP identity:
$$
\text{CA} + \text{KA} + \text{FA} = 0
$$

We will simulate and visualize:

• Changes in the CA due to trade balance adjustments.
• How these changes are offset by the KA and FA to maintain equilibrium.

Python Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44

import numpy as np
import matplotlib.pyplot as plt

# Parameters
time_steps = 20  # Number of periods
initial_CA = -2.0  # Initial Current Account balance (deficit, in $ billion)
initial_KA = 0.5   # Initial Capital Account balance (in $ billion)
initial_FA = 1.5   # Initial Financial Account balance (in $ billion)
trade_balance_changes = np.linspace(-0.5, 1.0, time_steps)  # Simulated trade balance changes

# BoP components
CA = [initial_CA]
KA = [initial_KA]
FA = [initial_FA]

# Simulate BoP dynamics
for i in range(1, time_steps):
    # Update Current Account (CA) due to trade balance changes
    new_CA = CA[-1] + trade_balance_changes[i]
    CA.append(new_CA)
    
    # Capital and Financial Accounts adjust to maintain BoP equilibrium
    new_FA = -new_CA - KA[-1]  # Financial Account offsets CA and KA
    FA.append(new_FA)
    KA.append(KA[-1])  # Assume KA remains constant for simplicity

# Visualization
time = np.arange(time_steps)

plt.figure(figsize=(12, 6))

# Plot BoP components
plt.plot(time, CA, label="Current Account (CA)", marker='o')
plt.plot(time, KA, label="Capital Account (KA)", marker='o')
plt.plot(time, FA, label="Financial Account (FA)", marker='o')
plt.axhline(y=0, color='r', linestyle='--', label="Equilibrium Line (BoP=0)")

plt.title("Balance of Payments Dynamics")
plt.xlabel("Time Period")
plt.ylabel("Balance ($ billion)")
plt.legend()
plt.grid()
plt.tight_layout()
plt.show()

Explanation

1. Current Account Dynamics:
   • Adjustments in trade balance affect the CA over time.
2. BoP Equilibrium:
   • The KA and FA adjust to ensure that the BoP identity $(\text{CA} + \text{KA} + \text{FA} = 0)$ holds.
3. Assumptions:
   • The KA remains constant for simplicity.
   • Changes in the CA are offset entirely by the FA.

Results

• The graph shows how the Current Account evolves due to trade balance changes.
• The Financial Account responds dynamically to maintain BoP equilibrium.
• The Capital Account remains constant, highlighting its lesser role in this simplified example.

This example illustrates the balancing act within the BoP framework, showing the interdependence of its components.

Example Problem: Inflation Targeting Simulation

In this example, we aim to simulate the inflation dynamics under a basic monetary policy rule, such as the $Taylor Rule$.

The $Taylor Rule$ adjusts the nominal interest rate based on deviations of inflation and output from their target levels.

We’ll model how inflation converges to the target over time under this rule.

Problem

Given:

1. Inflation starts at a level significantly higher than the target.
2. A central bank adjusts the interest rate to influence inflation.
3. The inflation equation is simplified to a function of past inflation, output gap, and the interest rate.

We will simulate how inflation evolves over time with:

• A fixed inflation target.
• Parameters for monetary policy responsiveness.

Python Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53

import numpy as np
import matplotlib.pyplot as plt

# Parameters
inflation_target = 2.0  # Target inflation (%)
initial_inflation = 5.0  # Initial inflation (%)
output_gap = 1.0  # Initial output gap (%)
beta_inflation = 1.5  # Responsiveness to inflation deviation
beta_output = 0.5  # Responsiveness to output gap
time_steps = 20  # Number of time periods
natural_rate_of_interest = 2.0  # Natural rate of interest

# Simulation variables
inflation = [initial_inflation]
interest_rate = []

# Simulate inflation dynamics
for t in range(time_steps):
    # Calculate interest rate using Taylor Rule
    interest = natural_rate_of_interest + beta_inflation * (inflation[-1] - inflation_target) + beta_output * output_gap
    interest_rate.append(interest)
    
    # Update inflation (simplified dynamics: next inflation is dampened by interest rate)
    next_inflation = inflation[-1] - 0.2 * (interest - natural_rate_of_interest)
    inflation.append(next_inflation)

# Visualization
time = np.arange(0, time_steps + 1)

plt.figure(figsize=(12, 6))

# Plot inflation dynamics
plt.subplot(2, 1, 1)
plt.plot(time, inflation, marker='o', label="Inflation")
plt.axhline(y=inflation_target, color='r', linestyle='--', label="Inflation Target")
plt.title("Inflation Dynamics under the Taylor Rule")
plt.xlabel("Time Period")
plt.ylabel("Inflation (%)")
plt.legend()
plt.grid()

# Plot interest rate
plt.subplot(2, 1, 2)
plt.plot(time[:-1], interest_rate, marker='o', color='orange', label="Interest Rate")
plt.axhline(y=natural_rate_of_interest, color='r', linestyle='--', label="Natural Rate of Interest")
plt.title("Interest Rate Adjustment")
plt.xlabel("Time Period")
plt.ylabel("Interest Rate (%)")
plt.legend()
plt.grid()

plt.tight_layout()
plt.show()

Explanation

1. Taylor Rule:
   • Sets the nominal interest rate based on inflation deviation and output gap.
2. Inflation Dynamics:
   • Inflation adjusts gradually in response to interest rate changes.
3. Visualization:
   • The top graph shows how inflation converges to the target.
   • The bottom graph displays the interest rate adjustments over time.

Results

• Inflation converges toward the target ($2$%) over time due to the central bank’s interest rate adjustments.
• Interest rates stabilize as inflation reaches the target.
  This illustrates the effectiveness of inflation targeting under the $Taylor Rule$.

Problem Statement: Debt Dynamics with Stochastic Shocks

Objective:
Analyze the debt-to-GDP ratio of a country over time under stochastic (random) variations in interest rates $(r_t)$ and GDP growth rates $(g_t)$.

This captures the uncertainty in economic conditions and provides insights into debt sustainability in a volatile environment.

Model

The debt-to-GDP ratio evolves as:
$$
b_{t+1} = b_t \cdot \frac{1 + r_t}{1 + g_t} - p
$$

• $b_t$: Debt-to-GDP ratio at time $t$.
• $r_t$: Random interest rate at $t$, drawn from a normal distribution.
• $g_t$: Random GDP growth rate at $t$, drawn from a normal distribution.
• $p$: Primary balance as a percentage of GDP ($p > 0$ for surplus, $p < 0$ for deficit).

The randomness in $r_t$ and $g_t$ represents economic uncertainties, such as financial market shocks or unexpected GDP contractions.

Python Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46

import numpy as np
import matplotlib.pyplot as plt

# Parameters
T = 50                  # Time horizon (years)
b0 = 0.8                # Initial debt-to-GDP ratio
p = -0.02               # Primary deficit (-2% of GDP)

# Stochastic parameters
mean_r = 0.03           # Mean interest rate (3%)
std_r = 0.01            # Standard deviation of interest rate
mean_g = 0.02           # Mean GDP growth rate (2%)
std_g = 0.015           # Standard deviation of GDP growth rate

# Simulation
np.random.seed(42)  # For reproducibility
r_t = np.random.normal(mean_r, std_r, T)  # Random interest rates
g_t = np.random.normal(mean_g, std_g, T)  # Random GDP growth rates

time = np.arange(T)
debt_to_gdp = np.zeros(T)
debt_to_gdp[0] = b0

for t in range(1, T):
    debt_to_gdp[t] = debt_to_gdp[t-1] * ((1 + r_t[t]) / (1 + g_t[t])) - p

# Visualization
plt.figure(figsize=(12, 6))

# Debt-to-GDP Ratio Over Time
plt.plot(time, debt_to_gdp, label="Debt-to-GDP Ratio", color="blue")
plt.axhline(1, color="black", linestyle="--", alpha=0.7, label="100% Debt-to-GDP")
plt.fill_between(time, 0, debt_to_gdp, color="blue", alpha=0.2)
plt.title("Debt-to-GDP Dynamics with Stochastic Shocks")
plt.xlabel("Time (Years)")
plt.ylabel("Debt-to-GDP Ratio")
plt.legend()
plt.grid()
plt.tight_layout()
plt.show()

# Summary
print(f"Initial Debt-to-GDP Ratio: {b0:.2f}")
print(f"Debt-to-GDP Ratio After {T} Years: {debt_to_gdp[-1]:.2f}")
print(f"Average Interest Rate: {np.mean(r_t):.2%}")
print(f"Average GDP Growth Rate: {np.mean(g_t):.2%}")

Explanation of Code

1. Stochastic Variables:

   • The interest rate $(r_t)$ and GDP growth rate $(g_t)$ are drawn from normal distributions with specified means and standard deviations.
   • This randomness models real-world economic uncertainties.

2. Debt Evolution:

   • At each time step, the debt-to-GDP ratio evolves based on the stochastic parameters.

3. Visualization:

   • The debt-to-GDP trajectory shows the impact of economic volatility.
   • The shaded region emphasizes the fluctuation over time.

Results

Initial Debt-to-GDP Ratio: 0.80
Debt-to-GDP Ratio After 50 Years: 2.31
Average Interest Rate: 2.77%
Average GDP Growth Rate: 2.03%

1. Debt Path:

   • The debt-to-GDP ratio fluctuates due to variations in interest rates and GDP growth rates.
   • Persistent high interest rates or low growth can lead to unsustainable debt.

2. Long-Term Trend:

   • The final debt-to-GDP ratio depends on the cumulative effect of the stochastic shocks.

3. Policy Implications:

   • Governments should account for uncertainty when setting fiscal targets.
   • Debt sustainability requires managing both expected trends and risks.

Problem Statement: Business Cycle Analysis with Real GDP Data

Objective:
Analyze the business cycle by identifying and visualizing the cyclical components of Real GDP using the $Hodrick$-$Prescott (HP) filter$.

This method separates the trend and cyclical components of GDP, allowing us to examine deviations from the long-term trend.

Steps:

1. Data:
   Use simulated or publicly available Real GDP time series data.

2. Hodrick-Prescott Filter:

   • Decompose GDP into its trend $(T_t)$ and cyclical $(C_t)$ components:
     $$
     GDP_t = T_t + C_t
     $$
   • The $HP filter$ minimizes the following loss function:
     $$
     \sum_t (GDP_t - T_t)^2 + \lambda \sum_t \left[(T_{t+1} - T_t) - (T_t - T_{t-1})\right]^2
     $$
     where $(\lambda)$ is a smoothing parameter (typically $1600$ for quarterly data).

3. Visualize:
   Plot the original GDP, trend, and cyclical components to analyze economic fluctuations.

Python Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45

import numpy as np
import matplotlib.pyplot as plt
from statsmodels.tsa.filters.hp_filter import hpfilter

# Simulated GDP Data (quarterly, 50 periods)
np.random.seed(42)
time = np.arange(50)
trend = 2.5 * time + 100  # Linear trend
cyclical = 15 * np.sin(0.3 * time)  # Cyclical component
noise = np.random.normal(0, 5, size=time.shape)  # Random noise
gdp = trend + cyclical + noise  # Simulated GDP

# Apply the HP filter
gdp_trend, gdp_cycle = hpfilter(gdp, lamb=1600)

# Plotting
plt.figure(figsize=(14, 8))

# Original and Trend
plt.subplot(2, 1, 1)
plt.plot(time, gdp, label="Original GDP", color="blue")
plt.plot(time, gdp_trend, label="Trend (HP Filter)", color="red", linestyle="--")
plt.title("Real GDP and Trend Component")
plt.xlabel("Time (Quarters)")
plt.ylabel("GDP")
plt.legend()
plt.grid()

# Cyclical Component
plt.subplot(2, 1, 2)
plt.plot(time, gdp_cycle, label="Cyclical Component", color="green")
plt.axhline(0, color="black", linestyle="--", linewidth=0.8)
plt.title("Cyclical Component of GDP (Deviations from Trend)")
plt.xlabel("Time (Quarters)")
plt.ylabel("GDP Deviations")
plt.legend()
plt.grid()

plt.tight_layout()
plt.show()

# Print summary statistics
print(f"Original GDP Mean: {np.mean(gdp):.2f}")
print(f"Trend Component Mean: {np.mean(gdp_trend):.2f}")
print(f"Cyclical Component Mean: {np.mean(gdp_cycle):.2f} (should be ~0)")

Explanation of the Code

1. Simulated GDP Data:

   • Real GDP is modeled as a combination of a long-term trend, a cyclical fluctuation, and random noise.

2. Hodrick-Prescott Filter:

   • The $HP filter$ separates the GDP into a smooth trend and short-term deviations (cyclical component).
   • The smoothing parameter $(\lambda)$ controls the smoothness of the trend; $1600$ is standard for quarterly data.

3. Visualization:

   • The first plot shows the original GDP and its trend.
   • The second plot highlights the cyclical component, indicating deviations from the long-term trend.

Results

Original GDP Mean: 161.77
Trend Component Mean: 0.00
Cyclical Component Mean: 161.77 (should be ~0)

1. Trend Component:

   • Represents the long-term economic growth path.

2. Cyclical Component:

   • Indicates business cycle fluctuations around the trend.
   • Peaks and troughs correspond to economic expansions and contractions, respectively.

3. Insights:

   • The cyclical component helps identify recessions and booms.
   • Policymakers and economists use this analysis to design countercyclical measures.

Problem Statement: Price Discrimination in a Monopoly

Objective:
A monopolist sells its product to two distinct markets with different demand elasticities.
The monopolist can price discriminate, setting different prices for each market to maximize overall profit.

The goal is to:

1. Simulate the optimal pricing strategy for the monopolist.
2. Compute the quantities sold and profits in each market.
3. Visualize the results to illustrate price discrimination.

Assumptions

1. Demand Functions:

   • Market $1$: $( Q_1 = a_1 - b_1 \cdot P_1 )$
   • Market $2$: $( Q_2 = a_2 - b_2 \cdot P_2 )$
     $ P_1 $ and $ P_2 $ are prices set in each market.

2. Profit Function:
   The monopolist’s total profit is:
   $$
   \pi = (P_1 - c) \cdot Q_1 + (P_2 - c) \cdot Q_2
   $$
   where $c$ is the marginal cost.

3. Goal:
   Maximize $\pi$ by choosing optimal $P_1$ and $P_2$.

Python Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66

import numpy as np
from scipy.optimize import minimize
import matplotlib.pyplot as plt

# Parameters
a1, b1 = 100, 2  # Market 1: Demand intercept and slope
a2, b2 = 120, 4  # Market 2: Demand intercept and slope
c = 20           # Marginal cost

# Demand functions
def demand1(p1):
    return max(0, a1 - b1 * p1)

def demand2(p2):
    return max(0, a2 - b2 * p2)

# Profit function
def total_profit(prices):
    p1, p2 = prices
    q1 = demand1(p1)
    q2 = demand2(p2)
    return -((p1 - c) * q1 + (p2 - c) * q2)  # Negative for minimization

# Initial guesses and bounds
initial_prices = [50, 50]
bounds = [(c, a1 / b1), (c, a2 / b2)]  # Prices must be at least marginal cost

# Solve for optimal prices
result = minimize(total_profit, initial_prices, bounds=bounds)
p1_opt, p2_opt = result.x
q1_opt = demand1(p1_opt)
q2_opt = demand2(p2_opt)
profit_opt = -(result.fun)

# Visualization
prices = np.linspace(0, 80, 500)
revenues1 = [(p - c) * demand1(p) for p in prices]
revenues2 = [(p - c) * demand2(p) for p in prices]

plt.figure(figsize=(12, 6))

# Revenue curves
plt.plot(prices, revenues1, label="Market 1 Profit", color="blue")
plt.plot(prices, revenues2, label="Market 2 Profit", color="orange")

# Optimal points
plt.scatter(p1_opt, (p1_opt - c) * q1_opt, color="blue", label=f"Opt. Price Market 1: ${p1_opt:.2f}")
plt.scatter(p2_opt, (p2_opt - c) * q2_opt, color="orange", label=f"Opt. Price Market 2: ${p2_opt:.2f}")

plt.axvline(p1_opt, color="blue", linestyle="--", alpha=0.7)
plt.axvline(p2_opt, color="orange", linestyle="--", alpha=0.7)

plt.title("Profit Maximization with Price Discrimination")
plt.xlabel("Price")
plt.ylabel("Profit")
plt.legend()
plt.grid()
plt.tight_layout()
plt.show()

# Print results
print(f"Optimal Price in Market 1: ${p1_opt:.2f}")
print(f"Optimal Price in Market 2: ${p2_opt:.2f}")
print(f"Quantity Sold in Market 1: {q1_opt:.2f}")
print(f"Quantity Sold in Market 2: {q2_opt:.2f}")
print(f"Total Profit: ${profit_opt:.2f}")

Explanation of Code

1. Demand Functions:

   • Market $1$ has a lower price sensitivity ($b_1 < b_2$), indicating it is less elastic.
   • Market $2$ is more elastic ($b_2 > b_1$), meaning customers are more price-sensitive.

2. Profit Maximization:

   • The monopolist maximizes profit by finding the best prices for each market using scipy.optimize.minimize.
   • The constraints ensure prices are at least equal to the marginal cost.

3. Visualization:

   • The revenue curves for both markets show how profit changes with price.
   • The optimal prices for both markets are highlighted.

Results

ptimal Price in Market $1$: $35.00
Optimal Price in Market $2$: $25.00
Quantity Sold in Market $1$: 30.00
Quantity Sold in Market $2$: 20.00
Total Profit: $550.00

1. Optimal Prices:

   • The monopolist sets a higher price in the less elastic market (Market $1$).
   • The price is lower in the more elastic market (Market $2$).

2. Quantities and Profit:

   • The quantity sold is higher in Market $2$ due to its larger demand intercept.
   • The monopolist achieves maximum profit by price discrimination.

3. Graph:

   • The revenue curves illustrate the profit-maximizing prices visually.
   • The points of tangency indicate the optimal strategy.