Simulazione campi in g.o. metallica rettangolare MATLAB

[deltax]

Salve, dovrei simulare i campi in una guida d'onda metallica rettangolare, al variare dei parametri m,n dei modi (devo simulare sia il caso TE che il caso TM) dove b=a/2. Dovrei plottarli come nella seguente figura:
http://it.wikipedia.org/wiki/Guida_d%27onda
Ho incontrato una serie di problemi, alcuni pratici, alcuni teorici. Tra quelli teorici ci sono: avendo l'espressione di Az e Fz, Le espressioni ricavate di Ex,Ey,Ex,Hx,Hy,Hz sono funzioni di (x,y,z) mentre nella figura si ha una figura bidimensionale. Tra I problemi pratici per esempio quello di non riuscire a fare la derivata seconda mista.
Qualcuno saprebbe aiutarmi?