We combined a new theory of the LAMBDA-effect with mixing-length descr
iption of the solar convection zone to develop a model of the Sun's di
fferential rotation. Only two adjustable parameters are involved: the
coefficient alpha(MLT) of the known mixing-length relation, l = alpha(
MLT)H(p), and the mean angular velocity, OMEGA. Both parameters are ne
arly fixed for the Sun. Nevertheless, the alphaOMEGA-model is successf
ul in reproducing both the helioseismology inversions of the internal
solar rotation and the observed angular velocity distribution at the s
urface. This encouraged us to proceed with simulating the rotation law
s for the rapidly rotating ''young Sun'' and slowly rotating ''old'' o
ne. The angular velocity distribution for the former case is close to
the present-day state. Our prediction for the old Sun is a negligible
latitudinal rotation profile and a pronounced super-rotation in radius
.