The numerical treatment of oscillatory integrals is a demanding problem in applied sciences, particularly for large-scale problems. The main concern of this work is on the approximation of oscillatory integrals having Bessel-type kernels with high frequency and large interpolation points. For this purpose, a modified meshless method with compactly supported radial basis functions is implemented in the Levin formulation. The method associates a sparse system matrix even for high frequency values and large data points, and approximates the integrals accurately. The method is efficient and stable than its counterpart methods. Error bounds are derived theoretically and verified with several numerical experiments.