Finding optimal convergence control parameter in the homotopy analysis method to solve integral equations based on the stochastic arithmetic