For my final project for my cosmology class, I numerically investigated single-field scalar inflation. The potential I used, quadratic inflation, has been ruled out by Planck but it was still nice to do the calculations of the spectral index and tensor-to-scalar ratio and see I get the same results. I then used Monte Carlo to show how non-Gaussianities arise from inflation and why they’re so hard to find. I’d like to extend this project to investigate eternal inflation, but that’s for the future. You can find my paper here (the code is TeXed in).