テイラー展開を求めるプログラムをかく

Question

前提・実現したいこと

sinx = Σ[k=1,n]{(-1)^(k-1)*x^(2k-1)}/(2k-1)!
のx=0におけるテーラ展開を求めたい。

試したこと

そこでak=(-1)^(k-1)/(2k-1)!としてコードを書きましたが、
テーラ展開された式が出力されません

from sympy import Symbol, pprint
import math

x = Symbol("x")
def ak(n):
ak = 0
for k in range(1, n+1):
ak = ak + x**(2k - 1)/math.factorial(2k - 1);
pprint(ak)
return ak
n = int(input())
ak(n)

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

Answer 1

テイラー展開なら、sympy.series() を使えばよいかと思います。

from sympy import Symbol, init_printing, series, sin

init_printing()

x = Symbol('x')
# x = 0 周りで12次まで展開 (マクローリン展開)
eqn = series(sin(x), x, 0, 12)
eqn

Jupyter Notebook だと Latex 形式で出力されます。

追記

試した質問者さんのコードは以下です。
実行して n=5 としたところ、
x**9/362880 + x**7/5040 + x**5/120 + x**3/6 + x と表示されました。

import math
from sympy import Symbol, pprint

x = Symbol("x")


def ak(n):
    ak = 0
    for k in range(1, n + 1):
        ak = ak + x**(2 * k - 1) / math.factorial(2 * k - 1)
    pprint(ak)
    return ak


n = int(input())
ak(n)

Comments

[emiime]

ありがとうございます。このコードがどうしてうまくいかないのか、教えていただけたら嬉しいです。すみません。

[tiitoi]

どういう結果を期待してうまくいかないと言っているのでしょうか？pprint() で式がでてくると思います。
latex 形式の式は Jupiter Notebook で実行すれば出てきます。
ご提示のコードでも符号が間違っていることを除いては n=6 として同じ式がでてきましたが
----
sinx の展開式
http://w3e.kanazawa-it.ac.jp/math/category/suuretu/maclaurin/henkan-tex.cgi?target=/math/category/suuretu/maclaurin/maclaurin_sinx.html

[emiime]

このコードで実行すると2・xという結果になってしまうのですが。
nをどんな数にしても2・xになってしまって、理由が分からず。
何度も質問申し訳ありません！

[tiitoi]

2xになったというのは質問欄のコードですか？
マークダウン記法を使ってないため、インデントが崩れていて元のコードがわかりませんが、インデントを直して実行したコードを回答欄に書きました。
これを実行しても2xとなるのですか？

[emiime]

回答欄に書いていただいたコードを貼り付けたら正しく実行できました。
インデントでこんなにも変わってしまうのだ、と反省しました。
本当にありがとうございました！！！！