diff --git a/SIRC.m b/SIRC.m
index c5cc97b0e251e840af803051643cff8cb8675faa..ba0d994987a94f7dca5d232d967639af0e9b76b4 100644
--- a/SIRC.m
+++ b/SIRC.m
@@ -1,5 +1,5 @@
 %% System initalization
-% Number of population
+% Population
 N = 1000;
 % Initial number of infected individuals
 I0 = 10;
@@ -25,7 +25,10 @@ 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)
+plot(t, y);
+xlabel('Time(days)');
+ylabel('Number of individuals');
+legend('S', 'I', 'R', 'C');
 
 
 %% Add noise