added fit_disk.py to model 2-d stellar distribution with sums of gaussians.