Added simple safeguards against erratic points

- If the points start too close to a star, then it has strange behavior.
- Simple safeguard in roche lobe solver to detect this as well.
- Needs lots of refinement - ideally values shouldn't be hard-coded, but
  should depend on the masses/mass ratio and h.