The iterative computation of singular points in parametrized nonlinear BVPs by so-called extended systems requires good starting values for the singular point itself and the corresponding eigenfunction. Using path-following techniques such starting values for the singular points should be generated automatically. However, pathfollowing does not provide approximations for the eigenfunctions. We ...