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.

Showing 1-2 of 2 results.

A162546 A Somos-4 variant: a(n) = (36*a(n-1)*a(n-3) - 68*a(n-2)^2)/a(n-4).

Original entry on oeis.org

1, 1, -16, -644, -40592, -4821056, 17059328, 2492895195136, 10659285907800064, 86296767700623425536, 1081586547380924161458176, -36649408809924048998874742784, -18144416387824430577315746611724288
Offset: 0

Views

Author

Paul Barry, Jul 05 2009

Keywords

Comments

Hankel transform of A162543, A162548.

Links

Crossrefs

Cf. A157005, A157101, A162547.

Programs

  • GAP
    a:=[1,1,-16,-644];; for n in [5..20] do a[n]:=(36*a[n-1]*a[n-3] - 68*a[n-2]^2)/a[n-4]; od; a; # G. C. Greubel, Feb 23 2019
  • Magma
    I:=[1,1,-16,-644]; [n le 4 select I[n] else (36*Self(n-1) *Self(n-3) - 68*Self(n-2)^2)/Self(n-4): n in [1..20]]; // G. C. Greubel, Feb 23 2019
  • Mathematica
    RecurrenceTable[{a[n]==(36*a[n-1]*a[n-3] - 68*a[n-2]^2)/a[n-4], a[0]==1, a[1]==1, a[2]==-16, a[3]==-644}, a, {n,20}] (* G. C. Greubel, Feb 23 2019 *)
  • PARI
    m=20; v=concat([1,1,-16,-644], vector(m-4)); for(n=5, m, v[n] = (36*v[n-1]*v[n-3] -68*v[n-2]^2)/v[n-4]); v \\ G. C. Greubel, Feb 23 2019
  • Sage
    def a(n):
    if (n==0): return 1
    elif (n==1): return 1
    elif (n==2): return -16
    elif (n==3): return -644
    else: return (36*a(n-1)*a(n-3) - 68*a(n-2)^2)/a(n-4)
[a(n) for n in (1..20)] # G. C. Greubel, Feb 23 2019

A162548 A Chebyshev transform of the little Schroeder numbers A001003.

Original entry on oeis.org

1, 1, 2, 9, 37, 156, 695, 3203, 15118, 72739, 355475, 1759624, 8804341, 44457125, 226256114, 1159387253, 5976713665, 30974296468, 161285018771, 843388543471, 4427120165182, 23319497761799, 123221525405447, 652989260163472
Offset: 0

Views

Author

Paul Barry, Jul 05 2009

Keywords

Comments

Hankel transform is Somos-4 variant A162546.

Links

Crossrefs

Cf. A001003, A162543.

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( (1+x+x^2 -Sqrt(1-6*x+3*x^2-6*x^3+x^4))/(4*x*(1+x^2)) )); // G. C. Greubel, Feb 26 2019
  • Mathematica
    CoefficientList[Series[(1+x+x^2 -Sqrt[1-6*x+3*x^2-6*x^3+x^4])/(4*x*(1+x^2)), {x,0,30}], x] (* G. C. Greubel, Feb 26 2019 *)
  • PARI
    my(x='x+O('x^30)); Vec((1+x+x^2 -sqrt(1-6*x+3*x^2-6*x^3+x^4))/( 4*x*(1+x^2))) \\ G. C. Greubel, Feb 26 2019
  • Sage
    ((1+x+x^2 -sqrt(1-6*x+3*x^2-6*x^3+x^4))/(4*x*(1+x^2))).series(x, 30).coefficients(x, sparse=False) # G. C. Greubel, Feb 26 2019

Formula

G.f.: (1/(1+x^2))*s(x/(1+x^2)), s(x) the g.f. of A001003.
G.f.: (1+x+x^2 - sqrt(1-6*x+3*x^2-6*x^3+x^4))/(4*x*(1+x^2)).
G.f.: 1/(1+x^2-x/(1-2*x/(1+x^2-x/(1-2*x/(1+x^2-x/(1-2*x/(1+x+x^2/(1-... (continued fraction);
a(n) = Sum_{k=0..floor(n/2)} (-1)^k*C(n-k,k)*A001003(n-2k).
Conjecture: (n+1)*a(n) +3*(-2*n+1)*a(n-1) +(4*n-5)*a(n-2) +12*(-n+2)*a(n-3) +(4*n-11)*a(n-4) +3*(-2*n+7)*a(n-5) +(n-5)*a(n-6)=0. - R. J. Mathar, Nov 15 2012. (Formula verified and used for computations. - Fung Lam, Feb 19 2014)
Showing 1-2 of 2 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.