for n =2: max time-1 iNner Loop E -Increments electric wave in space for k = 2: max_space E(k) E(k)+eta*(H(k-1)-H(k)); end %Hard Source- imposes a value on the grid pulse exp(-5*((to-n)/spread)"2); E(center__problem_space)= pulse; %Inner Loop H Increments magnetic wave in space for J =1: max_space-1 H(j)=H(j)+eta*(E(j)-E(j+1)); end end E2(k)=En2(k)-0.[”(k+1/2)-(k-12) 分(k+1/2)=B(+12)-05E(+1)-E2(k)( ) ( ) ( ) ( ) 1/2 1/2 0.5 1/ 2 1/ 2 n n n n x x y y k k k k + − = − + − −     E E H H ( ) ( ) ( ) ( ) 1 1/2 1/2 1/ 2 1/ 2 0.5 1 n n n n y y x x k k k k + + + + = + − + −   H H E E  