A=[-1,2,1;
   2,5,-3;
   1,-3,0];
[B,P]=GaussReduction(A);
disp('symmetric matrix of quadratic form:')
A
disp('orthogonal basis for quadratic form (columns of matrix):')
P
disp('diagonal matrix of reduced quadratic form:')
B
