Solution to Self-Evaluation Exercise 2: MATLAB code for the SIPDB model