Abstract: A program for calculating intense pulsed beams tran sp ort in axialsymmetrical electrostatic fields is designed. The program can calc ulate beam optical systems consisting of the following elements: drift spaces, t hreecylinder einzel lenses, threediaphragm einzel lenses, AC electrostatic a ccelerating tubes, as well as twocylinder lenses. In the calculations, the eff ective fields of each lens are divided into several small intervals, and each in terval is treated as a uniform accelerating field, each dividing point is consid ered as a thin lens, and iterations are used to get selfconsistent solutions. Iteration procedures are incorporated in the codes to get selfconsistent solut ions.