Removed input array "Lambda_star" from argument list of "calc_Lambda_star".
Removed extra input array "source_fn" from "calc_source_fn."

I forgot it was only there for the space and was abusing it terribly.
ALI looks exactly the way it does in OAB86 and OK87, but still doesn't give the correct answer...
Changed matrix elements of Lstar from individual arrays to a full matrix.

This is much less efficient storage-wise, but it makes the construction
of the full Lambda operator (or ALOs of variable bandwidth) much easier.
Changed Lstar to be the full Lambda operator so we can see if the matrix elements are correct.

Result: they are not correct. Hmm...