Recently, with the advance of the plasma diffusion theory of toroidal devices and the success of experimental results of the tokamak. octupole and spherator, quantitative comparison between the existing diffusion theory and the experimental results is necessary for an understanding of plasma diffusion in the magnetic trap.Here, the authors demonstrate the use of FORMAC – Formula Manipulation Compiler to carry out the analytical work and to prepare the corresponding portion of the FORTRAN program necessary to solve the equations numerically.