solve u''(r) + u'(r)/r == -lambda u(r)