  We find approximate expressions x(k,n) and y(k,n) for the real and imaginary
parts of the kth zero z_k=x_k+i y_k of the Bessel polynomial y_n(x). To obtain
these closed-form formulas we use the fact that the points of well-defined
curves in the complex plane are limit points of the zeros of the normalized
Bessel polynomials. Thus, these zeros are first computed numerically through an
implementation of the electrostatic interpretation formulas and then, a fit to
the real and imaginary parts as functions of k and n is obtained. It is shown
that the resulting complex number x(k,n)+i y(k,n) is O(1/n^2)-convergent to z_k
for fixed k
