Processing math: 100%

First time here? Check out the FAQ!

Ask Your Question
1

Coefficients of infinite polynomial products

asked 11 years ago

Blackadder gravatar image

updated 11 years ago

I would first like to ask if there are commands for computing infinite products,

eg.g=qn=1(1q8n)(1q16n).

If not, are there commands for finite products then? How do you compute g=qNn=1(1q8n)(1q16n)?

Basically, I'm only interested in the coefficients of gθ2 and gθ4 where g is as above and θt=qtn2.

This link to OEIS seems to have a code for generating gθ2 in PARI, I always thought that SAGE is able to call pari, but I'm unable to perform the command on SAGE.

I'll appreciate it if anyone can advise me how best to do this or direct me to any resources, thanks!

Preview: (hide)

1 Answer

Sort by » oldest newest most voted
1

answered 11 years ago

Luca gravatar image

For finite products you can use the function prod(), like this:

sage: prod(range(1,42))
33452526613163807108170062053440751665152000000000L

What do you mean by "compute K number of times"?

There certainly is not a generic function to compute infinite product, but there might be functions to do what you have in mind. Have a look at the documentation of eta and EtaProduct, you can use the ? syntax:

sage: eta?
...
sage: EtaProduct?

unfortunately, eta seems to accept only complex arguments.

You can use the function gp() to enter Pari code

sage: a = gp("a(n)=local(A); if(n<0, 0, A=x*O(x^n); polcoeff(eta(x^2+A)^5*eta(x^8+A)/(eta(x+A)^2*eta(x^4+A)), n))")
sage: [a(i) for i in range(10)]
[1, 2, 0, 0, 1, -2, 0, 0, -4, -2]
Preview: (hide)
link

Comments

Hi Luca, thanks for answering. This is what I was looking for. As to your question, I've edited what I meant. Hopefully that would clear up any confusion.

Blackadder gravatar imageBlackadder ( 11 years ago )

Your Answer

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

Add Answer

Question Tools

1 follower

Stats

Asked: 11 years ago

Seen: 1,256 times

Last updated: Aug 05 '13