diff --git a/deriv.m b/deriv.m
new file mode 100644
index 0000000000000000000000000000000000000000..6e7dad951d917f210626576f3f95e6c85e749591
--- /dev/null
+++ b/deriv.m
@@ -0,0 +1,12 @@
+function dydt = deriv(y, nu, mu, epsilon, beta, gamma, Gamma, q)
+    S = y(1);
+    I = y(2);
+    R = y(3);
+    C = y(4);
+    dSdt = nu - (beta*I + epsilon*beta*C)*S - mu*S;
+    dIdt = (beta*I + epsilon*beta*C)*S - gamma*I - mu*I;
+    dRdt = gamma*(1-q)*I + Gamma*C - mu*R;
+    dCdt = gamma*q*I - Gamma*C - mu*C;
+    dydt = [dSdt; dIdt; dRdt; dCdt];
+end
+