Natural orbitals implemented