Debugging

Modify your ``dft.m'' routine to repeat the simple-harmonic oscillator calculation with your minimal spectral representation, and verify that you reproduce your previous results to within 0.001 hartree!
Hints: The modification to ``dft.m'' should involve only changing the initialization of W to that of a random array of size length(gbl_active)$\times$Ns.
Hint: For debugging the above operators (most bugs will involve various ``nonconformant'' errors), make liberal use the size() function to print out the sizes of the relevant objects from within the functions that generate errors.