Initial two-body code

Derivatives are wrong.
Two body derivatives work now.