<< Calculus`VectorAnalysis`
w = 10; a = 1; b = 0.5; H = 1; \
kc = Sqrt[\[Pi]^2*（m^2/a^2 + n^2/b^2 ）]; kz = (w^2 - kc^2)^(1/2);
Hz[x_, y_, z_, t_] :=
H*Cos[m*\[Pi]*x/a]*Cos[n*\[Pi]*y/b]*Cos[kz*z - w*t];
Hx[x_, y_, z_, t_] := -kz/kc^2 *H*m*\[Pi]/a*Sin[m*\[Pi]*x/a]*
Cos[n*\[Pi]*y/b] *
Sin[kz*z -
w*t] ; \
Hy[x_, y_, z_, t_] := -kz/kc^2*H*n*\[Pi]/b*Cos[m*\[Pi]*x/a]*
Sin[n*\[Pi]*y/b]*
Sin[kz*z -
w*t]; \
\
Ex[x_, y_, z_, t_] := -kz/w *
CrossProduct[{0, 0, 1}, {Hx[x, y, z, t], Hy[x, y, z, t], 0}][[
1]];
Ey[x_, y_, z_, t_] := -kz/w *
CrossProduct[{0, 0, 1}, {Hx[x, y, z, t], Hy[x, y, z, t], 0}][[
2]];
Ez[x_, y_, z_, t_] := 0
Ea = Sqrt[Ex[x, y, z, t]^2 + Ey[x, y, z, t]^2] ;
H = Sqrt[
Hz[x, y, z, t]^2 + Hy[x, y, z, t]^2 + Hz[x, y, z, t]^2] ; p1 =
ContourPlot[
Ea /. {m\[LongRightArrow]1, n\[LongRightArrow]0, t -> 0,
z -> 2}, {x, 0, a} , {y, 0, b}, PlotPoints -> 30,
AspectRatio -> Automatic,
DisplayFunction ->
Identity] \
p2 = Plot3D[
Ea /. {m\[LongRightArrow]1, n\[LongRightArrow]0,
t\[LongRightArrow]0, z -> 2}, {x, 0, a} {y, 0, b} ,
PlotPoints\[LongRightArrow]30,
AspectRatio\[LongRightArrow]Automatic,
DisplayFunction\[LongRightArrow]Identity]

第一，把中文（）改成英文的，箭头->也要改；
第二，Plot3D里面少一个逗号；
第三，不需要<< Calculus`VectorAnalysis` 以及其中的CrossProduct，9版本后有现成的Cross。