We give two distinct approaches to finding bounds, as functions of the order i/, for the smallest real or purely imaginary zero of Bessel and some related functions. One approach is based on an old method due to Euler, Rayleigh, and others for evaluating the real zeros of the Bessel function Jv(x) when v > -1. Here, among other things, we extend this method to get bounds for the two purely imaginary zeros which arise in the case -2 < v < -1. If we use the notation j u i for the smallest positive zero, which approaches 0 as v -> -1 + , we can think of j^ as continued to-2 < v < -1, where it has negative values. We find an infinite sequence of successively improving upper and lower bounds for j^ in this interval. Some of the weakest, but simplest, lower bounds in this sequence are given by 4(z/ + 1) and 2 5 / 3 (i/ + l)[(i/ + 2)(u + 3)] 1 / 3 while a simple upper bound is 4(i/ + l)(v + 2) 1 / 2 . The second method is based on the representation of Bessel functions as limits of Lommel polynomials. In this case, the bounds for the zeros are roots of polynomials whose coefficients are functions of v. The earliest bounds found by this method already are quite sharp. Some are known in the literature though they are usually found by ad hoc methods. The same ideas are applied to get bounds for purely imaginary zeros of other functions such as J' v {x), J"(x), and OLJ V (X) + xJl(x).