The Lie algebraic method is applied to the analysis of the nonlinear transport of an intense pulsed beam in cylindrically symmetrical electrostatic lenses, and particle orbits in a six-dimensional phase space (x, px, y, py, τ, pτ)are obtained in the second order approximation. They can also be acquired in the third or higher order approximation if needed. In the analysis, we divide the electrostatic lenses into several segments. Each segment is considered as a uniform accelerating field, and each dividing point is treated as a thin lens. The particle distribution in a three-dimensional ellipsoid is of Gaussian type.