  A set of Mathematica functions is presented to model classically two of the
most important processes in strong field physics, namely high-order harmonic
generation (HHG) and above-threshold ionization (ATI). Our approach is based on
the numerical solution of the Newton-Lorentz equation of an electron moving on
an electric field and takes advantage of the symbolic languages features and
graphical power of Mathematica. Similarly as in the Strong Field Approximation
(SFA), the effects of atomic potential on the motion of electron in the laser
field are neglected. The SFA has proven to be an essential tool in strong field
physics in the sense that it is able to predict with great precision the
harmonic (in the HHG) and energy (in the ATI) limits. We have extended
substantially the conventional classical simulations, where the electric field
is only dependent on time, including spatial nonhomogeneous fields and spatial
and temporal synthesized fields. Spatial nonhomogeneous fields appear when
metal nanosystems interact with strong and short laser pulses and temporal
synthesized fields are routinely generated in attosecond laboratories around
the world. Temporal and spatial synthesized fields have received special
attention nowadays because they would allow to exceed considerably the
conventional harmonic and electron energy frontiers. Classical simulations are
an invaluable tool to explore exhaustively the parameters domain at a cheap
computational cost, before massive quantum mechanical calculations, absolutely
indispensable for detailed analysis, are performed.
