Dva bodové zdroje

DensityPlot[(1/(1000+(x^2)+(y^2)))((Cos[15(x/(Sqrt[1000+(x^2)+(y^2)]))])^2),{x,-75,75},{y,-75,75},PlotRange->All,PlotPoints->100,ImageSize->600]

Obdélníkový otvor

DensityPlot[(1/(1000+(x^2)+(y^2)))(((((Sin[((a x)/(Sqrt[1000+(x^2)+(y^2)]))])/((a x)/(Sqrt[1000+(x^2)+(y^2)])))))^2)(((((Sin[((b y)/(Sqrt[1000+(x^2)+(y^2)]))])/((b y)/(Sqrt[1000+(x^2)+(y^2)])))))^2)/.{a->80,b->25},{x,-15,15},{y,-15,15},PlotRange->{0,0.000025},PlotPoints->100,ImageSize->600]

Difrakční mřížka

DensityPlot[(1/(1000+(x^2)+(y^2)))(((((Sin[((a n x)/(Sqrt[1000+(x^2)+(y^2)]))])/(Sin[((a x)/(Sqrt[1000+(x^2)+(y^2)]))]))))^2)/.{a->20,n->5},{x,-15,15},{y,-10,10},PlotRange->All,PlotPoints->200,ImageSize->600,AspectRatio->2/3]

Kruhový otvor

GraphicsRow[{
DensityPlot[Evaluate[(1/(1000+(x^2)+(y^2)))(((((BesselJ[1,((a(Sqrt[(x^2)+(y^2)]))/(Sqrt[1000+(x^2)+(y^2)]))])/((a(Sqrt[(x^2)+(y^2)]))/(Sqrt[1000+(x^2)+(y^2)])))))^2)/.{a->25,n->5}],{x,-15,15},{y,-10,10},PlotPoints->200,PlotRange->All,ImageSize->600,AspectRatio->2/3,PlotLabel->"Difrakce na kruhovém otvoru - dobře patrné centrální maximum"],
DensityPlot[Evaluate[(1/(1000+(x^2)+(y^2)))(((((BesselJ[1,((a(Sqrt[(x^2)+(y^2)]))/(Sqrt[1000+(x^2)+(y^2)]))])/((a(Sqrt[(x^2)+(y^2)]))/(Sqrt[1000+(x^2)+(y^2)])))))^2)/.{a->25,n->5}],{x,-15,15},{y,-10,10},PlotPoints->200,ImageSize->600,AspectRatio->2/3,PlotLabel->"Difrakce na kruhovém otvoru - zvýrazněná maxima vyššího řádu"]
}]