Use a Matplotlib gridspec to plot the Fraunhofer pattern of a 2 rectangle aperture.

The used solution is

$$ u_z(x_z, y_z)\propto \exp\left(j\frac{\pi}{\lambda Z}(x_z^2+y_z^2)\right) 4hw \mbox{ sinc}\left(\dfrac{2 x_z w}{\lambda Z} \right) \cos\left(\dfrac{2\pi x_z d}{\lambda Z}\right) \mbox{sinc}\left(\dfrac{y_z h}{\lambda Z} \right) $$

where $x_z$ and $y_z$ are the coordinates on the screen, $Z$ is the distance from the object to the screen and $\lambda$ is the wavelength of the light.

The plotted pattern will be given by

$$ I_z (x_z, y_z) \propto \left| u_z(x_z, y_z) \right|^2 $$

When executed, the Python 3 script Fraunhofer-difraction.py will produce the following output:

The parameters used in this case are:

$$ \begin{array}{rcr} Z&=&1\mbox{ m},\\ w&=&0.1\times 10^{-3}\mbox{ m},\\ h&=&5\times 10^{-3}\mbox{ m},\\ d&=&0.7\times 10^{-3}\mbox{ m},\\ \lambda&=&500\times 10^{-9}\mbox{ m}.\\ \end{array} $$