A numerical method is presented to analyze the elastohydrodynamic lubrication of finite line contacts. The problem is analyzed using finite difference method and Newton-Raphson method. Pressure distribution and film shape are compared with the results for infinite line contacts. Near the center line of the roller, pressure distribution and film shape resemble those of infinite line contacts. However, near the start of the profiling, pressure and film shape are considerably different from those of infinite line contacts. The maximum pressure and minimum film thickness occur near the start of the profiling. The method used makes it possible to design an optimum axial profile of the roller to increase the life of rolling bearing.