At present, some scientific calculation software can get approximate solution of complex integration through numerical integration. However, scientific calculation software can only solve triple integration at most. To solve the problem of limited multiplicity in integral calculation, an adaptive recursive fast calculation method for arbitrary multiple integrals was proposed based on the analytic calculation idea of multiple integrals and recursive algorithm in this paper. The algorithm and recursion process of the method were introduced in principle, and the source code of the method was given in Matlab. The method was compared with the existing methods through examples. The results show that the effectiveness of the method, which can reduce the calculation time greatly. It can be solved only by providing integral interval and integrand function, eliminating the programming steps of matching program code with integral multiplicity. The method is simple in operation, high in calculation accuracy and speed, which can meet the actual engineering needs.