Use the methods of superposition and thinning to write an algorithm for simulating a non-homogeneous Poisson process with arrival rate (t) = 0.8+0.4sin(0.2t) over the time interval (0, 600]. Do this by dividing each period into 5 equal subintervals, compute the maximum and minimum values for for each of these 5 regions, and make a table of these two values for every region.

Preferably in Python, but other languagues welcome.