Das Kommando GramSchmidt berechnet die Orthonormalbasis zu einer gegebenen Liste reeller Vektoren. Es muß zunächst aus einem der Mathematica-Packages geladen werden.
(* Beispiel *) Remove["Global`*"] Needs["LinearAlgebra`Orthogonalization`"] v1:={a1,b1,c1}; v2:={a2,b2,c2}; {on1,on2}=GramSchmidt[{v1,v2}]