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

@author: sam
"""



'''Scalar Function with Parametric Surface'''

import numpy as np
import FinalProject as FP

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

def surface(u,v): # cylinder in y-axis
    return np.array([np.cos(u),v,np.sin(u)])

# where 0<u<2π, 0<v<3
# expected result 18π ≈ 56.54867



print(FP.integral(function,surface,[0,2*np.pi],[0,3]))