2D_SH_matlab (Forked from Luca's code)