Copied and pasted in C++ code for orbital calculations
Need to convert this to Python!