  We discuss a recently proposed approach termed Q-BOR-FDTD method and develop
its main equations in a clearer and more rigorous way. We show that it is
unsuitable for the calculation of the eigenfunctions in the case of degenerate
states and propose an improvement to overcome such limitation.
