fixed classify
plot orbit name
rudimentary tests for core dynamics
add hernquist potential
use c version
subclass tests
wrong public defs
was using numpy sqrt