An effective optimization strategy has been developed to construct highly
accurate bound state wave functions in various three-body systems. Our
procedure appears to be very effective for computations of weakly bound states
and various excited states, including rotationally excited states, i.e. states
with L≥1. The efficiency of our procedure is illustrated by computations
of the excited P∗(L=1)−states in the ddμ,dtμ and ttμ muonic
molecular ions, P(L=1)−states in the non-symmetric pdμ,ptμ and
dtμ ions and 21P(L=1)− and 23P(L=1)−states in He atom(s)