Numerical Modelling of Stellar Differential Rotation and Dynamos