Abstract: The main objective of this study is to derive two new rules numerically for finding the values of the continuous double integrals in region of the integration by using mid-point rule and Simpson’s rule on the x and y dimensions and how to find the correction terms (error formula) and we use Romberg acceleration to improve the results when the number of subintervals of on the two dimensions are not equal.