In order to research multi-component ion beam transport process and improve transport efficiency, a special simulating program for ECR beam is becoming more and more necessary. We have developed a program written by Visual Basic to be dedicated to numerical simulation of the highly charged ion beam and to optimization of beam dynamics in transport line. In the program the exchange of electrons between highly charged ions and low charged ions or neutral atoms (residual gas in transport line) is taken into account, adopting classical molecular over-barrier model and Monte Carlo method, so the code can easily give the change of charge state distribution along the transmission line. The main advantage of the code is the ability to simultaneously simulate a large quantity of ions with different masses and charge states, and particularly, to simulate the loss of highly charged ions and the increase of low charged ions due to electron exchange in the whole transport process. Some simulations have been done to study the transmission line of LECR3^[1] which is an ECR ion source for highly charged ion beam at IMP. Compared with experimental results, the simulations are considered to be successful.