fdtd matlab codes for rods in free space

hi
does anybody know how to write matlab codes fdtd for structure like square with on the corners there are rods with special dielctric constant but out of these rods is aire(free space)?
please help me. i really nedd these codes.

https://www.edaboard.com/viewtopic.p...c172a1ecee2869