Program to calculate the least-squares estimates of the spherical harmonic expansion coefficients of an equally angular-gridded scalar field