Here we present a new approach to compute symmetries of rational second order ordinary differential equations (rational 2ODEs). This method can compute Lie symmetries (point symmetries, dynamical symmetries and nonlocal symmetries) admitted by a rational 2ODE if it presents a Liouvillian first integral. The procedure is based on an idea arising from the formal equivalence between the total derivative operator and the vector field associated with the 2ODE over its solutions (Cartan vector field). Basically, from the formal representation of a Lie symmetry it is possible to extract information that allows to use this symmetry practically (in the 2ODE integration process) even in cases where the formal operation cannot be performed, i.e., in cases where the symmetry is nonlocal. Furthermore, when the 2ODE in question depends on parameters, the procedure allows an analysis that determines the regions of the parameter space in which the integrable cases are located. We finish by applying the method to the 2ODE that represents a forced modified Duffing–Van der Pol oscillator (which presents chaos for arbitrary values of the parameters) and determining a region of the parameter space that leads to integrable cases, some not previously known.