A numerical method based on the the method of characteristics for hyperbolic systems of partial differential equations in four independent variables is developed and used for solving time domain Maxwell's equations. Th ...