Trouble when computing A006743

Thomas Baruchel baruchel at bluebottle.com
Sat Oct 15 19:44:39 CEST 2005


Hi, Simon Plouffe mentionned the A006743 sequence to me in a separate mail.
I tried to compute it according to the formula, but I can't figure out
what is wrong with my computation. Though I use Maxima syntax, you will
understand very easely this piece of code even if you never used this
software:

block([a6:1,a5:0,a4:3,a3:2,a2:10,a1:14],for n:9 thru 20 do block([a:(-6*
a1+2*a1*n-33*a2+7*a2*n+60*a3-12*a3*n+84*a4-16*a4*n-96*a5+16*n*a5-96*a6+16*n*a6)/
(n-3)],sprint(a,","),a6:a5,a5:a4,a4:a3,a3:a2,a2:a1,a1:a));

output being (for a_9, a_10, ...):
40 , 74 , 177 , 360 , 809 , 1694 , 3709 , 7824 , 16895 , 35650 , 76282 , 160638
,

instead of:
1,0,3,2,10,14,40,74,176,358,798,1670,3626,7638,16366,34462,


-- 
Thomas Baruchel
  Home Page: http://baruchel.free.fr/~thomas/





More information about the SeqFan mailing list