求助mathematica程序问题,那个flux代表了什么?module那块是做什么用的?
Clear["Global`*"]
tu :=
E^(I \[CurlyPhi]) v1^2 ((0.` + 0.` I) +
v0 ((0.` + 16.` I) e^3 +
e^2 ((0.` - 24.` I) \[CurlyEpsilon]l - (0.` +
24.` I) \[CurlyEpsilon]u) + \[CurlyEpsilon]l ((0.` -
8.` I) \[CurlyEpsilon]l - (0.` +
8.` I) \[CurlyEpsilon]u) \[CurlyEpsilon]u + ((0.` +
0.5832` I) \[CurlyEpsilon]l + (0.` +
0.5832` I) \[CurlyEpsilon]u) \[CurlyPhi]^2 +
e ((0.` + 8.` I) \[CurlyEpsilon]l^2 + (0.` +
32.` I) \[CurlyEpsilon]l \[CurlyEpsilon]u + (0.` +
8.` I) \[CurlyEpsilon]u^2 - (0.` +
1.1664` I) \[CurlyPhi]^2)) Cos[\[CurlyPhi]/
2] Sin[\[Theta]])
Tu[en_,\[Phi]n_] := Module[{},
e := en;
\[CurlyPhi] := 2*Pi*\[Phi]n/4;
\[Theta] := ArcCos[- (((e - \[CurlyEpsilon]0)/v0)/2.0)];
t1 := tu;
Return[(t1)]
]
flux := .0
e := 0.
\[CurlyEpsilon]0 = 0.0
v0 = 1
\[CurlyEpsilon]u = 0.3
\[CurlyEpsilon]l = 0.3
v1 = 0.07
t1u[en_,\[Phi]n_] := Abs[Tu[en,\[Phi]n]]^2;
Plot[(Abs[Tu[e,flux]]^2),{e,-2,2},
PlotRange -> {{-2,2.},{0,1.1}},Axes -> False,Frame -> True,
FrameLabel -> {"\[CapitalPhi]/\[CapitalPhi]o","Tu"},
BaseStyle -> {"Times New Roman",14},
PlotStyle -> {AbsoluteThickness[2.8]}]