def test_gosper_normal():
eq = 4 * n + 5, 2 * (4 * n + 1) * (2 * n + 3), n
assert gosper_normal(*eq) == (
Poly(Rational(1, 4), n),
Poly(n + Rational(3, 2)),
Poly(n + Rational(1, 4)),
)
assert gosper_normal(*eq, polys=False) == (
Rational(1, 4),
n + Rational(3, 2),
n + Rational(1, 4),
)
def test_gosper_normal():
assert gosper_normal(4*n + 5, 2*(4*n + 1)*(2*n + 3), n) == \
(Poly(S(1)/4, n), Poly(n + S(3)/2), Poly(n + S(1)/4))
def test_gosper_normal():
assert gosper_normal(4*n + 5, 2*(4*n + 1)*(2*n + 3), n) == \
(Poly(S(1)/4, n), Poly(n + S(3)/2), Poly(n + S(1)/4))
def test_gosper_normal():
eq = 4 * n + 5, 2 * (4 * n + 1) * (2 * n + 3), n
assert gosper_normal(*eq) == \
(Poly(S(1)/4, n), Poly(n + S(3)/2), Poly(n + S(1)/4))
assert gosper_normal(*eq, polys=False) == \
(S(1)/4, n + S(3)/2, n + S(1)/4)
from sympy import S from sympy import * from sympy.concrete.gosper import gosper_sum, gosper_normal from sympy.functions import factorial from sympy.abc import z, n f1 = z**2 / ((z - 1) * (z - 0.5)) F1_z = gosper_normal(f1, (z, 0, S.Infinity)) print(F1_z)