#!/usr/bin/python3
import math
import numpy as np
import matplotlib.pyplot as plt
from scipy.special import p_roots #for pulling zeroes of Legendre orthogonal polynomials
from scipy.stats import qmc
import myquad


print(monty2d(f,region,a,b,M))