plot values not deviation
add integrate orbit function to potential
plots
mean not median