(Magma) Z<x>:=PolynomialRing(Integers()); N<r>:=NumberField(x^2-7); S:=[ ((3+r)^n-(3-r)^n)/(2*r): n in [1..21] ]; [ Integers()!S[j]: j in [1..#S] ]; // Klaus Brockhaus, Jan 07 2009
(Magma) I:=[1, 6]; [n le 2 select I[n] else 6*Self(n-1)-2*Self(n-2): n in [1..50]]; // Vincenzo Librandi, Feb 02 2012
(SageMath) [lucas_number1(n, 6, 2) for n in range(1, 22)] # Zerinvary Lajos, Apr 22 2009
(Maxima) a[1]:1$ a[2]:6$ a[n]:=6*a[n-1]-2*a[n-2]$ makelist(a[n], n, 1, 21); /* Bruno Berselli, May 30 2011 */