The method can be used to calculate the equilibrium orbit of particles in those circular accelerators with very large magnetic flutter,say,the fixed field alternating gradient(FFAG) accelerator.