FDTD Simulation in Matlab,C++ or FORTRAN

Does anyone know how to code an FDTD code in Matlab for a dielectric slab? I don't mind being pointed in the helpful resources so that I do my own reading. The slab is like three horizontal slabs (like a sandwidch) and the EM waves propagate in the middle slab

Your starting points should be two textbooks, one by Taflove and Hagness and another by Sullivan, which are full of hints.

Happy reading!

Hi, I went through the books and managed to write simple codes in C and MATLAB. But my problem is: all the simulations have the dielectric vertical. I would like to have the same scenario, but with waves propagation through the dielectric that is bounded by layers that have a much higher conductivity. How do I 'rotate' the computational domain (in 1D first and the I can do the rest)?