In this paper an algorithm is developed where Reynolds' equation, equilibrium equations and non-negativity of pressure are formulated as a system of equations, which are not differentiable in the usual sense. This system is then solved using Pang's Newton method for B-differentiable equations.