To simulate the intense bunched beam transport, a computer program LEADS-3D has been developed. The particle trajectories are analysed with the Lie algebraic method. The third order approximation of the trajectory solutions is made with space charge forces off, and the second order approximation is made with space charge forces on. The particle distribution in the 3D ellipsoid is uniform or Gaussian. Most of the conventional beam optical elements are incorporated in the code. The optimization procedures are provided to fit the beam lines to satisfy the given optical conditions.
