Added a check routine (that calls gauss_parse.py)
Fixed SCF (works for everything but CO)
Started on a new density guess
Added test cases with corresponding calls to test_out.m