cp's OEIS Frontend

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.

A051028 Ramanujan's a-series: expansion of (1+53x+9x^2)/(1-82x-82x^2+x^3).

Original entry on oeis.org

1, 135, 11161, 926271, 76869289, 6379224759, 529398785665, 43933719985479, 3645969360009049, 302571523160765631, 25109790452983538281, 2083810036074472911735, 172931123203728268135681, 14351199415873371782349831, 1190976620394286129666900249
Offset: 0

Views

Author

Keywords

Comments

The "amazing" identity of Ramanujan is a(n)^3 + b(n)^3 = c(n)^3 + (-1)^n, where a(n)=A051028(n), b(n)=A051029(n) and c(n)=A051030(n). - Emeric Deutsch, Oct 14 2006

Crossrefs

Programs

  • Maple
    g:=(1+53*x+9*x^2)/(1-82*x-82*x^2+x^3): gser:=series(g,x=0,20): seq(coeff(gser,x,n),n=0..12); # Emeric Deutsch, Oct 14 2006
  • Mathematica
    CoefficientList[Series[(1 + 53 x + 9 x^2)/(1 - 82 x - 82 x^2 + x^3), {x, 0, 33}], x] (* Vincenzo Librandi, Jul 22 2015 *)
  • PARI
    Vec((1+53*x+9*x^2)/(1-82*x-82*x^2+x^3) + O(x^30)) \\ Michel Marcus, Feb 29 2016

Formula

G.f.: (1+53*x+9*x^2)/((1+x)*(1-83*x+x^2)).
X(n+1) = A*X(n), where X(n) = transpose(A051028(n), A051029(n), A051030(n)) and A = matrix(3,3,[63,104,-68; 64,104,-67; 80,131,-85]). - Emeric Deutsch, Oct 14 2006

Extensions

Minor edits (g.f. and name) by M. F. Hasler, May 08 2016