Taking Newton-Raphson approach to solve elastohydrudynamic lubrication (EHL) problem, a complete numerical solution with typical EHL characteristics of oil film thickness and pressure profile is obtained. Numerical calculating expressions of traction regarding Newton model. Maxwell model, J-T model and B - W model are given herein, with the first two pointedly numerically calculated. Complete numerical case solution to isothermal EHL with line contact is then worked act, perfectly coinciding with expectation.