A new numerical model is proposed herein that predicts the space-time profiles of an electromagnetic wave in an anisotropic warm plasma. The model is based upon a second-order spatial and temporal discretization of the coupled equations of Maxwell and Euler. In order to comply with the basic nature of the curl, divergence, and gradient operators, an uncollocated nodal scheme is used to approximate the electric field, magnetic field, velocity, and pressure throughout the computational domain. The theory is postulated in a three-dimensional domain but is verified in two dimensions. A magnetic line source is used to excite the medium, and the numerical data are tested against closed-form solutions. As demonstrated herein, the agreement is shown to be excellent.