Simpson's rule for a function of two variables

Good day respected geeks! Has any one of you ever written this thing (Simpson’s rule for a function of two variables) in Python?

I’m trying to find it but in vain. I’ve already coded it but it works really slow.

Many thanks in advance!