The ability to simulate mixed lubrication problems has greatly improved, especially in concentrated lubricated contacts. A mixed lubrication simulation method was developed by utilizing the semi-system approach which has been proven to be highly useful for improving stability and robustness of mixed lubrication simulations. Then different variants of the model were developed by varying the discretization schemes used to treat the Couette flow terms in the Reynolds equation, varying the evaluation of density derivatives and varying the contribution of terms in the coefficient matrix. The resulting pressure distribution, film thickness distribution, lambda ratio, contact ratio, and the computation time were compared and found to be strongly influenced by the choice of solution scheme. This indicates that the output from mixed lubrication solvers can be readily used for qualitative and parametric studies, but care should be taken when making quantitative predictions.