Abstract: Introducing chaotic factor as a local operator of cuckoo search algorithm,a method was presented to calculate two dimensional numerical integration based on chaotic cuckoo search (CCS) optimization algorithm.CCS algorithm selected some initial partition points randomly on the both direction range in rectangular integral domain and then to optimize these points.Points optimized would be as the segmentation nodes to solve two dimensional numerical integration.The experimental results show that our algorithm is easy to implement and has fast convergence speed and high integral value accuracy.Therefore,CCS is an effective method for solving the two dimensional numerical integration.