Added a get_mu(state) function which returns mu.

Fixed the lagrangian point solver. I think that works now.