In order to solve the double integral approximation problem,a high order numerical integral formula was applied to extend it to the double integrals.Then,analyze the remainder and error of this formula.Finally,the high validities of this formula are proved by several typical examples.