Advanced topics for further programming development

To ensure good numerical conditioning of the matrix He and to prevent overflow or underflow in evaluating the inverse of He, the introduction of a local non-dimensional coordinates system is usually suggested [6, 9].