Double slit, Yee FDTD/Maxwell's equations in 2D (MATLAB)

for t = 1:512 for k = 1:N-1 for l = 2:N-1 ex(k, l) = ex(k, l) + (hz(k, l) - hz(k, l-1))/2; end end for k = 2:N-1 for l = 1:N-1 ey(k, l) = ey(k, l) + (hz(k-1, l) - hz(k, l))/2; end end for k = 1:N-1 for l = 1:N-1 hz(k, l) = hz(k, l) + (ex(k, l+1) - ex(k, l) + ey(k, l) - ey(k+1,l))/2; end end end