How to simulate the posterior π(θ|y) in Python

Question

I am trying to simulate a posterior using Python. The posterior should get updated as data being processed. Here is the set up of the question:

The likelihood distributio