About calculation of the EFS of Photonic Crystals

Anyone have the matlab code for equifrequency contour (EFS) of photonic crystals? Would you please paste it here?
Thanks very much!

The code is no different than the standard PWEM code. Search the discussion groups for code from Guo (or Gou?). From there, rather than calculate the eigen-frequencies just around the perimeter of the Brillouin zone, you will calculate the eigen-frequencie inside the entire Brillouin zone. The isofrequency contours in this space are your EFSs.