Infinite parallel-plates Simulation

Hi,
does anyone can help me? I need do make a matlab simulator for plot the field distribution (TE and TM modes) inside a parallel plates waveguide.

I've found the lines of field distribution equations (they are corect), but when I make it in matlab, the field does not correspond to the expected.

Please, I realy need to solve it.

Thanks.