Simulations written in python, using the library QuTiP. The equation for time evolution of a photon under parametric down conversion was adapted from arXiv:1608.05005. Code for plotting the average photon number and variances were adapted from qutip-lectures.
matplotlib, scipy, numpy, qutip