How is the orthogonalization symmetric in the "bst_henv.m" function?

I see, thank you very much for the clarification!