find_roche_lobe takes one step size

- Step size in each direction is calculated from the derivative of V
before equipotential evolution
- Function call changed in main.py
Fixed conflicts in lib.py