We define two important arrays E zeros(max_space, 1) %Initialize Electric array H = E iNitialize Magnetic array Two things are important in understanding the algorithm compared to the equations. The first is that the values of array e are updated at time n+ 1/2 whereas the values of array H are updated at time n+ 1 The other factor in the time equation is the placement of the pulse. It happens between the updates of array E and array H. So the pulse is factored into the arrays H and e at time step n and n+1/2 respectively. A visualization of this process is given in the figure below E(k+1)at n Array h outer Pulse added to array e loop n Trav E E(k′+1) from n-1 kk+ k+1 k E(k)H(j)E(k+1)H(j+1)9 Two things are important in understanding the algorithm compared to the equations. The first is that the values of array E are updated at time n+1/2 whereas the values of array H are updated at time n + 1: The other factor in the time equation is the placement of the pulse. It happens between the updates of array E and array H. So the pulse is factored into the arrays H and E at time step n and n+1/2 respectively. A visualization of this process is given in the figure below