import numpy as np
import matplotlib.pyplot as plt
 
# フーリエ係数
a = [0, 2/np.pi**2, 0, 2/(9*np.pi**2), 0]
b = [2/np.pi, 4/np.pi**2, 2/(3*np.pi), 4/(9*np.pi**2), 2/(5*np.pi)]
 
# 関数を定義
def f(x, n):
    result = np.pi/2  # 直流成分
    for i in range(len(n)):
        result += np.sqrt(a[i]**2 + b[i]**2) * np.sin(n[i] * x)
    return result
 
# x軸の値
x_values = np.linspace(-np.pi, np.pi, 1000)
 
# 関数の値を計算
y_values = f(x_values, [1, 2, 3, 4, 5])
 
# グラフを描画
plt.figure(figsize=(8, 6))
plt.plot(x_values, y_values, label='f(x)')
plt.xlabel('x')
plt.ylabel('f(x)')
plt.title('Graph of f(x)')
plt.grid(True)
plt.legend()
plt.show()# your code goes here

aW1wb3J0IG51bXB5IGFzIG5wCmltcG9ydCBtYXRwbG90bGliLnB5cGxvdCBhcyBwbHQKCiMg44OV44O844Oq44Ko5L+C5pWwCmEgPSBbMCwgMi9ucC5waSoqMiwgMCwgMi8oOSpucC5waSoqMiksIDBdCmIgPSBbMi9ucC5waSwgNC9ucC5waSoqMiwgMi8oMypucC5waSksIDQvKDkqbnAucGkqKjIpLCAyLyg1Km5wLnBpKV0KCiMg6Zai5pWw44KS5a6a576pCmRlZiBmKHgsIG4pOgogICAgcmVzdWx0ID0gbnAucGkvMiAgIyDnm7TmtYHmiJDliIYKICAgIGZvciBpIGluIHJhbmdlKGxlbihuKSk6CiAgICAgICAgcmVzdWx0ICs9IG5wLnNxcnQoYVtpXSoqMiArIGJbaV0qKjIpICogbnAuc2luKG5baV0gKiB4KQogICAgcmV0dXJuIHJlc3VsdAoKIyB46Lu444Gu5YCkCnhfdmFsdWVzID0gbnAubGluc3BhY2UoLW5wLnBpLCBucC5waSwgMTAwMCkKCiMg6Zai5pWw44Gu5YCk44KS6KiI566XCnlfdmFsdWVzID0gZih4X3ZhbHVlcywgWzEsIDIsIDMsIDQsIDVdKQoKIyDjgrDjg6njg5XjgpLmj4/nlLsKcGx0LmZpZ3VyZShmaWdzaXplPSg4LCA2KSkKcGx0LnBsb3QoeF92YWx1ZXMsIHlfdmFsdWVzLCBsYWJlbD0nZih4KScpCnBsdC54bGFiZWwoJ3gnKQpwbHQueWxhYmVsKCdmKHgpJykKcGx0LnRpdGxlKCdHcmFwaCBvZiBmKHgpJykKcGx0LmdyaWQoVHJ1ZSkKcGx0LmxlZ2VuZCgpCnBsdC5zaG93KCkjIHlvdXIgY29kZSBnb2VzIGhlcmU=