We propose an undergraduate numerical project for simulating the results of the second-order correlation function as obtained by an intensity interference experiment for two kinds of light, namely bunched light with Gaussian or Lorentzian power density spectrum and antibunched light obtained from single-photon sources. While the algorithm for bunched light was based on a previous publication, the one for antibunched light was developed by master students attending a Quantum Optics course at our University. (Contains 7 figures.)