Вычислить коэффициенты Фурье с помощью Python?

Я пытаюсь вычислить следующие коэффициенты Фурье

где V_{pot} — предыдущая функция определения этой формы.

Я действительно не знаю, какой численный метод я могу использовать, однако я начал с правила Симпсона библиотеки scipy.

import numpy as np
from scipy.integrate import simps

Nf = 200 
IVp = np.zeros(2*Nf)
snn = np.zeros(NP)
def f(k):
    for i in range(0,NP):
        sn = (i-1)*H
        snn[i] = sn
        return (1/SF) * np.cos(np.pi*k*sn/SF) * Vpot(sn)

 for k in range(0,2*Nf):
      Func = f(k)
      y1 = np.array(Func,dtype=float)
      I = simps(y1,snn)

У меня была эта ошибка:

IndexError: tuple index out of range