The first-order optical properties are computed by numerical solution of the paraxial ray equation, and the computed paraxial rays are displayed on the screen. The primary aberrations are computed by evaluating appropriate aberration integrals involving the paraxial rays.  The computed geometric aberrations are second and third-order aberrations, if hexapole lenses are present, and third-order aberrations, otherwise.  The program can handle both Gaussian round-beam and shaped-beam systems.  The output is the computed optical properties, including a table of the first-order optical properties and a table of the aberration coefficients.  MULTIPOLE also has an integrated aberration spot generator to visually show the aberration effects using the computed aberration coefficients and given initial image conditions, such as beam half angle, beam size, deflection parameter at either object or image plane.