10 January 2008 Corrected algebraic error in the expressions for Axy(E), Axy(W), Ayx(N) and Ayx(S). The earlier version was only correct in the special case when e=w and n=s. (I.e., when the mesh is uniform in the vicinity of the dielectric edges.) Many thanks to yuchunlu_china@hotmail.com, who pointed out this bug. The corrected expressions now agree with the isotropic equations of Lusse et al, in limit that the materials are isotropic. 8 July 2008 Corrected small error in external plotting functions: replaced "transp" with "transpose", which works for all versions of MATLAB.