#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Tue Oct 29 18:42:09 2024

@author: sam
"""



'''Scalar Function with Explicit Surface'''

import numpy as np
import FinalProject as FP

def function(x,y,z):
    return y**2 + z

def surface(x,y): # paraboloid
    return x**2 + y**2

# where -1<x<1, -1<y<1
# expected result 8.56991



print(FP.integral(function,surface,[-1,1],[-1,1]))