Galaxy modelling code revision

Revision of the galaxy source modelling code, including the following changes:
+ Implementation of the spatial and velocity distribution models as `astropy.modeling.Fittable1DModel`
+ Use of `astropy.units.Quantity` throughout the code
+ Added tests
+ Added minimal docs
This comment is out of place here ... the `fsolve` call is further up.