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

@author: sam
"""



'''Vector Function with Explicit Surface'''

import numpy as np
import FinalProject as FP

def function(x,y,z):
    return np.array([x,0,-z**2])

def surface(x,y): # plane independent of x
    return 2-y

# where 0<x<2, 0<y<2
# expected result 16/3 ≈ 5.33333



print(FP.integral(function,surface,[0,2],[0,2]))