rho t f(t) delta(r/(a(t))) + rho t