python バーゼル級数が出力されない

Question

前提・実現したいこと

大学の課題でバーゼル級数のプログラムを組んでいた

発生している問題・エラーメッセージ

print(b(10))としても１が出力される。

import numpy as np

def b(N):
    x = np.arange(1,N+1)
    a = np.sum(1/(x**2))
    return a
print(str(b(10)))

ソースコード


### 試したこと


### 補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

Comments

[meg_]

「1.5497677311665408」になりませんか？

[Etori]

解決しました！ありがとうございました！

Answer 1

python のバージョンが書かれていませんが、python2 をお使いですよね？ python2 では、整数を整数で除算した結果は整数になります。小数点以下は切り捨てられます。
例：

> print(5/3)
1

したがって、1/(x**2)の結果もすべて切り捨てられた整数となり、最初の要素は１でそれ以外は０になります。
したがって、その和は１となります。
被除数を浮動小数点にすることで、浮動小数点で演算できます。以下を試してください。

import numpy as np

def b(N):
    x = np.arange(1,N+1)
    a = np.sum(1.0/(x**2))
    return a
print(str(b(10)))

ちなみに python3 であれば、被除数が整数の場合でも結果は浮動小数点数になります。

Comments

[Etori]

解決しました！　ありがとうございます！