In this paper, we discuss the existence and uniqueness of the solution of the second kind nonlinear Volterra-Fredholm integral equations (NV-FIEs) which appear in mathematical modeling of many phenomena, using Picard’s method. In addition, we use Banach fixed point theorem to show the solvability of the first kind NV-FIEs. Moreover, we utilize the homotopy analysis method (HAM) to approximate the solution and the convergence of the method is investigated. Finally, some examples are presented and the numerical results are discussed to show the validity of the theoretical results.