new function to calculate c(t) of Buxton

Implementation of eq [5] of Chapell's QUASAR paper