Elastohydrodynamic point contacts part I: Formulation and numerical solution