(1f)  d C g    d t   = ρ ϵ H g   − ψ C g