Added hamiltonian jacobi work through in mathematica to automagically work out ns and plot the potential

couldnt get the third order eqn to solve!