Нахождение н-ой частичной суммы ряда Фурье функции с использованием ядер Дирихле

import matplotlib.pyplot as plt
from math import *
 
print('Нахождение н-ой частичной суммы ряда Фурье функции с использованием ядер Дирихле')
print("Введите L")
L = int(input())
print("Введите N")
N = int(input())
ArgF = []
SND = []
ResF = []
step = pi/L
x = -pi
t = -pi 
while (x <= pi) and (t <= pi):
    # f(t) -  функция под знаком интеграла
    def f(t, x):
        if (t >= 0) and (t <= pi):
            if (sin(t - x) == 0):
              return ((1/pi)*(((pi-t)/2)*(N+1)) / 2)
            else:
               return ((1/pi)*((pi-t)/2) * sin((N + 0.5)*(t - x))) / (2 * pi * sin((t - x) / 2))
        elif (t >= - pi) and (t < 0):
            if (sin(t - x) == 0):
              return (((1/pi)*((-pi-t)/2)*(N+1)) / 2)
            else:
               return ((1/pi)*(-pi-t)/2 * sin((N + 0.5)*(t - x))) / (2 * pi * sin((t - x) / 2))
    t += step
    x += step 
 
# по формуле трапеции считаем интеграл integral
# значение SN(f,x) - y
while(x <= pi):
    while (t <= pi):
        integral = (f(-pi, -x) + f(pi, x))/2
        while (t <= pi):
            if (t <= 0):
                integral += f(t, -x)
            else:
                integral += f(t, x)
            t += step
        integral *= step
        y = integral
        ArgF.append(x)
        SND.append(y)
        y0 = pi * x - x * x
        ResF.append(y0)
        t += step
    x += step
    
    
 
print('')
print('Значения функции')
print(ArgF)
print(ResF)
print('')
print('Значения  частичной суммы')
print(SND)
plt.plot(ArgF, ResF)
plt.plot(ArgF, SND)
plt.grid()
plt.show()