I Started the code.  Currently, it has two contained functions, the initialization function, and the function that uses the redshifts to compute the critical density.  More changes to come.