Mixing of molecular orbitals in polyatomic molecules
The standard way is to solve for the eigenvectors of your Fock matrix. Since you have a single electron, that's your Hamiltonian. Any number of programming languages have algorithms for this implemented. Each of your coefficients would correspond to an entry in the eigenvector.