simulateNHPP function

Simulates event times from an inhomogenous Poisson process on [0,T]