In[]:=
ClearAll["Global`*"];​​ϵ=0.5;​​β=45;​​λ=1;​​F=2000;​​Nu=
829440
103680-8640
2
ϵ
+3960
4
ϵ
+297
6
ϵ
+40
2
F
6
ϵ
2
λ
+40
2
F
6
ϵ
2
λ
Cos[2β]
;​​T1=Nu
1
4
(1-
2
r
)+ϵ
1
8
(
3
r
-r)Sin[ζ]+
2
ϵ
1
24
Cos[2ζ](
4
r
-
2
r
)+
2
r
16
-
4
r
32
-
1
32
;​​ContourPlot[T1,{r,-1,1},{ζ,0,2Pi},ContourLabelsTrue,PlotTheme->"Scientific",PlotPoints10]
Out[]=