The linear stability of equilibria of charged particles moving near a compact object with a dipole magnetic field and a pseudo-Newtonian potential is analyzed detailedly. An optimal fourth-order force gradient symplectic method, as a global symplectic integrator that can simultaneously solve both the equations of motion and the variational equations, is used to calculate fast Lyapunov indicators. In this way, dynamical structures are described, and parameter domains for causing chaos are found. PACS numbers: 05.45.-a, 04.70.-s, 05.10.-a, 95.10.Fh