We investigate numerical techniques to locate the event points of the differential system x′ = f(x), characterized by a vector field f which is a discontinuous function along an event surface Σ = {x ∈ R| h(x) = 0} splitting the state space into two different regions R1 and R2, particularly f(x) = fi(x) when x ∈ Ri, for i = 1, 2 and f1(x) 6= f2(x) when x ∈ Σ. We propose event location techniques...