matlab菲涅尔衍射光强分布,求助！！！用颜色表示的菲涅尔衍射的光强分布图

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

a = 0.001

x0 = 0

y0 = 0

z0 = -0.5

z = 0.1

\[Lambda] = 5*10^-7

k = 2*10^6

A = 1000

f[\[Xi]_, \[Eta]_, x_,

y_] := ((-(x0/Sqrt[z0^2 + x0^2 + y0^2]) - x/Sqrt[

z^2 + x^2 + y^2]) \[Xi] + (-(y0/Sqrt[z0^2 + x0^2 + y0^2]) - y/

Sqrt[z^2 + x^2 + y^2]) \[Eta] +

1/2*((1/Sqrt[z0^2 + x0^2 + y0^2] + 1/Sqrt[

z^2 + x^2 + y^2])*(\[Xi]^2 + \[Eta]^2) - (-(x0/Sqrt[

z0^2 + x0^2 + y0^2])*\[Xi] + -(y0/Sqrt[

z0^2 + x0^2 + y0^2])*\[Eta])^2/Sqrt[

z0^2 + x0^2 +

y0^2] - (x/Sqrt[z^2 + x^2 + y^2]*\[Xi] +

y/Sqrt[z^2 + x^2 + y^2]*\[Eta])^2/Sqrt[

z^2 + x^2 + y^2]))*(-I*A)/\[Lambda]*

Exp[I*k (Sqrt[z0^2 + x0^2 + y0^2] + Sqrt[z^2 + x^2 + y^2])]/(

Sqrt[z0^2 + x0^2 + y0^2] + Sqrt[z^2 + x^2 + y^2])

DensityPlot[

Re[Integrate[

Exp[I*k*f[\[Xi], \[Eta], x]], {\[Xi], -a, a}, {\[Eta], -a,

a}]]^2, {x, 0, 0.1}, {y, 0, 0.1},

ColorFunction -> "BlueGreenYellow"]

这是我的代码，但是一直报错，求大佬看看有什么问题