The aim of this article is focused on developing an efficient algorithm for simulating Cellular Neural Network arrays (CNNs) using numerical integration techniques. The role of the simulator is that it is capable of performing raster simulation for any kind as well as any size of input image. It is a powerful tool for researchers to investigate the potential applications of CNN. This article proposes an efficient pseudo code for exploiting the latency properties of Cellular Neural Networks along with well known numerical integration algorithms. Simulation results and comparison have also been presented to show the efficiency of the numerical integration algorithms. It is observed that the Runge-Kutta (RK) sixth order algorithm outperforms well in comparison with the Explicit Euler, RK-Gill and RK-fifth order algorithms.