In the LCAO model, the molecular orbitals are approximated as a weighted sum of atomic orbitals. The coefficients of this sum are determined by solving the Schrödinger equation for the system. This involves constructing and diagonalizing the Hamiltonian matrix, which represents the total energy of the system.