Airy disk derivation ($A(u)=2J_1(u)/u$)

Write $a=D/2$ and wavenumber $k=2\pi /\lambda$. The pupil Fourier transform, normalized by pupil area, gives the Fraunhofer diffraction wave amplitude. In pupil polar coordinates,The angular average is $J_0$ because $J_0(z)=(2\pi {)}^{-1}{\int }_0^{2\pi }{e}^{-iz\cos \varphi }d\varphi$. The Bessel function of the first kind identity $d\{z{J}_1(z)\}/dz=z{J}_0(z)$ then reduces the radial integral to $A(u)=2J_1(u)/u$. Taking the squared magnitude yields the Airy disk light intensity profile.^[ref]

The apparent division by zero is removable. The power seriesgives $A(u)=1-u^2/8+u^4/192+O(u^6)$ and $P(u)=1-u^2/4+5u^4/192+O(u^6)$. At the centre, every pupil contribution has the same wave phase, so their normalized sum is one. The first dark ring instead corresponds to $J_1(u)=0$ with $u>0$, rather than to the removable zero of $J_1$ at the origin.