Introduction. - In Mathematical Analysis, integrals are sometimes difficult to solve by integration methods, for this reason it is necessary to think of numerical methods to solve them, if these are definite integrals, even more so when it comes to double integrals. Objective: The ability to solve these integrals depends a lot on the solution knowledge and experience, for this reason the present study presents an alternative solution by numerical methods. Methodology: A double integral conceptually calculates the volume bounded by a surface over a region, the region may be rectangular or non-rectangular. The present study solves the type of indifferently, for which partitions are made along x and y of the region, thus showing a mesh within points (x,y) of the region, which are evaluated in the function f(x,y), which represents the surface, Results: with these values the integral is solved horizontally using the Simpson method and with the result of these it is solved vertically with the same method, obtaining the result of the double integral with excellent precision, Conclusions: a method of calculating double integrals  is then proposed with numerical calculation for rectangular or non-rectangular regions.Introducción: En análisis matemático en ocasiones las integrales resultan complicadas resolverlas por métodos de integración, por esta razón es necesario pensar en métodos numéricos para resolverlas, siempre que estas sean integrales definidas, más aún cuando se trata de integrales dobles. El objetivo: La habilidad para resolver estas integrales depende mucho del conocimiento de solución y de la experiencia, por esta razón el presente estudio presenta una alternativa de solución por métodos numéricos. Metodología: Una integral doble conceptualmente calcula el volumen limitado por una superficie sobre una región, la región puede ser rectangular o no rectangular. El presente estudio resuelve indiferentemente el tipo de región, para lo cual se realiza particiones a lo largo de  e  de la región, generando de esta manera una malla de puntos  dentro de la región, los cuales se evalúan en la función , que representa la superficie, Resultados: con estos valores se resuelve en forma horizontal la integral mediante el método de Simpson. Con el resultado de estos se resuelve en sentido vertical con el mismo método, obteniéndose el resultado de la integral doble con excelente precisión. Conclusiones: se propone entonces un método de cálculo de integrales dobles con cálculo numérico para regiones de rectangulares o no rectangulares.