(1d)  d I g    d t   = η δ E g   − μ I g