  This paper describes the development of beam dynamic simulation code for
cyclotron. Starting from a description of beam dynamics in the cyclotron,
lattice functions were determined and the solutions for the 2nd-order nonlinear
Hamiltonian were revised. Based on the description of beam dynamics in the
cyclotron, simulation code was also developed for cyclotron design.
