Initialize W

It is helpful to at least start with orthonormal wave functions. Thus, immediately after the call to fdtest in your ``sch.m'' file, orthonormalize $W$ according to equation Eq. (6), being sure to store your result back in $W$.

Hint: In octave, inv() and sqrtm(), respectively, return the matrix inverse and square root of a matrix, as opposed to simply inverting or taking the square root of each matrix element separately.