%funkcija DY=crvi(t,Y,param) izracuna desne strani ODE
%param=[N,eta,gama,mu*N,beta0*N]
%Y = [D,Z,C]
function DY = crvi(t,Y,param)
N = param(1);
eta = param(2);
gama = param(3);
mu = param(4)/N;
I = (param(1)-sum(Y));
beta = param(5)/N*(1-I/N)^eta;
DY = [Y(1)*(-beta*I-mu*Y(1)*(I+Y(2)));gama*I;mu*Y(1)*(I+Y(2))];
