Numerical solutions are obtained for the thermoelastohydrodynamic lubrication of a tilting pad based on the three-dimensional flow of lubricant. For elastic deformation, the pad is idealized as a uniform plate with free boundaries so that a displacement coefficient matrix can be established by the conventional finite difference method. For pressure distribution, a fast Newton-Raphson's approach is used in solving the generalized Reynolds equation, film thickness equation, inlet pressure buildup equation, and the force balance and moment balance equations of the pad, simultaneously. For temperature distribution, the finite difference method is carried out with the help of a sweeping scheme.