How to compute lambda(t) for a self-exciting process in R using the thinning algorithm?

Question

Here is my code for the thinning algorithm for a self exciting process. I\'ve already finished the algorithm for a standard poisson process and a non-homogeneous.

I\'