Using G as (bandwidth-1)^2 array and whitespace changes.

These results match the MATLAB code when using Cholesky. However, these
results do not match the backward-elimination / forward substitution
results.
Finished coding of compute_omega_prime().

NOTE: This routine still needs to be heavily tested. I suspect the linear
solve is still incorrect for some reason. These results do not exactly
match those of the MATLAB code.