Как разобраться, чтобы преобразовать «сумму продуктов» в «рациональное выражение» p (x) / q (x)

Question

Я пытаюсь получить симпатию к производимому и выводимому, который является частным двух полиномов: p / q.

import sympy as sp
sp.init_printing()

a,b,c = sp.symbols("a b c")

N=a*b*100 - (a**2) * (b**2)
D=2*(a-b)

V = N / D
print(V)
#Output: (-a**2*b**2 + 100*a*b)/(2*a - 2*b)

# HERE"s WHERE I GET THE RESULT:
g = V.diff(a)
print(g)
#Output: (-2*a*b**2 + 100*b)/(2*a - 2*b)   
             -2*(-a**2*b**2 + 100*a*b)/(2*a - 2*b)**2
# The problem here is that its the sum of two terms

# So I try simplifying it to get it as p/q
h= g.simplify()
print(h)
# Output: b*(a*(a*b - 100) + 2*(a - b)*(-a*b + 50)) / (2*(a - b)**2)
#
# It works to get the function as "p/q", except now
# it didn't expand the numerator and denominator into 
# into a sum of polynomial terms.  how to undo the factoring
# of numerator and denominator while still maintaining the whole
# function as a rational function of the form p/q? 

Что я хочу, чтобы это выглядело так:

(- 2 * a 2 * b 2 - 4 * a * b 3) / (4a + 4b 2)

Answer 1

Вы можете использовать ключевое слово numer = True или denom = True с расширением:

>>> q
(x + 1)*(x + 2)/((x - 2)*(x - 1))

>>> q.expand(numer=True)
(x**2 + 3*x + 2)/((x - 2)*(x - 1))

>>> _.expand(denom=True)
(x**2 + 3*x + 2)/(x**2 - 3*x + 2)

Как исправить выше:

import sympy as sp
sp.init_printing()

a,b,c = sp.symbols("a b c")

N=a*b*100 - (a**2) * (b**2)
D=2*(a-b)

N / D
_.diff(a)
_.simplify()
_.expand(numer=True)
_.expand(denom=True)
V = _

Другим способом является использование метода cancel (), который в основном делает то же самое (хотя название несколько нелогично):

import sympy as sp
sp.init_printing()

a,b,c = sp.symbols("a b c")

N=a*b*100 - (a**2) * (b**2)
D=2*(a+b)
V = N / D
V.diff(a).cancel()