I have attached a copy of my program. I have a numerical solution generator for my function as well as an analytic solution that I solved out by hand. For some reason at particular initial conditions and values of my constants the results diverge heavily. Any suggestions on methods to clean up the comparison for wider ranges of values? I want the two plots to match as close as possible for as many ranges of my variables as possible.
Attachment: Cosmology_Plots_New.nb, URL: ,