function?Uz?=?frenel_diff（T1z1lambdadim）
EM?=?size（T11）;
X?=?dim;
Y?=?dim;
x?=?linspace（-dimdimEM）;
y?=?linspace（-dimdimEM）;

k?=?2*pi/lambda;

[XY]?=?meshgrid（xy）;
u1?=