目的是基于mathematics仿真杨氏双缝干涉实验的干涉图样，以下是鄙人所写的代码：
(*定义参量，单位均为m*)
\[Lambda] = 632.8*10^-9; (*入射光波长*)
d = 1*10^-3; (*双缝间距*)
L = 1; (*双缝到屏的垂直距离*)
(*定义干涉函数T*)
Subscript[r, 1] = Sqrt[L^2 + (x - d/2)^2]; (*S1到观察点P的距离*)
Subscript[r, 2] = Sqrt[L^2 + (x + d/2)^2]; (*S2到观察点P的距离*)
\[Delta] =
Subscript[r, 2] - Subscript[r, 1] = (
Subscript[r, 2]^2 - Subscript[r, 1]^2)/(
Subscript[r, 2] + Subscript[r, 1]) = (2 xd)/(
Sqrt[L^2 + (x - d/2)^2] + Sqrt[L^2 + (x + d/2)^2]);(*两束光在P点处的光程差*)
T == 4 Subscript[T,
0]*(Cos[(
2 \[Pi]xd)/(\[Lambda] Sqrt[L^2 + (x - d/2)^2] + Sqrt[
L^2 + (x + d/2)^2])])^2;(*P点的光强的准确公式*)
(*绘制干涉条纹分布图样*)
DensityPlot[T, {x, -0.4, 0.4}, {y, -1, 1}]
问题是现在一直出不来图，错误显示标签被保护，图三是应该出来的效果图，拜托各位大佬有时间帮帮我🥹
我会怀有感恩的心，等学成归来后将这份善意传递出去🥹

\[Lambda] = 632.8*10^-9;
d = 1*10^-3;
L = 1;
T = 4 (Cos[(2 \[Pi] x d)/(\[Lambda] Sqrt[L^2 + (x - d/2)^2] +
Sqrt[L^2 + (x + d/2)^2])])^2;
DensityPlot[T, {x, -0.01, 0.01}, {y, -1, 1}, PlotRange -> Full]
你的πxd连在一起，作为了一个新的变量而不是变量的乘积π x d
T0没有定义

\[Lambda] = 632.8/10^9; d = 1/10^3; Subscript[T, 0] = 1; d = 1; L = 1;
T = 4 Subscript[T, 0]
Cos[(2 \[Pi] x d)/(\[Lambda] Sqrt[L^2 + (x - d/2)^2] + Sqrt[
L^2 + (x + d/2)^2])]^2;
DensityPlot[T, {x, -1, 1}, {y, -1, 1}]