The transport equations of the bumpy torus EBT-I and EBT-P are solved numerically. The transport equations are derived from the drift kinetic model. In the case of the EBT-I, the boundary values from experiments are used, but in the case of the planned EBT-P, since there are no experimental values available, the boundary values which give the stable states are obtained simultaneously with the solutions of transport equations. The solutions give the profiles of the temperatures and densities of charged particles.