Matlab codes to compute Quasi-Static Thermoelastic Deformation by Spherical Heat Sources. Details were shown in Furuya (2005, Geophysical Journal International).
I acknowledge Dr Yuki Nakashima who corrected a minor bug and reproduced the Figure 5 in Furuya (2005).