Mathematica Notebook that generated the initial data for the fortran programs. Each fortran program is for different values of m and k and solves the mode equation for the Bunch-Davies state for that value of m and k. For this figure m = 10. The starting values for the modes which come from the Mathematica notebook were copied and pasted into the Fortran programs. Note that the Fortran programs all use the Numerical Recipies routine odebs and its subroutines to numerically solve the mode equation. This can easily be replaced by a call to any differential equation solver.

Some of the work referenced on this page was supported in part by the National Science Foundation
under grant numbers PHY-0856050 and PHY-1308325. Any opinions, findings and conclusions or recomendations expressed in this material are those of the author(s) and do not necessarily reflect the views
of the National Science Foundation (NSF).