Eigenvalue problems arise in many areas of physics, from solving a classical electromagnetic problem to calculating the quantum bound states of the hydrogen atom. In textbooks, eigenvalue problems are defined for linear problems, particularly linear differential equations such as time-independent Schrödinger equations. Eigenfunctions of such problems exhibit several standard features independent of the form of the underlying equations. As discussed in Bender et al (2014 J. Phys. A: Math. Theor. 47 235204), separatrices of nonlinear differential equations share some of these features. In this sense, they can be considered eigenfunctions of nonlinear differential equations, and the quantized initial conditions that give rise to the separatrices can be interpreted as eigenvalues. We introduce a first-order nonlinear eigenvalue problem involving a general class of functions and obtain the large-eigenvalue limit by reducing it to a random walk problem on a half-line. The introduced general class of functions covers many special functions such as the Bessel and Airy functions, which are themselves solutions of second-order differential equations. For instance, in a special case involving the Bessel functions of the first kind, i.e. for y′(x) = J
ν
(xy), we show that the eigenvalues asymptotically grow as 241/42
n
1/4. We also introduce and discuss nonlinear eigenvalue problems involving the reciprocal gamma and the Riemann zeta functions, which are not solutions to simple differential equations. With the reciprocal gamma function, i.e. for y′(x) = 1/Γ(−xy), we show that the nth eigenvalue grows factorially fast as
(
1
−
2
n
)
/
Γ
(
r
2
n
−
1
)
, where r
k
is the kth root of the digamma function.