Compute the PSFs for each wavelength

$$PF \left( \lambda \right) = AT \times \exp \left( \frac{i2 \pi W}{\lambda} \right)$$ (1)

$$ASF \left( \lambda \right) = ifft\left( PF \left( \lambda \right)\right)$$ (2)

$$OTF \left( \lambda \right) = fft\left( PSF \left( \lambda \right)\right)$$ (3)

$$RI \left( \lambda \right) = I\left( \lambda \right) \otimes PSF\left( \lambda \right)$$ (4)

$$PSF \left( \lambda \right) = Mod \left( ASF \left( \lambda \right... ...F \left( \lambda \right) \times conj \left( ASF \left( \lambda \right) \right)$$ (5)

$$PTF\left( \lambda \right) = Phase \left( OTF \left( \lambda \right) \right)$$ (6)

$$MTF\left( \lambda \right) = Modulus \left( OTF \left( \lambda \right) \right)$$ (7)

PF	Pupil function
AT	Amplitude Transmittance (Stiles Crawford effect)
W	Wave aberration = f(Pupil diameter)
ASF	Amplitude Spread Function
OTF	Optical Transfer Function
RI	Retinal Image
I	Displayed image (CRT)

Default Parameters :
Pupil Diameter	3 mm,
Focus Wavelength	570 nm,
Zernike order	15,
Pupil shift X, Pupil Shift Y	0mm,
Wavelengths	400:10:750,
Displayed image size	512 pixels,
Pupil function size	256.

Images given for wavelegths: 400, 450, 500, 550, 600, 650, 700, 750.

Full images, Zoom-in (110 $\times$ 110 pixels), Surface plot (110 $\times$ 110).

Table 3: PSFs

			400nm
			450nm
			500nm
			550nm

Table 4: PSFs

			600nm
			650nm
			700nm
			750nm