diff --git a/SIRC.m b/SIRC.m
new file mode 100644
index 0000000000000000000000000000000000000000..c5cc97b0e251e840af803051643cff8cb8675faa
--- /dev/null
+++ b/SIRC.m
@@ -0,0 +1,33 @@
+%% System initalization
+% Number of population
+N = 1000;
+% Initial number of infected individuals
+I0 = 10;
+% Initial number of recovered individuals
+R0 = 0;
+% Initial number of carrier individuals
+C0 = 0;
+% Everyone else is susceptible
+S0 = N - I0 - R0 - C0;
+% Initial parameteres
+nu = 14;
+mu = 0.015;
+epsilon = 1.2;
+beta = 0.0002;
+gamma = 0.05;
+Gamma = 0.03;
+q = 0.8;
+
+time = 200;
+%t = linspace(0, time, time);
+
+y0 = [S0, I0, R0, C0];
+[t, y] = ode45(@(t, y) deriv(y, nu, mu, epsilon, beta, gamma, Gamma, q), [0 time], y0);
+
+figure(1)
+plot(t, y)
+
+
+%% Add noise
+
+