Pythonで書くフーリエ級数の検算用コード - Qiita

import numpy as np import matplotlib.pyplot as plt # x軸を定義する x = np.arange(-np.pi,np.pi,0.01) # シグマの中身の関数 def func(k): return -1*(((-1)**k)/k)*np.sin(k*x) # シグマ演算の関数. 引数で和の範囲を指定. def sigma(func,frm,to): ret = np.zeros_like(x) for i in range(frm,to): ret += func(i) return ret # フーリエ級数展開の関数 f = sigma(func,1,100) # 結果の表示 fig, ax = plt.subplots(1,1) ax.set_aspect('equal') ax.plot(x,f) plt.show() # こっち