これの続き。無限級数の一般項を外部関数化して与えると次のようになる。
import numpy as np
from matplotlib import pyplot as plt
def func(n):
return ((-1) ** (n - 1)) / n
N = 100
f = np.zeros(N)
for n in range(1, N):
f[n] = f[n - 1] + func(n)
print(f[n], np.log(2))
plt.figure()
plt.plot(f)
plt.show()
これで、funcのみを変更すれば好きな級数を試せるようになった。
コメント