fun[x_?NumberQ,
y_?NumberQ, \[Epsilon]_] := (Sqrt[\[Epsilon]] -
1/Sqrt[\[Epsilon]])*2*
Tan[Sqrt[\[Epsilon]]*x] - (-I*Sqrt[y] + 1/(-I*Sqrt[y]))*
Tanh[-I*Sqrt[y]*x] - (\[Epsilon]/(-I*Sqrt[y]) -
I*Sqrt[y]/\[Epsilon])*(Tan[Sqrt[\[Epsilon]]*x])^2*
Tanh[-I*Sqrt[y]*x];epssol[x_?NumberQ, y_?NumberQ] := \[Epsilon] /.
NSolve[fun[x, y, \[Epsilon]] == 0, \[Epsilon]]；ContourPlot[epssol[x, y], {x, 0.02, 0.08}, {y, -9000, -1000}]
就是我想干这样一件事，，用x，y画一个二维的图，然后Epsilon是x，y的隐函数，然后画一张强度分布的图，用颜色深浅表示大小，，
最终目标是这样：

1，你的代码里有中文的分号，必须改为英文的分号；
2，你的epssol函数有误，等号右边完全是个莫名其妙的东西，恐怕是转录代码时忘了其它重要部分。