This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A133475 #16 Sep 08 2022 08:45:31
%S A133475 -4,-3,-1,0,1,3,5,11,15,28,47,55,81,549,1799,8361
%N A133475 Integers n such that n^3 + n^2 - 9*n + 16 is a square.
%C A133475 The set of x values of integral points on the elliptic curve y^2 = x^3 + x^2 - 9*x + 16.
%e A133475 0^3 + (-5)^2 + (-9) = 4^2, 1^3 + (-4)^2 + (-8) = 3^2, 3^3 + (-2)^2 + (-6) = 5^2
%t A133475 ok[x_] := Reduce[{y^2 == x^3 + x^2 - 9*x + 16, y >= 0}, y, Integers] =!= False; Select[Table[x, {x, -4, 10000}], ok ] (* _Jean-François Alcover_, Sep 07 2011 *)
%o A133475 (Magma) P<n> := PolynomialRing(Integers()); {x: x in Sort([ p[1] : p in IntegralPoints(EllipticCurve(n^3 + n^2 - 9*n + 16)) ])};
%o A133475 (Sage) EllipticCurve([0,1,0,-9,16]).integral_points()
%o A133475 (PARI) is(n)=issquare(n^3+n^2-9*n+16) \\ _Charles R Greathouse IV_, Sep 06 2016
%Y A133475 Cf. A117950, A132411, A132414, A002522, A028872.
%K A133475 sign,full,fini
%O A133475 1,1
%A A133475 _Mohamed Bouhamida_, Nov 29 2007
%E A133475 Edited by _Max Alekseyev_, Nov 13 2009