Poincare Section MATLAB code (Pullen-Edmonds Hamiltonian)