Das Kommando GramSchmidt berechnet die Orthonormalbasis zu einer gegebenen Liste reeller Vektoren. Es muß zunächst aus einem Mathematica-Package geladen werden.
(* Beispiel *) Remove["Global`*"] Needs["LinearAlgebra`Orthogonalization`"] mat={{1,2,3,4},{0,5,6,7},{4,5,6,7}} (* Drei 4-dim. Vektoren *) {on1,on2,on3}=GramSchmidt[mat]; mat1={on1,on2,on3} //MatrixForm (* Probe *) {on1 . on2, on1 . on3} {on1 . on1, on2 . on2, on3 . on3}