% % SIR_rhs_complex % function dy = SIR_rhs_complex(t,y,params); global N beta = params(1); gamma = params(2); mu = params(3); S = y(1); I = y(2); % % dy = [mu*(N-S) - beta*I*S; beta*I*S - (gamma+mu)*I];