new function to calculate c(t)

equation [9] of QUASAR paper