Roche lobe solver works, but sensitive to initial conditions and step sizes
Merge branch 'master' of https://github.com/dkong1796/ay190-final