Determining Liouvillian First Integrals for Dynamical Systems in the Plane

Here we present/implement an algorithm to find Liouvillian first integrals of dynamical systems in the plane. In \cite{JCAM}, we have introduced the basis for the present implementation. The particular form of such systems allows reducing it to a single rational first order ordinary differential equation (rational first order ODE). We present a set of software routines in Maple 10 for solving rational first order ODEs. The package present commands permitting research incursions of some algebraic properties of the system that is being studied.