[en] We are developing a series of computer codes (MAK0, MAK1, MAK2) which realistically model the performance of tokamak reactors in 0, 1, and 2 dimensions. Versions of MAK0 and MAK1 are operational and described in this report. They solve the time evolution of the density and energy conservation equations in D-T plasmas with alphas, neutrals, and two impurity species. MAK1 solves the plasma equilibrium equation on the Lagrangian flux-coordinate grid with transport coefficients computed as surface average quantities which include Pfirsh--Schulter, neoclassical, and trapped particle effects. Plasma recycle and impurity injection are self-consistent with charged particle and neutron outflow. An average ion coronal model determines impurity charges. Major radius control and plasma scrapeoff by a limiter are modeled. These codes are applicable to a wide variety of physics or engineering studies of tokamak reactors