hexagon logo

How to minimize the computational error of a particular measure which is solved through a subroutine?

The contact subroutine CNFSUB is giving slightly different result for the same set of parameters that is used in ADAMS VIEW contact module. What is the reason behind this difference?