The equilibrium equation for a rotating plasma is constructed supposing the thermal Mach number is much less than unity. The canonical profile of angular rotation velocity is defined as the profile which minimizes the total plasma energy while conserving toroidal current and obeying the equilibrium condition. The transport model based on this canonical profile, with stiffness calibrated by JET data, reasonably describes the velocity of the forced toroidal rotation. The RMS deviations of the calculated rotation profiles from the experimental ones do not exceed 10-15%. The developed model is also applied to the modeling of MAST rotation.