printf("GNU Octave, version 3.4.3\nOctave was configured for \"x86_64-redhat-linux-gnu\".\nFedora 16. Running Octave in QtOctave.\n");
X=[0,1,0,0;1,0,0,0;0,0,1E-30,0;0,0,0,5]
printf("X(3,3):");
X(3,3)
printf("rref(X):");
rref(X)
printf("Is not equal to:\n");
Correct=[1,0,0,0;0,1,0,0;0,0,1,0;0,0,0,1;]
printf("Which it should be. For all values of X(3,3) not equal to 0, NaN, or +-Inf this \nmatrix should row reduce to [Correct].\n");
