A321118 -id:A321118 - cp's OEIS frontend

A321119 a(n) = ((1 - sqrt(3))^n + (1 + sqrt(3))^n)/2^floor((n - 1)/2); n-th row common denominator of A321118.

4, 2, 8, 10, 28, 38, 104, 142, 388, 530, 1448, 1978, 5404, 7382, 20168, 27550, 75268, 102818, 280904, 383722, 1048348, 1432070, 3912488, 5344558, 14601604, 19946162, 54493928, 74440090, 203374108, 277814198, 759002504, 1036816702, 2832635908, 3869452610

Offset: 0

Franck Maminirina Ramaharo, Nov 01 2018

			a(0) = ((1 - sqrt(3))^0 + (1 + sqrt(3))^0)/2^floor((0 - 1)/2) = 2*(1 + 1) = 4.

Harold J. Ahlberg, Edwin N. Nilson and Joseph L. Walsh, The Theory of Splines and Their Applications, Academic Press, 1967. See p. 47, Table 2.5.2.

Encyclopedia of Mathematics, Quadrature formula
John C. Holladay, A smoothest curve approximation, Math. Comp. Vol. 11 (1957), 233-243.
Peter Köhler, On the weights of Sard's quadrature formulas, CALCOLO Vol. 25 (1988), 169-186.
Leroy F. Meyers and Arthur Sard, Best approximate integration formulas, J. Math. Phys. Vol. 29 (1950), 118-123.
Arthur Sard, Best approximate integration formulas; best approximation formulas, American Journal of Mathematics Vol. 71 (1949), 80-91.
Isaac J. Schoenberg, Spline interpolation and best quadrature formulae, Bull. Amer. Math. Soc. Vol. 70 (1964), 143-148.
Frans Schurer, On natural cubic splines, with an application to numerical integration formulae, EUT report. WSK, Dept. of Mathematics and Computing Science Vol. 70-WSK-04 (1970), 1-32.
Index entries for linear recurrences with constant coefficients, signature (0,4,0,-1).

Cf. A002176 (common denominators of Cotesian numbers).
Cf. A321118, A321120.
Cf. A001834, A002530, A002531, A003500, A005246, A048788, A083336, A125905.

LinearRecurrence[{0, 4, 0, -1}, {4, 2, 8, 10}, 50]

a(n) := ((1 - sqrt(3))^n + (1 + sqrt(3))^n)/2^floor((n - 1)/2)$
makelist(ratsimp(a(n)), n, 0, 50);

a(n) = (((sqrt(2) - sqrt(6))/2)^n + ((sqrt(6) + sqrt(2))/2)^n)*((2 - sqrt(2))*(-1)^n + 2 + sqrt(2))/2.
a(-n) = (-1)^n*a(n).
a(n) = 2*A000034(n+1)*A002531(n).
a(2*n) = 2*A001834(n).
a(2*n+1) = 2*A003500(n).
a(n) = 4*a(n-2) - a(n-4) with a(0) = 4, a(1) = 2, a(2) = 8, a(3) = 10.
a(2*n+3) = a(2*n+1) + a(2*n+2).
a(2*n+2) = a(2*n) + 2*a(2*n+1).
G.f.: 2*(1 - x)*(2 + 3*x - x^2)/(1 - 4*x^2 + x^4).
E.g.f.: (1 + exp(-sqrt(6)*x))*((2 - sqrt(2))*exp(sqrt(2 - sqrt(3))*x) + (2 + sqrt(2))*exp(sqrt(2 + sqrt(3))*x))/2.
Lim_{n->infinity} a(2*n+1)/a(2*n) = (1 + sqrt(3))/2.

A321121 Triangle read by rows: T(n,k) is the unreduced numerator of the k-th weight in the quadrature rule for parabolic runout spline with respect to a mesh of n + 1 points.

Original entry on oeis.org

0, 1, 1, 1, 4, 1, 3, 9, 9, 3, 13, 44, 30, 44, 13, 35, 115, 90, 90, 115, 35, 16, 53, 40, 46, 40, 53, 16, 131, 433, 330, 366, 366, 330, 433, 131, 179, 592, 450, 504, 486, 504, 450, 592, 179, 163, 539, 410, 458, 446, 446, 458, 410, 539, 163, 668, 2209, 1680, 1878, 1824, 1842, 1824, 1878, 1680, 2209, 668

Offset: 0

Franck Maminirina Ramaharo, Nov 16 2018

The weights in this quadrature rule are T(n,k)/A321122(n), 0 <= k <= n. For n = 1, 2, 3, we obtain the trapezoid rule, Simpson's rule, and Simpson's 3/8 rule, respectively.

			Triangle begins (denominator is factored out):
    0;                                                 1/4
    1,   1;                                            1/2
    1,   4,   1;                                       1/3
    3,   9,   9,   3;                                  1/8
   13,  44,  30,  44,  13;                             1/36
   35, 115,  90,  90, 115,  35;                        1/96
   16,  53,  40,  46,  40,  53,  16;                   1/44
  131, 433, 330, 366, 366, 330, 433, 131;              1/360
  179, 592, 450, 504, 486, 504, 450, 592, 179;         1/492
  163, 539, 410, 458, 446, 446, 458, 410, 539, 163;    1/448
  ...

Harold J. Ahlberg, Edwin N. Nilson and Joseph L. Walsh, The Theory of Splines and Their Applications, Academic Press, 1967. See p. 47, Table 2.5.3.

Franck Maminirina Ramaharo, Rows n = 0..150 of triangle, flattened
Harold J. Ahlberg, Edwin N. Nilson and Joseph L. Walsh, Chapter II. The Cubic Spline, Mathematics in Science and Engineering Volume 38 (1967), pp. 9-74.
Wikipedia, Newton-Cotes formulas

Cf. A321122 (Common denominators).

Cf. A093735/A093736 (Newton-Cotes formulas), A100640/A100641 (Cotesian numbers), A321118/A321119 (Holladay-Sard best quadrature formulas).

s = -2 + Sqrt[3];
e[n_] := s*(2 + s)*(-1 + s^n)/(2*(1 - s)*(-s + s^n));
f[n_, k_] := 6*s^(1 - k)*(s^(2*k) + s^n)/((1 - s)*(-s + s^n));
w[n_, k_] := If[k == 0 || k == n, 1/4 + e[n]/6, If[k == 1 || k == n - 1, 2 - (1 + 1/6)*e[n], 1 + f[n, k]/4]];
a321122[n_] := LCM @@ Table[Denominator[FullSimplify[w[n, k]]], {k, 0, n}]
Join[{0, 1, 1, 1, 4, 1}, Table[FullSimplify[a321122[n]*w[n, k]], {n, 3, 12}, {k, 0, n}]] // Flatten

s : -2 + sqrt(3)$
e(n) := s*(2 + s)*(-1 + s^n)/(2*(1 - s)*(-s + s^n))$
f(n, k) := 6*s^(1 - k)*(s^(2*k) + s^n)/((1 - s)*(-s + s^n))$
w(n, k) := if k = 0 or k = n then 1/4 + e(n)/6 else if k = 1 or k = n - 1 then 2 - (1 + 1/6)*e(n) else 1 + f(n, k)/4$
a321122(n) := lcm(makelist(denom(fullratsimp(w(n, k))), k, 0, n))$
append([0, 1, 1, 1, 4, 1], create_list(fullratsimp(a321122(n)*w(n, k)), n, 3, 12, k, 0, n));

T(n,k) = T(n,n-k).
T(0,0) = 0 and T(n,k) = A093735(n,k) for n = 1, 2, 3.
Let s = -2 + sqrt(3), and define e(n) = s*(2 + s)*(-1 + s^n)/(2*(1 - s)*(-s + s^n)), f(n,k) = 6*s^(1 - k)*(s^(2*k) + s^n)/((1 - s)*(-s + s^n)), and w(n,0) = 1/4 + e(n)/6, w(n,1) = 2 - (1 + 1/6)*e(n), w(n,k) = 1 + f(n,k)/4 for 2 <= k <= n - 2. Then T(n,k) = A321122(n)*w(n,k) for 0 <= k <= n, n >= 3.

A321120 Decimal expansion of (3 + sqrt(3))/12.

Original entry on oeis.org

3, 9, 4, 3, 3, 7, 5, 6, 7, 2, 9, 7, 4, 0, 6, 4, 4, 1, 1, 2, 7, 2, 8, 7, 1, 9, 5, 1, 2, 5, 4, 8, 9, 3, 6, 3, 9, 1, 1, 9, 0, 0, 4, 3, 7, 8, 1, 7, 5, 3, 1, 7, 1, 9, 0, 0, 4, 6, 5, 0, 5, 8, 1, 6, 2, 0, 9, 9, 4, 4, 1, 8, 0, 7, 5, 7, 3, 3, 3, 3, 6, 4, 2, 3, 4, 2, 8

Offset: 0

Franck Maminirina Ramaharo, Nov 09 2018

The smallest weight in Holladay-Sard's quadrature formula for semi-infinite integrals.

			0.3943375672974064411272871951...

Harold J. Ahlberg, Edwin N. Nilson and Joseph L. Walsh, The Theory of Splines and Their Applications, Academic Press, 1967.

John C. Holladay, A smoothest curve approximation, Math. Comp. Vol. 11 (1957), 233-243.
Leroy F. Meyers and Arthur Sard, Best approximate integration formulas, J. Math. Phys. Vol. 29 (1950), 118-123.
Arthur Sard, Best approximate integration formulas; best approximation formulas, American Journal of Mathematics Vol. 71 (1949), 80-91.
Frans Schurer, On natural cubic splines, with an application to numerical integration formulae, EUT report. WSK, Dept. of Mathematics and Computing Science Vol. 70-WSK-04 (1970), 1-32.

Cf. A020805, A165663.

Cf. A321118, A321119.

Digits := 1000; evalf((3 + sqrt(3))/12);

RealDigits[(3 + Sqrt[3])/12, 10, 100][[1]]

PARI
```
(3 + sqrt(3))/12
```

Equals lim_{n->infinity} A321118(0,n)/A321119(n).
Irrational number represented by the periodic continued fraction [0, 2, 1, 1; [6, 2]].
Largest real root of 1 - 12*x + 24*x^2.

cp's OEIS Frontend

A321119 a(n) = ((1 - sqrt(3))^n + (1 + sqrt(3))^n)/2^floor((n - 1)/2); n-th row common denominator of A321118.

Views

Author

Keywords

Examples

References

Links

Crossrefs

Programs

Mathematica

Maxima

Formula

A321121 Triangle read by rows: T(n,k) is the unreduced numerator of the k-th weight in the quadrature rule for parabolic runout spline with respect to a mesh of n + 1 points.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Maxima

Formula

A321120 Decimal expansion of (3 + sqrt(3))/12.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula