Abstract: The method presented is based on a similarity transformation of the coefficients matrix allowing the transformed variables to be relaxed independently in a recursive fashion.