By defining the Hamiltonian function in an electron optical system the canonical aberration theory has been developed in the up to fifth‐order approximation. Therefore, the canonical position and momentum aberrations can be expressed by eikonal functions and equivalently by Poisson brackets. Consequently, all fifth‐order canonical aberration coefficients for rotationally symmetrical electron optical systems have been derived at any arbitrary observation planes.