Delaminating quadrature method for multi-dimensional highly oscillatory integrals

A delaminating quadrature method for multi-dimensional highly oscillatory integrals is put forward. In comparison with the Levin method, the new method has the benefits of being numerically stable and having a high computational speed. Numerical experiments testified the advantages of this method.