The algorithm and results of calculations of the relative energy density in an image point on a square area element are described. The energy distribution in the image point is calculated for aberration-free optical systems with circular and annular pupils in the Fraunhofer diffraction approximation. Longitudinal and transverse displacements of the center of the square area element relative to the paraxial image of the point are considered.