In this paper, utilizing the technique of Petryshyn’s fixed point theorem in Banach algebra, we analyze the existence of solution for functional integral equations, which includes as special cases many functional integral equations that arise in various branches of non-linear analysis and its application. Finally, we introduce the numerical method formed by modified homotopy perturbation approach to resolving the problem with acceptable accuracy.