import numpy as np
+import matplotlib.pyplot as plt
+import scipy.stats as stats
+
+# Time variable
+t = np.linspace(-3, 3, 500)
+
+# Normal distribution parameters
+mu = 0
+sigma = 1
+
+# Functions for the normal distribution
+f_t = stats.norm.pdf(t, mu, sigma) # PDF of the normal distribution
+S_t = 1 - stats.norm.cdf(t, mu, sigma) # Survival function (1 - CDF)
+h_t = f_t / S_t # Hazard function
+
+# Plotting the three panels
+fig, ax = plt.subplots(1, 3, figsize=(15, 5))
+
+# f(t)
+ax[0].plot(t, f_t, label='f(t)')
+ax[0].set_title('Probability Density Function (f(t))')
+ax[0].set_xlabel('Time (t)')
+ax[0].set_ylabel('f(t)')
+ax[0].legend()
+ax[0].grid(True, linestyle='--', linewidth=0.5)
+ax[0].spines['right'].set_visible(False)
+ax[0].spines['top'].set_visible(False)
+
+# S(t)
+ax[1].plot(t, S_t, label='S(t)', color='orange')
+ax[1].set_title('Survival Function (S(t))')
+ax[1].set_xlabel('Time (t)')
+ax[1].set_ylabel('S(t)')
+ax[1].legend()
+ax[1].grid(True, linestyle='--', linewidth=0.5)
+ax[1].spines['right'].set_visible(False)
+ax[1].spines['top'].set_visible(False)
+
+# h(t)
+ax[2].plot(t, h_t, label='h(t)', color='green')
+ax[2].set_title('Hazard Function (h(t))')
+ax[2].set_xlabel('Time (t)')
+ax[2].set_ylabel('h(t)')
+ax[2].legend()
+ax[2].grid(True, linestyle='--', linewidth=0.5)
+ax[2].spines['right'].set_visible(False)
+ax[2].spines['top'].set_visible(False)
+
+# Adjust layout
+plt.tight_layout()
+
+# Show plot
+plt.show()
+