# -*- coding: utf-8 -*-
"""
Created on Fri Sep  9 09:59:13 2022

@author: Chris
Lab7_Q6
"""

import myquad as myq
    
def region(x,y,z):
    return (x/2+y/2+z/3 < 1)&(x>0)&(y>0)&(z>0)

def f(x,y,z):
    f = z-x*y**2
    return f

a,b = -5,5
M = 20
I = myq.monty3d(f,region,a,b,M)

print(I)