In the current investigation, we present a numerical technique to solve Fredholm-Hammerstein integral equations of the second kind. The method utilizes the free shape parameter radial basis functions (RBFs) constructed on scattered points as a basis in the discrete Galerkin method to estimate the solution of integral equations. The accuracy and stability of the classical RBFs heavily depend on the selection of shape parameters. But on the other hand, the choice of suitable value for shape parameters is very difficult. Therefore, to get rid of this problem, the free shape parameter RBFs are used in the new method which establish an effective and stable method to estimate an unknown function. We utilize the composite Gauss-Legendre integration rule and employ it to estimate the integrals appeared in the method. Since the scheme does not need any background meshes, it can be identified as a meshless method. The error analysis of the method is provided. The convergence accuracy of the new technique is tested over several Hammerstein integral equations and obtained results confirm the theoretical error estimates.