Debugging

The main program in Appendix D solves for the total energy of the oxygen atom using your getg() routine. It initializes with zero density, Rho[]=0, and then iterates Rho$_{i+1}$[]=Rho[]+alpha*g[], as discussed in class, with alpha=0.25. This program should converge to a value of -74.473076 hartree in about 30 iterations.

Hint: If your return value from getg() is acting wierd, don't forget to #include ``physics.h'' at the top of your main programs (and include all of the function prototypes from physics.c into physics.h)!!! Several students in previous years have had this problem and this fixed it, so it may help you too!