Abstract: In this work we present computational methods for examining dynamical systems. We focus on those systems being characterized by slow–fast vector fields or corresponding differential algebraic equations that commonly occur in physical applications. In the latter ones scientists usually consider a manifold of admissible physical states and a vector field describing the time evolution of the physical system. The manifold is typically implicitly defined within a higher-dimensional space by a system of equations. Certain physical systems, such as relaxation oscillators, perform sudden jumps in their state evolution when they are forced into an unstable state. The main contribution of the present work is to model the dynamical evolution incorporating the jumping behavior from a perspective of computational geometry which not only provides a qualitative analysis but also produces quantitative results. We use geodesic polar coordinates (GPC) to numerically obtain explicit parametrizations of the implicitly defined manifold and of the relevant jump and hit sets. Moreover, to deal with the possibly high co-dimension of the considered implicitly defined manifold we sketch how GPC in combination with the cut locus concept can be used to numerically obtain an essentially injective global parametrization. This allows us to parametrize and visualize the dynamical evolution of the system including the aforementioned jump phenomena. As main tools we use homotopy approaches in conjunction with concepts from differential geometry. We discuss how to numerically realize and how to apply them to several examples from mechanics, electrical engineering and biology.