LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUkjbWlHRiQ2I1EhRic=
with(plots):
a:=0.4*10^(-4);
d:=1.80*10^(-4);
N:=2.0;
lambda:=530.0*10^(-9);
L:=13.0;
Pi*a/lambda/L;
Pi*d/lambda/L;
xmax:=0.20;
densityplot(exp(-y^2/(2*0.03^2))*(sin(Pi*a*x/(lambda*L))^2/(Pi*a*x/(lambda*L))^2)*(sin(Pi*N*d*x/(lambda*L))^2/sin(Pi*d*x/(lambda*L))^2),x=-xmax...xmax,y=-0.5...0.5,style=PATCHNOGRID,scaletorange=0..0.2,grid=[100,100]);
plot(sin(Pi*a*y/(lambda*L))^2/(Pi*a*y/(lambda*L))^2*(sin(Pi*N*d*y/(lambda*L))^2/sin(Pi*d*y/(lambda*L))^2),y=-xmax...xmax);