NOT WORKING YET:

need to sort out the a_dependence and get it in terms of conformal time
need to work on the case when you get logs and how you pass that through etc
implemented the w>1/3 condition to swap the limits
need to change the a values for the pow spectrum to math.pow(time,a_pow)
should theoretically work after that.....