from scipy.integrate import quad from scipy.stats import norm import numpy as np # Define the unnormalized posterior as a function of unnormalized_posterior = lambda : (1/np.sqrt(2*np.pi)) * np.exp(-0.5 * (1 - )**2) * ( 0.5/np.sqrt(2*np.pi) * np.exp(-0.5 * ( + 1)**2) + 0.5/np.sqrt(2*np.pi) * np.exp(-0.5 * ( - 1)**2) ) # Calculate the normalization constant normalization_constant, _ = quad(unnormalized_posterior, -np.inf, np.inf) # Print the normalization constant print("Normalization Constant:", normalization_constant)