Suppose that I know all intensity functions lambda(t) during given period [0,t], how can I compute the mean function m(t) for non-homogeneous Poisson process?
Basically, m(t) in the integral of lambda(t) from 0 to t. Through simulation, I can actually compute all lambda(t). Is there an efficient algorithm to compute m(t)?