In python please

Write a program that demonstrates the Fourier Series method as defined on the MAE 3013 equation sheet. You must write and call at least the following 3 functions: def FourierCoeffs (func, L, nterms): #calculates and returns the Fourier a's and b's as numpy arrays of length nterms. func is the name of the function. L is the half-period. This function must call scipy.integrate.quad def Plot Fourier (a, b, L, xmin, xmax, npoints=5000): #which plots the Fourier Series def main(): def Sharkfin(x): if x