Updated the function so that err_matrix in spherical harmonics has a 4pi
less (though this seems to change the error in the wrong direction...)
Slightly update to the code (fixing the 4pi integration in the spherical
harmonics error)