A new method for finding solutions to ordinary differential equation boundary value problems is introduced, in which Sobolev gradient steepest descent is used to determine eigenfunctions and eigenvalues simultaneously in ...