[seqfan] Re: A008583

[israel at math.ubc.ca]

The Maple program is not a program, it is an expression. If the generating 
function is g(x), this expression is g(x^2) (I guess that's what "aerated 
with zeros" is supposed to mean, but I have no idea why anyone would want 
to do that). Here's a Maple program for getting the sequence up to a(n):

A008583:= proc(n) local G,j;
  G:= 
series(1/((1-x)*(1-x^3)*(1-x^4)*(1-x^5)*(1-x^6)*(1-x^7)*(1-x^9)),x,n+1);
  [seq(coeff(G,x,j),j=0..n)];
  end proc;

For example:

> A008583(100);
[1, 1, 1, 2, 3, 4, 6, 8, 10, 14, 18, 22, 29, 36, 44, 55, 67, 80, 98, 117, 
138, 165, 194, 226, 266, 309, 356, 413, 475, 542, 622, 708, 802, 911, 1029, 
1157, 1304, 1462, 1633, 1827, 2036, 2261, 2514, 2785, 3077, 3402, 3750, 
4123, 4536, 4977, 5449, 5967, 6520, 7110, 7755, 8441, 9172, 9966, 10810, 
11707, 12678, 13707, 14799, 15976, 17222, 18541, 19958, 21455, 23038, 
24732, 26520, 28406, 30420, 32541, 34776, 37155, 39658, 42291, 45087, 
48024, 51110, 54379, 57809, 61407, 65212, 69198, 73375, 77782, 82395, 
87222, 92307, 97622, 103178, 109020, 115121, 121491, 128179, 135155, 
142432, 150060, 158010]

Robert Israel
University of British Columbia

On Mar 26 2012, Harvey P. Dale wrote:

>            Is the Maple program for this sequence correct?
>
> 
>
>            Best,
>
>            Harvey
>
>
>_______________________________________________
>
>Seqfan Mailing list - http://list.seqfan.eu/
>