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. (5), being sure to store your result back in $W$.

Hint: In octave, inv() and sqrtm(), respectively, invert and take the square root of a square matrix.