Оценка параметров для эпидемических моделей

Question

У меня есть модель эпидемии:

function ypsirid = ypsirid(t,y)
global q m b r
ypsirid(1) = q - m*y(1) - b*y(1)*y(2);
ypsirid(2) = (b*y(1)*y(2)) - (m+r)*y(2);
ypsirid(3) = (r*y(2)) - m*y(3);
ypsirid = ypsirid(:)

Может ли кто-нибудь любезно предложить самый простой из возможных методов для оценки параметров с использованием matlab. Спасибо!