  We discuss the automatic solution of the multichannel Schr\"odinger equation.
The proposed approach is based on the use of a CP method for which the step
size is not restricted by the oscillations in the solution. Moreover, this CP
method turns out to form a natural scheme for the integration of the Riccati
differential equation which arises when introducing the (inverse) logarithmic
derivative. A new Pr\"ufer type mechanism which derives all the required
information from the propagation of the inverse of the log-derivative, is
introduced. It improves and refines the eigenvalue shooting process and implies
that the user may specify the required eigenvalue by its index.
