In this paper we present an operational method for solving linear Volterra-Fredholm integral equations (VFIE). The method is con- structed based on three matrices with simple structures which lead to a simple and high accurate algorithm. We also present an error estimation and demonstrate accuracy of the method by numerical examples.