|
PROGRAM
|
(PARI) /* Generate using matrix power method: */ {a(n)=local(A=Mat(1), B); for(m=1, n+3, B=matrix(m, m); for(i=1, m, for(j=1, i, if(j==i, B[i, j]=1, B[i, j]=(A^(fibonacci(i-1)))[i-1, j]); )); A=B); return( ((A)[n+3, 3]))} (PARI) /* Generate using partial sums method: */ {a(n)=local(A=vector(n+1), p); A[1]=1; for(j=1, n, p=fibonacci(n+4)-fibonacci(n-j+4)-j; A=Vec((Polrev(A)+x*O(x^p))/(1-x))); A[ #A]}
|
