Sqrt[x^2 + (y + d)^2]; node[x_, y_] = r1[x, y] - r2[x, y]; Export["Diffraction_beats.gif", Table[ContourPlot[{node[x, y] == 1 + a, node[x, y] == -1 + a,...
{10, 10}]; p0 = Table[ GraphicsRow[{ ArrayPlot[cd, PlotLabel -> "Diffraction limited image", LabelStyle -> {Black, Bold}], ArrayPlot[RotateRight[kernel...
