Taylor series methods compute a solution to an initial value problem in ordinary differential equations by expanding each component of the solution in a long Taylor series. The series terms are generated recursively using the techniques of automatic differentiation. The ATOMFT system includes a translator to transform statements of the system of ODEs into a FORTRAN 77 object program that is compiled, linked with the ATOMFT runtime library, and run to solve the problem. We review the use of the ATOMFT system for nonstiff and stiff ODEs, the propagation of global errors, and applications to differential algebraic equations arising from certain control problems, to boundary value problems, to numerical quadrature, and to delay problems.