The article is devoted to common methods for creating formulas of the evaluating 
of two dimensions highly oscillating integrals with using spline-interlination. 
̉he new formulas are optimal by order of exactness. The functions with high exactness 
of approximation were fined.