How to use event solver to compute numerical solution to impulsively forced van der pol oscillator?

I want to use the event solver to find numerical solution to the van der pol equation with impulsive forcing on the right hand side.

I was told that I could use event solver to incorporate impulsive force on the system by first integrating using ode45 and by stopping the integration at time t when i want impulse force to exert and looping it over for multiple n.

How should I write a code on matlab if I want it to compute this for me? Should I start with unforced van der pol system and add the integration step then use the event solver? Or is there any other way to incorporate dirac delta function so that there won't be any problems in the integration process?