外部関数で定義された無限級数を計算する

これの続き。無限級数の一般項を外部関数化して与えると次のようになる。

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のみを変更すれば好きな級数を試せるようになった。

コメント

タイトルとURLをコピーしました