Using Grad’s theory of generalized differential equations, the temporal evolution from low to high &bgr; due to ’’adiabatic’’ and nonadiabatic (i.e., neutral beam injection) heating of a cylindrical tokamak plasma with circular cross section and peaked current profiles is calculated analytically. The influence of shaping the initial safety factor profile and the beam deposition profile and the effect of minor radius compression on the equilibrium is analyzed.