Paraxial expansion of a static magnetic field in a ring accelerator

In the paraxial approximation, trajectories of beam particles in a ring accelerator are computed expanded in powers of lateral displacements and slopes from a closed reference orbit. To do this, one needs first the expanded expressions of the magnetic field and potentials producing the particle motion. This is derived here in a most general form