function yprime = epidemi(t,y,c,d)
    yprime = [-c*y(1)*y(2);...
             c*y(1)*y(2) - d*y(2);...
             d*y(2)];

end