This just provides a test that the NFW potential is giving the correct orbit 
for an elliptical orbit as calculated by Jo Bovy's galpy package. If 
galpy is not installed on your system you can install it by using:
pp install galpy --user