Ask Your Question

Finding the coefficient of multivariable generating function

asked 2023-08-10 16:06:32 +0200

jacop gravatar image

I want to find the coefficient of $[a^3b^3c^3]$ in the expansion of $$\bigg[1-\bigg(\frac{a+a^2}{1-(a+a^2)}+\frac{b+b^2}{1-(b+b^2)}+\frac{c+c^2}{1-(c+c^2)}\bigg)\bigg]^{(-1)}$$

I have made research on internet about coefficient of this multivariable generating function, but could not find something valuable. Can you help me to calculate the coefficent using Sagemath ?

edit retag flag offensive close merge delete

1 Answer

Sort by ยป oldest newest most voted

answered 2023-08-10 18:23:14 +0200

Max Alekseyev gravatar image

updated 2023-08-10 18:36:07 +0200

For example:

myprec = 10
def f(z):
    t = (z + z^2).add_bigoh(myprec)
    return t/(1-t)

R.<a,b,c> = PowerSeriesRing(QQ,default_prec=myprec)
g = (1-(f(a)+f(b)+f(c)))^(-1)
print( g.coefficients()[a^3*b^3*c^3] )

For small degrees, default_prec= and .add_bigoh() are not really needed but they will speed up things / save memory when degrees of interest are large.

edit flag offensive delete link more


thanks a lot

jacop gravatar imagejacop ( 2023-08-10 19:39:44 +0200 )edit

@Max Alekseyev , i want to add something. Assume that i want to calculate g = (1-x(f(a)+f(b)+f(c)))^(-1) where print( g.coefficients()[x^9a^3b^3c^3] ). It did not work, how can i handle it ?

jacop gravatar imagejacop ( 2023-08-10 20:49:28 +0200 )edit

You need to define x first. E.g.:

K.<x> = QQ[]
R.<a,b,c> = PowerSeriesRing(K,default_prec=myprec)
g = (1-x*(f(a)+f(b)+f(c)))^(-1)
Max Alekseyev gravatar imageMax Alekseyev ( 2023-08-10 22:43:36 +0200 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools


Asked: 2023-08-10 16:02:26 +0200

Seen: 97 times

Last updated: Aug 10 '23