how to make w-k diagram with Matlab

Hi, everyone, I have some problem.
I want make w-k diagram (with n) about wave propagating to infinite plane of glass, with Matlab. (like attchment file)

these days, I studied about it.. but still there is no result.

please. Is there Matlab code or site or manual for help?

did you checked on http://www.mathworks.com/matlabcentral/fileexchange/

i cheked there..
but coudn't find right code..
i found there are 2d, 3d PhC wave propagating code