Finding the coefficient of multivariable generating function

Loading [MathJax]/jax/output/HTML-CSS/jax.js

First time here? Check out the FAQ!

jacop
17 ●6

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 ?

add a comment

1 Answer

Sort by » oldest newest most voted

answered 1 year ago

Max Alekseyev

6932 ●9 ●52 ●162

updated 1 year ago

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.

link

Comments

thanks a lot

jacop ( 1 year ago )

@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 ( 1 year ago )

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 ( 1 year ago )

add a comment

Your Answer

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

Add Answer

Question Tools

Stats

Asked: 1 year ago

Seen: 222 times

Last updated: Aug 10 '23

Finding the coefficient of multivariable generating function savecancel

1 Answer

Comments