%funkcija DY=crvi(t,Y,param) izracuna desne strani ODE
%param=[N,eta,gama,mu*N,beta0*N]
%Y = [I,J,C]
function dY=crvi(t,Y,param)
N=param(1);
eta = param(2);
gama = param(3);
mu = param(4)/N;

beta=param(5)/N*(1-Y(1)/N)^eta;
D=N-Y(2)-Y(3);
dY=[(beta*D-gama)*Y(1);beta*D*Y(1);mu*D*Y(2)];
