%I #20 Sep 08 2022 08:45:18
%S 1,2,6,15,38,94,231,564,1372,3329,8064,19512,47177,114010,275430,
%T 665247,1606534,3879302,9366735,22615356,54601628,131825377,318263328,
%U 768369744,1855031473,4478479058,10812064614,26102729679,63017720390
%N A Fibonacci-Pell convolution.
%C Diagonal sums of A106513.
%H G. C. Greubel, <a href="/A106515/b106515.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (3,0,-3,-1).
%F G.f.: (1-x)/((1-x-x^2)*(1-2*x-x^2)).
%F a(n) = Sum_{k=0..n} Fibonacci(n-k-1)*Pell(k+1).
%F a(n) = Sum_{k=0..floor(n/2)} Sum_{j=0..floor((n-k+1)/2)} binomial(n-k+1, 2*j+k+1)*2^j.
%F a(n) = Pell(n) + Pell(n+1) - Fibonacci(n). - _Ralf Stephan_, Jun 02 2007
%F a(n) = 3*a(n-1) - 3*a(n-3) - a(n-4). - _Wesley Ivan Hurt_, May 27 2021
%t Table[Fibonacci[n, 2] + Fibonacci[n+1, 2] - Fibonacci[n], {n, 0, 30}] (* _Vladimir Reshetnikov_, Sep 27 2016 *)
%o (Magma)
%o Pell:= func< n | Round(((1+Sqrt(2))^n - (1-Sqrt(2))^n)/(2*Sqrt(2))) >;
%o [Pell(n) + Pell(n+1) - Fibonacci(n): n in [0..30]]; // _G. C. Greubel_, Aug 05 2021
%o (Sage) [lucas_number1(n+1, 2, -1) + lucas_number1(n, 2, -1) - lucas_number1(n, 1, -1) for n in (0..30)] # _G. C. Greubel_, Aug 05 2021
%Y Cf. A000045, A000129, A106513.
%K easy,nonn
%O 0,2
%A _Paul Barry_, May 05 2005