Hi,
I understand that you are facing an error while solving the symbolic differential equations using dsolve. This is due to dsolve not being able to find a symbolic solution to the equations.
This can be solved by using a numeric solver like ode45 as:
% Parameters for the SIR model
beta = 0.3; gamma = 0.1; mu = 0.05; nu = 0.02; eps = 0.1;
% Define the system of differential equations
sir_ode = @(t, y, eps) [
mu - beta * y(1) * y(2) - mu * y(1) + eps * (mu - beta * y(1) * y(2) - mu * y(1));
beta * y(1) * y(2) - gamma * y(2) - nu * y(2) + eps * (beta * y(1) * y(2) - gamma * y(2) - nu * y(2));
gamma * y(2) - mu * y(3) + eps * (gamma * y(2) - mu * y(3))
];
% Initial conditions
y0 = [0.9; 0.1; 0]; tspan = [0 100];
% Solve the system using ode45
[t, y] = ode45(@(t, y) sir_ode(t, y, eps), tspan, y0);
You can refer to the below documentation for more information on dsolve function:
Hope it helps!
