level 4
r[i_] := h + (2 i - 1)/(2 n)*L2 Sin[\[Alpha]] Cot[\[Theta]];
\[Alpha] =
10 Degree; \[CurlyPhi]1 = 0; \[CurlyPhi]2 = 0; L2 = 65; L1 = 65; \
Wmax = -1.6; h = 900; n = 10; \[Theta] = 85 Degree;
W[x_, y_] :=
Sum[Wmax/(2 (1 +
Tan[\[Alpha]] Cot[\[Theta]]) (Erf[
Sqrt[Pi]/
r[i] (y -
h Cot[\[Theta]] - ((i - 1)/
n) L2 Cos[\[Alpha]] Cos[\[CurlyPhi]1] (1 +
Tan[\[Alpha]] Cot[\[Theta]]))] -
Erf[Sqrt[Pi]/
r[i] (y -
h Cot[\[Theta]] - (i/n)*
L2 Cos[\[Alpha]] Cos[\[CurlyPhi]1] (1 +
Tan[\[Alpha]] Cot[\[Theta]]))]) ) ((Erf[
Sqrt[Pi] (x + (i - 1)/n*L2 Sin[\[CurlyPhi]1])/r[i]] -
Erf[Sqrt[
Pi] (x + (i - 1)/n*
L2 Sin[\[CurlyPhi]1] - (L1 +
L2 (2 i - 1)/(2 n)) (Sin[\[CurlyPhi]1] +
Cos[\[CurlyPhi]1] + Tan[\[CurlyPhi]2]))/r[i]]))/2, {i,
1, 100}]; Plot3D[W[x, y], {x, -1000, 2000}, {y, -1000, 1000},
PlotRange -> All, AxesLabel -> {"x", "y", "W[x, y]"},
PlotLabel -> "W[x, y]的三维图"]
Print["在 x = 0 且 y = 500 处的具体值为: ", result]
2023年12月08日 16点12分



