MATLAB code for SSIIR model