modify r(t) and m(t)

new r(t) and m(t) is based on eq [5] of Chappell’s QUASAR paper
change dispersion function

implementation of eq [6] of Chappell’s QUASAR paper