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 = 1. 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.

Fortran program for m = 1, k = 1

Fortran program for m = 1, k = 10

Fortran program for m = 1, k = 100

Fortran program for m = 10, k = 1

Fortran program for m = 10, k = 10

Fortran program for m = 10, k = 100

Data file for the top plot in Fig. 6

Data file for the bottom plots in Fig. 6

Data file for the top plots in Fig. 7

Data file for the middle plots in Fig. 7

Data file for the bottom plots in Fig. 7
For the plots in Figure 7, I multiplied the contributions to the energy density by the 3rd power of the scale factor. The first column is the dimensionless time u and the scale factor is Cosh[u]/H.

Data file for the top plots in Fig. 8

Data file for the middle plots in Fig. 8

Data file for the bottom plots in Fig. 8
For the plots in Figure 8, I multiplied the contributions to the energy density by the 3rd power of the scale factor. The first column is the dimensionless time u and the scale factor is Cosh[u]/H.

Data file for the top plots in Fig. 10

Data file for the middle plots in Fig. 10

Data file for the bottom plots in Fig. 10
For the left hand plots in Figure 10, I multiplied the contributions to the energy density by the 3rd power of the scale factor. For the right hand plots, I multiplied the contributions to the energy density by the 4th power of the scale factor. The first column in the data files is the dimensionless time u and the scale factor is Cosh[u]/H.

Data file for the top plots in Fig. 11

Data file for the middle plots in Fig. 11

Data file for the bottom plots in Fig. 11
For the left hand plots in Figure 11, I multiplied the contributions to the energy density by the 3rd power of the scale factor. For the right hand plots, I multiplied the contributions to the energy density by the 4th power of the scale factor. The first column in the data files is the dimensionless time u and the scale factor is Cosh[u]/H.