Commiting MKS 2nd order
4/11/14 end of week commit
Trying to debug so that I only have to find independent rows/columns for
the first 2 frequencies
Here I separate into an error checking version and a working version of
the 2nd order MKS code. It appears as though there are several different
sets of columns that rref or QR decomposition may fine as independent,
only the set with the lowest indices works for frequencies 1:el^3
(ie. [0,1,2,3], not [0,1,2,4] or [0,1,2,5].
The errorcheck version contains all of the different debugging methods I
was using for future reference.