Построение график частичных сумм ряда Фурье для двух функций на Python

import numpy as np
import matplotlib.pyplot as plt
from scipy import integrate 
from math import pi
 
f = lambda x: abs(x)
 
print("Введите количество членов частичной суммы ряда")
n = int(input())
a0 = integrate.quad(f, -1*pi, pi)
F = a0/2.00
for i in range (1, n):
    an = integrate.quad(f*np.cos(x*n), -1*pi, pi)
    bn = integrate.quaf(f*np.sin(x*n), -1*pi, pi)
    F += an*np.cos(n*x)+bn*np.sin(n*x)
plt.axis([-4*pi, 4*pi, -2, 2])
plt.title('$S(x)$')
plt.plot(x, F)
plt.show

import numpy as np
import matplotlib.pyplot as plt
N=1024 # число точек
# интегрирование производится методом прямоугольников
x=np.linspace(-np.pi, np.pi, N) 
h=2*np.pi/N # шаг
f=np.abs(x)
n = int(input("Введите количество членов частичной суммы ряда: \n"))
a0 = h*np.sum(f)/np.pi
F = a0/2.00
for i in range (1, n):
    an = h*np.sum(f*np.cos(x*i))/np.pi
    bn = h*np.sum(f*np.sin(x*i))/np.pi
    F += an*np.cos(i*x)+bn*np.sin(i*x)
plt.axis([-np.pi, np.pi, -0.2, 4])
plt.title('$S(x)$')
plt.plot(x, F) #  график частичной суммы
plt.plot(x, f)  # график функции
plt.show

F = a0/2.00*np.ones_like(x)

import numpy as np
import matplotlib.pyplot as plt
N=1024 # число точек
# интегрирование производится методом прямоугольников
x=np.linspace(-np.pi, np.pi, N) 
x_extended=np.linspace(-4*np.pi, 4*np.pi, 4*N)
h=2*np.pi/N # шаг
f=np.abs(x)
n = int(input("Введите количество членов частичной суммы ряда: \n"))
a0 = h*np.sum(f)/np.pi
F = a0/2.00*np.ones_like(x_extended)
for i in range (1, n):
    an = h*np.sum(f*np.cos(x*i))/np.pi
    bn = h*np.sum(f*np.sin(x*i))/np.pi
    F += an*np.cos(i*x_extended)+bn*np.sin(i*x_extended)
plt.axis([-4*np.pi, 4*np.pi, -0.2, 4])
plt.title('$S(x)$')
plt.plot(x_extended, F) #  график частичной суммы
plt.plot(x, f)  # график функции
plt.show