function hogerled = epidemiODE(t, y, N, b, d, beta, u, v)

S = y(1);
I = y(2);
R = y(3);
hogerled = [ b*N - d*S - beta*I*S/N - v*S;
             beta*I*S/N - (u+d)*I;
             u*I - d*R + v*S ];

%              ELLER

% function hogerled = epidemiODE(t, y, N, b, d, beta, u, v)

% hogerled = [ b*N - d*y(1) - beta*y(2)*y(1)/N - v*y(1);
%              beta*y(2)*y(1)/N - (u+d)*y(2);
%              u*y(2) - d*y(3) + v*y(1) ];
