an algorithm for numerical double integration over non-rectangular regions.

[okkvlt]

is there one that is stable and accurate?

[HallsofIvy]

What are you referring to? A numerical algorithm? No, there certainly will not be a numerical algorithm that will work over every smooth region, much less every region.