A Matlab based line contact EHL solver with artificial diffusion for highly loaded cases. The degree of the diffusion is controlled. Fully coupled, differential deflection based numerical scheme. Finite difference discretization, central difference approximation of the pressure term and 2-nd order backward difference approximation of the wedge term. Roelands pressure-viscosity relation. It is designed for solving highly loaded cases. Newtons numerical solver is used. Launch run_numerical.m.