First order ODEs: Mathematica and symbolic-numerical methods