Используя y(t,m,i) = q(t,m,i) + eps/MI
и ln(a/b)=ln(a)-ln(b)
, нужно беспокоиться о двух терминах:
f(y(t,m,i)) = y(t,m,i) ln y(t,m,i)
и
g(y(t,m,i),y(t-1,m,i)) = y(t,m,i) ln y(t-1,m,i)
Первый прост.Второй не так много, см. ссылка для некоторых заметок.