Add titles

SIRC.m

 %% 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