Added a new file for the functions that generalize the coefficient calculation to any spherical harmonic and do stuff with the full Q matrix. It's not yet reproducing the results for l,m=0,0, so it's still in progress.