A numerical method for evaluating the Hankel transform of order of function is given. The method is based upon B-polynomial multi-wavelets forms an orthonormal bases for we expand the part of the integrand in its wavelet series reducing the Hankel transform integral reduces a series of Bessel function multiplied by the wavelet coefficients of the input function. The proposed method has been illustrated by examples.