WRITE ME A PYTHON CODE THAT DOES THIS

Using Python for Monte Carlo simulation for predicting protein nucleation aggregation A simple strategy for protein aggregation is given below:

a$)$ Simulate the above system using MC and the kinetic rates as inputs. Calculate the probability

distribution of $($histogram plots obtained from many runs of the simulation$)$ for the number of

stable nucleation complexes at different time points. $20+5=25$

b$)$ Find out the distribution of sizes of stable nucleation complexes. $7$

c$)$ Vary ke and kf to show i$)$ how the number of stable nucleation complexes vary with time ii$)$

how the maximum size of the stable nucleation complex vary with ke and kf$.4+4=8$