In this article, a novel and efficient approach based on Lucas polynomials is introduced for solving three-dimensional mixed Volterra–Fredholm integral equations for the two types (3D-MVFIEK2). This method transforms the 3D-MVFIEK2 into a system of linear algebraic equations. The error evaluation for the suggested scheme is discussed. This technique is implemented in four examples to illustrate the efficiency and fulfillment of the approach. Examples of numerical solutions to both linear and nonlinear integral equations were used. The Lucas polynomial method and other approaches were contrasted. A collection of tables and figures is used to present the numerical results. We observe that the exact solution differs from the numerical solution if the exact solution is an exponential or trigonometric function, while the numerical solution is the same when the exact solution is a polynomial. The Maple 18 program produced all of the results.