(MAGMAMagma) Z<x>:=PolynomialRing(Integers()); N<r>:=NumberField(x^2-3); S:=[ ((3+r)*(4+r)^n+(3-r)*(4-r)^n)/6: n in [0..20] ]; [ Integers()!S[j]: j in [1..#S] ]; // Klaus Brockhaus, Jul 13 2009
(MAGMAMagma) I:=[1, 5]; [n le 2 select I[n] else 8*Self(n-1)-13*Self(n-2): n in [1..40]]; // Vincenzo Librandi, Aug 30 2016