Add noise to the system

SIRC.m

@@ -32,5 +32,20 @@ legend('S', 'I', 'R', 'C');
 
 
 %% Add noise
+noise = wgn(size(t,1), 4, 0, 42);
+y_noise = y + noise;
+
+figure(2)
+plot(t, y_noise);
+xlabel('Time(days)');
+ylabel('Number of individuals');
+legend('S', 'I', 'R', 'C');
 
 
+[t_noisy, y_noisy] = ode45(@(t, y) noisy_deriv(y, nu, mu, epsilon, beta, gamma, Gamma, q), [0 time], y0);
+
+figure(3)
+plot(t_noisy, y_noisy);
+xlabel('Time(days)');
+ylabel('Number of individuals');
+legend('S', 'I', 'R', 'C');