In order to model the flow of nematic crystals, the theoretical framework according to Ericksen and Leslie is applied. The essentials of the theory are compiled and then specialized to Couette flow. The profiles for linear velocity and orientation angle will be computed and, in particular, we shall also study the rise in temperature due to viscous dissipation, which is frequently ignored by mechanicians. Analytical and numerical solutions for the fields are derived for different boundary conditions and will subsequently be discussed.