MATLAB code for Figure 4 from Indentation into a plastic fluid layer

MATLAB code to reproduce figure 4. Code solves system of equations (3.21)-(3.25) in section 3(d) for a Herschel Bulkley substrate with n=1/3.