We develop a numerical formulation to calculate the classical motion of charges in strong electromagnetic fields, such as those occurring in high-intensity laser beams. By reformulating the dynamics in terms of SL(2, C) matrices representing the Lorentz group, our for-