The stationary Reynolds equation is solved over a rectangular region. The problem is linearized by Picard linearization. The ADI method is used to solve the resulting set of linear equations. A set of parameters is introduced to speed up convergence as well for the Picard linearization as for the ADI method. A comparison is made with Booy-Coleman's method. Results are given for bearing numbers 10 to 1000.