Matlab code for Hopenhayn-Nicolini's optimal unemployment insurance model