Et al, I do not agree with the generating function in A008651 as stated in the Maple program line. What I do get is (1 + x - x^3 - x^4 - x^5 + x^7 + x^8)/((-1 + x)^2 (1 + 3x + 5x^2 + 6x^3 + 6x^4 + 5x^5 + 3x^6 + x^7)). Bob.