First-order ODEs: Mathematica and Symbolic-Numerical Methods