Reorganised methods - now includes new 'coefficients' object
to deal with the shapelet coefficients

Currently writing methods to calculate total flux and centroid from the coefficients

Total flux works, but centroids are off by ~3-4 (?)