Pop_fun function

Function that solves the population dynamics equations including discontinuities