Solution to 7.9 #32 with Maple.

> with(LinearAlgebra):

> GramSchmidt([<1,2,1,0>,<1,1,2,0>,<0,1,1,-2>,<1,0,3,1>]);

Notice it did not normalize the vecotrs. But, we can add this as an option.

> GramSchmidt([<1,2,1,0>,<1,1,2,0>,<0,1,1,-2>,<1,0,3,1>],normalized);

>