Real-rooted atlas

A086645 → Veronese 2 → D+1

nonnegative coefficients

Proved for all n all roots real, in (−∞, 0]coefficients of one signHadamard kernel

Recurrence

\(P_{n}(x) = \left(127 - 189\,x - 3\,x^{2} + 65\,x^{3}\right)\,P_{n-6}(x) + \left(-64 + 192\,x - 192\,x^{2} + 64\,x^{3}\right)\,P_{n-6}'(x) + \left(3214 - 988\,x - 1186\,x^{2} - 16\,x^{3}\right)\,P_{n-5}(x) + \left(-128 + 1264\,x - 1120\,x^{2} - 16\,x^{3}\right)\,P_{n-5}'(x) + \left(10193 + 5330\,x + 61\,x^{2}\right)\,P_{n-4}(x) + \left(1152 + 5184\,x + 64\,x^{2}\right)\,P_{n-4}'(x) + \left(-188 + 468\,x + 32\,x^{2}\right)\,P_{n-3}(x) + \left(-2048 + 416\,x + 32\,x^{2}\right)\,P_{n-3}'(x) + \left(1361 + 131\,x\right)\,P_{n-2}(x) + \left(1472 + 128\,x\right)\,P_{n-2}'(x) + \left(-370 - 16\,x\right)\,P_{n-1}(x) + \left(-384 - 16\,x\right)\,P_{n-1}'(x)\)
\(P_{1}(x) = 1\)
\(P_{2}(x) = 1\)
\(P_{3}(x) = 2 + x\)
\(P_{4}(x) = 16 + 15\,x\)
\(P_{5}(x) = 71 + 72\,x + x^{2}\)
\(P_{6}(x) = 211 + 300\,x + 45\,x^{2}\)

Rows ?

nx0x1x2x3x4x5x6x7Pn(1)
111
211
3213
4161531
571721144
621130045556
7496148549812480
81002700732769111376
9182127560183301824149536
103061905769945019176153212416
11484625678550388014535048501915712
1273166468552259708945098757682313934976
131062714815688847939564604078860110632116785408
141495131395003053537530904640646288523218032571289856
1520476623668594373370152108775459622803230955204821301953024
Technical details
Node id
n-6edce57e1ba6e1f4 · revision 0 · JSON
Definition source
derived:one_plus_derivative
Fingerprint
f6-f20248c20d6da8d2

Observed on the stored rows ?

  • Degrees: 0, 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7 (half affine: slope 1, intercept -1)
  • Palindromic rows: no
  • γ-positive (gamma_positive): no
  • Nonnegative coefficients (nonn): yes
  • Only rational roots (rr_by_factors: every nonzero row is a product of linear factors over ℚ): no

Matches ?

No OEIS matches.

Submission text ?

Atlas rows n = 1, …, ascending powers of x
%e Example
Triangle begins:
   1:    1;
   2:    1;
   3:    2,     1;
   4:   16,    15;
   5:   71,    72,     1;
   6:  211,   300,    45;
   7:  496,  1485,   498,     1;
   8: 1002,  7007,  3276,    91;
   9: 1821, 27560, 18330,  1824,   1;
  10: 3061, 90576, 99450, 19176, 153;
%F Formula
The row polynomials satisfy P(n,x) = (127-189*x-3*x^2+65*x^3)*P(n-6,x) - (64-192*x+192*x^2-64*x^3)*d/dx P(n-6,x) + (3214-988*x-1186*x^2-16*x^3)*P(n-5,x) - (128-1264*x+1120*x^2+16*x^3)*d/dx P(n-5,x) + (10193+5330*x+61*x^2)*P(n-4,x) + (1152+5184*x+64*x^2)*d/dx P(n-4,x) - (188-468*x-32*x^2)*P(n-3,x) - (2048-416*x-32*x^2)*d/dx P(n-3,x) + (1361+131*x)*P(n-2,x) + (1472+128*x)*d/dx P(n-2,x) - (370+16*x)*P(n-1,x) - (384+16*x)*d/dx P(n-1,x) for n >= 7, with P(1,x) = 1, P(2,x) = 1, P(3,x) = 2 + x, P(4,x) = 16 + 15*x, P(5,x) = 71 + 72*x + x^2, P(6,x) = 211 + 300*x + 45*x^2.
%t Mathematica
P[1, x_] := 1;
P[2, x_] := 1;
P[3, x_] := 2 + x;
P[4, x_] := 16 + 15*x;
P[5, x_] := 71 + 72*x + x^2;
P[6, x_] := 211 + 300*x + 45*x^2;
P[n_, x_] := P[n, x] = Expand[(127 - 189*x - 3*x^2 + 65*x^3)*P[n - 6, x] - (64 - 192*x + 192*x^2 - 64*x^3)*D[P[n - 6, x], x] + (3214 - 988*x - 1186*x^2 - 16*x^3)*P[n - 5, x] - (128 - 1264*x + 1120*x^2 + 16*x^3)*D[P[n - 5, x], x] + (10193 + 5330*x + 61*x^2)*P[n - 4, x] + (1152 + 5184*x + 64*x^2)*D[P[n - 4, x], x] - (188 - 468*x - 32*x^2)*P[n - 3, x] - (2048 - 416*x - 32*x^2)*D[P[n - 3, x], x] + (1361 + 131*x)*P[n - 2, x] + (1472 + 128*x)*D[P[n - 2, x], x] - (370 + 16*x)*P[n - 1, x] - (384 + 16*x)*D[P[n - 1, x], x]];
Table[CoefficientList[P[n, x], x], {n, 1, 10}] // Flatten
%o PARI/GP
P(N) = {my(v = vector(N)); v[1] = 1; v[2] = 1; v[3] = 2 + x; v[4] = 16 + 15*x; v[5] = 71 + 72*x + x^2; v[6] = 211 + 300*x + 45*x^2; for(n = 7, N, my(i = n); v[i] = (127 - 189*x - 3*x^2 + 65*x^3)*v[i - 6] - (64 - 192*x + 192*x^2 - 64*x^3)*deriv(v[i - 6], x) + (3214 - 988*x - 1186*x^2 - 16*x^3)*v[i - 5] - (128 - 1264*x + 1120*x^2 + 16*x^3)*deriv(v[i - 5], x) + (10193 + 5330*x + 61*x^2)*v[i - 4] + (1152 + 5184*x + 64*x^2)*deriv(v[i - 4], x) - (188 - 468*x - 32*x^2)*v[i - 3] - (2048 - 416*x - 32*x^2)*deriv(v[i - 3], x) + (1361 + 131*x)*v[i - 2] + (1472 + 128*x)*deriv(v[i - 2], x) - (370 + 16*x)*v[i - 1] - (384 + 16*x)*deriv(v[i - 1], x)); v};
row(n, p) = {my(c = Vecrev(p)); c};
T(N) = {my(v = P(N)); concat(vector(#v, i, row(i, v[i])))};
T(10)

polytool c48629374597

Properties

Proof chains ? JSON

References

  1. [ASW52] Michael Aissen, I. J. Schoenberg and A. M. Whitney. On the generating functions of totally positive sequences I. Journal d’Analyse Mathématique, 2(1):93–103, 1952. doi:10.1007/BF02786970
  2. [AW24] Christos A. Athanasiadis and David G. Wagner. Veronese sections and interlacing matrices of polynomials and formal power series. arXiv:2404.12989. arXiv
  3. [Edr53] Albert Edrei. Proof of a Conjecture of Schoenberg on the Generating Function of a Totally Positive Sequence. Canadian Journal of Mathematics, 5:86–94, 1953. doi:10.4153/CJM-1953-010-3
  4. [Fav35] Jean Favard. Sur les polynômes de tchebicheff. C. R. Acad. Sci. Paris, 200:2052–2053, 1935.
  5. [LW07] Lily L. Liu and Yi Wang. A unified approach to polynomial sequences with only real zeros. Advances in Applied Mathematics, 38(4):542–560, May 2007. doi:10.1016/j.aam.2006.02.003
  6. [SP14] J Schur and G Polya. Über zwei Arten von Faktorenfolgen in der Theorie der algebraischen Gleichungen. J. Reine Angew. Math., 1914.
  7. [Sze39] G. Szegő. Orthogonal Polynomials. Colloquium publications, American Mathematical Society, 1939. doi:10.1090/coll/023