Add diffraction command

[griffin28]

Procedure: Compute a diffraction pattern for a circular aperature.
Usage: diffraction radius [n-points]

[griffin28]

diffraction 5 100