add in new package for oblate star calculation, see test_Zeipel.py