A098603 -id:A098603 - cp's OEIS frontend

A097480 Positive integers n such that 2n-15 is prime.

9, 10, 11, 13, 14, 16, 17, 19, 22, 23, 26, 28, 29, 31, 34, 37, 38, 41, 43, 44, 47, 49, 52, 56, 58, 59, 61, 62, 64, 71, 73, 76, 77, 82, 83, 86, 89, 91, 94, 97, 98, 103, 104, 106, 107, 113, 119, 121, 122, 124, 127, 128, 133, 136, 139, 142, 143, 146, 148, 149, 154, 161, 163

Offset: 1

Douglas Winston (douglas.winston(AT)srupc.com), Sep 19 2004

Shawn A. Broyles, Table of n, a(n) for n = 1..1000

Cf. A000040, A098602, A098603.
Numbers n such that 2n+k is prime: A005097 (k=1), A067076 (k=3), A089038 (k=5), A105760 (k=7), A155722 (k=9), A101448 (k=11), A153081 (k=13), A089559 (k=15), A173059 (k=17), A153143 (k=19).
Numbers n such that 2n-k is prime: A006254 (k=1), A098090 (k=3), A089253 (k=5), A089192 (k=7), A097069 (k=9), A097338 (k=11), A097363 (k=13), this sequence (k=15), A098605 (k=17), A097932 (k=19).

(Prime[Range[2,100]]+15)/2 (* Vladimir Joseph Stephan Orlovsky, Feb 08 2010 *)
Select[Range[9,200],PrimeQ[2#-15]&] (* Harvey P. Dale, Apr 04 2021 *)

is(n)=isprime(2*n-15) \\ Charles R Greathouse IV, Jul 12 2016

Half of p+15 where p is a prime greater than 2.

A098605 Positive integers n such that 2n - 17 is prime.

Original entry on oeis.org

10, 11, 12, 14, 15, 17, 18, 20, 23, 24, 27, 29, 30, 32, 35, 38, 39, 42, 44, 45, 48, 50, 53, 57, 59, 60, 62, 63, 65, 72, 74, 77, 78, 83, 84, 87, 90, 92, 95, 98, 99, 104, 105, 107, 108, 114, 120, 122, 123, 125, 128, 129, 134, 137, 140, 143, 144, 147, 149, 150, 155, 162

Offset: 1

Douglas Winston (douglas.winston(AT)srupc.com), Sep 20 2004

Shawn A. Broyles, Table of n, a(n) for n = 1..1000

Cf. A000040, A098602, A098603.
Numbers n such that 2n+k is prime: A005097 (k=1), A067076 (k=3), A089038 (k=5), A105760 (k=7), A155722 (k=9), A101448 (k=11), A153081 (k=13), A089559 (k=15), A173059 (k=17), A153143 (k=19).
Numbers n such that 2n-k is prime: A006254 (k=1), A098090 (k=3), A089253 (k=5), A089192 (k=7), A097069 (k=9), A097338 (k=11), A097363 (k=13), A097480 (k=15), this sequence (k=17), A097932 (k=19).

(Prime[Range[2,100]]+17)/2 (* Vladimir Joseph Stephan Orlovsky, Feb 08 2010 *)

is(n)=isprime(2*n-17) \\ Charles R Greathouse IV, Jul 12 2016

Half of p+17 where p is a prime greater than 2.

A158405 Triangle T(n,m) = 1+2*m of odd numbers read along rows, 0<=m

Original entry on oeis.org

1, 1, 3, 1, 3, 5, 1, 3, 5, 7, 1, 3, 5, 7, 9, 1, 3, 5, 7, 9, 11, 1, 3, 5, 7, 9, 11, 13, 1, 3, 5, 7, 9, 11, 13, 15, 1, 3, 5, 7, 9, 11, 13, 15, 17, 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23

Offset: 1

Paul Curtz, Mar 18 2009

Row sums are n^2 = A000290(n).
The triangle sums, see A180662 for their definitions, link this triangle of odd numbers with seventeen different sequences, see the crossrefs. The knight sums Kn14 - Kn110 have been added. - Johannes W. Meijer, Sep 22 2010
A208057 is the eigentriangle of A158405 such that as infinite lower triangular matrices, A158405 * A208057 shifts the latter, deleting the right border of 1's. - Gary W. Adamson, Feb 22 2012
T(n,k) = A099375(n-1,n-k), 1<=k<=n. [Reinhard Zumkeller, Mar 31 2012]

			The triangle contains the first n odd numbers in row n:
  1;
  1,3;
  1,3,5;
  1,3,5,7;
From _Seiichi Manyama_, Dec 02 2017: (Start)
    |       a(n)        |                               | A000290(n)
   -----------------------------------------------------------------
   0|                                                      (=  0)
   1|                 1 = 1/3 * ( 3)                       (=  1)
   2|             1 + 3 = 1/3 * ( 5 +  7)                  (=  4)
   3|         1 + 3 + 5 = 1/3 * ( 7 +  9 + 11)             (=  9)
   4|     1 + 3 + 5 + 7 = 1/3 * ( 9 + 11 + 13 + 15)        (= 16)
   5| 1 + 3 + 5 + 7 + 9 = 1/3 * (11 + 13 + 15 + 17 + 19)   (= 25)
(End)

Seiichi Manyama, Rows n = 1..140, flattened
Daniel Erman, The Josephus Problem, Numberphile video (2016)
Boris Putievskiy, Transformations [of] Integer Sequences And Pairing Functions arXiv:1212.2732 [math.CO], 2012.
Index entries for sequences related to the Josephus Problem

Cf. A129326, A099375, A005408.
Triangle sums (see the comments): A000290 (Row1; Kn11 & Kn4 & Ca1 & Ca4 & Gi1 & Gi4); A000027 (Row2); A005563 (Kn12); A028347 (Kn13); A028560 (Kn14); A028566 (Kn15); A098603 (Kn16); A098847 (Kn17); A098848 (Kn18); A098849 (Kn19); A098850 (Kn110); A000217 (Kn21. Kn22, Kn23, Fi2, Ze2); A000384 (Kn3, Fi1, Ze3); A000212 (Ca2 & Ze4); A000567 (Ca3, Ze1); A011848 (Gi2); A001107 (Gi3). - Johannes W. Meijer, Sep 22 2010
Cf. A063656, A063657, A208057, A292610, A292611.

a158405 n k = a158405_row n !! (k-1)
a158405_row n = a158405_tabl !! (n-1)
a158405_tabl = map reverse a099375_tabl
-- Reinhard Zumkeller, Mar 31 2012

Table[2 Range[1, n] - 1, {n, 12}] // Flatten (* Michael De Vlieger, Oct 01 2015 *)

a(n) = 2*(n-floor((-1+sqrt(8*n-7))/2)*(floor((-1+sqrt(8*n-7))/2)+1)/2)-1;
vector(100, n, a(n)) \\ Altug Alkan, Oct 01 2015

a(n) = 2*i-1, where i = n-t(t+1)/2, t = floor((-1+sqrt(8*n-7))/2). - Boris Putievskiy, Feb 03 2013

a(n) = 2*A002262(n-1) + 1. - Eric Werley, Sep 30 2015

Edited by R. J. Mathar, Oct 06 2009

A120071 a(n) = n*(n+20).

Original entry on oeis.org

0, 21, 44, 69, 96, 125, 156, 189, 224, 261, 300, 341, 384, 429, 476, 525, 576, 629, 684, 741, 800, 861, 924, 989, 1056, 1125, 1196, 1269, 1344, 1421, 1500, 1581, 1664, 1749, 1836, 1925, 2016, 2109, 2204, 2301, 2400, 2501, 2604, 2709, 2816, 2925, 3036, 3149, 3264

Offset: 0

Wolfdieter Lang, Jul 20 2006

Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

a(n-10), n >= 11, tenth column (used for the n=10 series of the hydrogen atom) of triangle A120070.
For n*(n+18) see A098850.
Cf. A001008, A001477, A005563, A028347, A028552, A028557, A028560, A028563, A028566, A028569, A098603, A098847, A098848, A098849, A102928, A120061, A132759, A132760, A132761, A132762, A132765.

seq(sum(4*k-n, k=7..n), n=6..51); # Zerinvary Lajos, Feb 17 2008

s=0;lst={};Do[s+=n;AppendTo[lst,s],{n,21,6!,2}];lst (* Vladimir Joseph Stephan Orlovsky, Feb 26 2009 *)

a(n)=n*(n+20) \\ Charles R Greathouse IV, Jan 21 2015

a(n) = (n+10)^2 - 10^2 = n*(n+20), n >= 0.
G.f.: x*(21-19*x)/(1-x)^3.
a(n) = 2*n + a(n-1) + 19 (with a(0)=0). - Vincenzo Librandi, Nov 13 2010
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(20)/20 = A001008(20)/A102928(20) = 11167027/62078016, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 155685007/4655851200. (End)
From Elmo R. Oliveira, Jan 12 2025: (Start)
E.g.f.: exp(x)*x*(21 + x).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n >= 3. (End)

A132765 a(n) = n*(n + 23).

Original entry on oeis.org

0, 24, 50, 78, 108, 140, 174, 210, 248, 288, 330, 374, 420, 468, 518, 570, 624, 680, 738, 798, 860, 924, 990, 1058, 1128, 1200, 1274, 1350, 1428, 1508, 1590, 1674, 1760, 1848, 1938, 2030, 2124, 2220, 2318, 2418, 2520, 2624, 2730, 2838, 2948, 3060, 3174, 3290, 3408

Offset: 0

Omar E. Pol, Aug 28 2007

G. C. Greubel, Table of n, a(n) for n = 0..5000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, ResearchGate, 2014.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A001477, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566, A028569, A056126, A098603, A098847, A098848, A098849, A098850, A102928, A120071, A132759, A132760, A132761, A132762, A132763, A132764, A132765, A132766.

Table[n (n + 23), {n, 0, 50}] (* Bruno Berselli, Sep 03 2018 *)

a(n)=n*(n+23) \\ Charles R Greathouse IV, Jun 17 2017

[n*(n+23) for n in (0..50)] # G. C. Greubel, Mar 14 2022

a(n) = n*(n + 23).
a(n) = 2*n + a(n-1) + 22 for n>0, a(0)=0. - Vincenzo Librandi, Aug 03 2010
From Chai Wah Wu, Dec 17 2016: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2.
G.f.: 2*x*(12 - 11*x)/(1-x)^3. (End)
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(23)/23 = A001008(23)/A102928(23) = 444316699/2736605872, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/23 - 3825136961/123147264240. (End)
E.g.f.: x*(24 + x)*exp(x). - G. C. Greubel, Mar 14 2022

A132762 a(n) = n*(n + 19).

Original entry on oeis.org

0, 20, 42, 66, 92, 120, 150, 182, 216, 252, 290, 330, 372, 416, 462, 510, 560, 612, 666, 722, 780, 840, 902, 966, 1032, 1100, 1170, 1242, 1316, 1392, 1470, 1550, 1632, 1716, 1802, 1890, 1980, 2072, 2166, 2262, 2360, 2460, 2562, 2666, 2772, 2880, 2990, 3102, 3216

Offset: 0

Omar E. Pol, Aug 28 2007

G. C. Greubel, Table of n, a(n) for n = 0..5000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566, A028569, A051942, A098603, A098847, A098848, A098849, A098850, A102928, A120071, A132759, A132760, A132761, A132763, A132764, A132765, A132766, A132767.

Table[n (n + 19), {n, 0, 50}] (* Bruno Berselli, Sep 05 2018 *)
LinearRecurrence[{3,-3,1},{0,20,42},60] (* Harvey P. Dale, Jun 03 2021 *)

a(n)=n*(n+19) \\ Charles R Greathouse IV, Jun 17 2017

[n*(n+19) for n in (0..50)] # G. C. Greubel, Mar 14 2022

a(n) = 2*n + a(n-1) + 18 for n > 0, a(0) = 0. - Vincenzo Librandi, Aug 03 2010
From Chai Wah Wu, Dec 17 2016: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2.
G.f.: 2*x*(10 - 9*x)/(1-x)^3. (End)
a(n) = 2*A051942(n+9). - R. J. Mathar, Sep 05 2018
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(19)/19 = A001008(19)/A102928(19) = 275295799/1474352880, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/19 - 33464927/884611728. (End)
E.g.f.: x*(20 + x)*exp(x). - G. C. Greubel, Mar 14 2022

A132764 a(n) = n*(n+22).

Original entry on oeis.org

0, 23, 48, 75, 104, 135, 168, 203, 240, 279, 320, 363, 408, 455, 504, 555, 608, 663, 720, 779, 840, 903, 968, 1035, 1104, 1175, 1248, 1323, 1400, 1479, 1560, 1643, 1728, 1815, 1904, 1995, 2088, 2183, 2280, 2379, 2480, 2583, 2688, 2795, 2904, 3015, 3128, 3243, 3360

Offset: 0

Omar E. Pol, Aug 28 2007

			a(1)=2*1+0+21=23; a(2)=2*2+23+21=48; a(3)=2*3+48+21=75. - _Vincenzo Librandi_, Aug 03 2010

G. C. Greubel, Table of n, a(n) for n = 0..5000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566, A028569, A098603, A098847, A098848, A098849, A098850, A102928, A120071, A132759, A132760, A132761, A132762, A132763, A132765, A132766, A132767.

Table[n(n+22),{n,0,50}] (* or *) LinearRecurrence[{3,-3,1},{0,23,48},50] (* Harvey P. Dale, May 02 2012 *)

a(n)=n*(n+22) \\ Charles R Greathouse IV, Oct 07 2015

[n*(n+22) for n in (0..50)] # G. C. Greubel, Mar 14 2022

a(n) = n*(n + 22).
a(n) = 2*n + a(n-1) + 21 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
a(0)=0, a(1)=23, a(2)=48, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, May 02 2012
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(22)/22 = A001008(22)/A102928(22) = 19093197/113809696, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 156188887/5121436320. (End)
From G. C. Greubel, Mar 14 2022: (Start)
G.f.: x*(23 - 21*x)/(1-x)^3.
E.g.f.: x*(23 + x)*exp(x). (End)

A132763 a(n) = n*(n+21).

Original entry on oeis.org

0, 22, 46, 72, 100, 130, 162, 196, 232, 270, 310, 352, 396, 442, 490, 540, 592, 646, 702, 760, 820, 882, 946, 1012, 1080, 1150, 1222, 1296, 1372, 1450, 1530, 1612, 1696, 1782, 1870, 1960, 2052, 2146, 2242, 2340, 2440, 2542, 2646, 2752, 2860, 2970, 3082, 3196, 3312

Offset: 0

Omar E. Pol, Aug 28 2007

G. C. Greubel, Table of n, a(n) for n = 0..5000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566, A028569, A098603, A098847, A098848, A098849, A098850, A102928, A120071, A132759, A132760, A132761, A132762, A132764, A132765, A132766, A132767.

Table[n(n+21),{n,0,50}] (* or *) LinearRecurrence[{3,-3,1},{0,22,46},50] (* Harvey P. Dale, May 25 2014 *)

a(n)=n*(n+21) \\ Charles R Greathouse IV, Oct 07 2015

[n*(n+21) for n in (0..50)] # G. C. Greubel, Mar 14 2022

a(n) = n*(n + 21).
a(n) = 2*n + a(n-1) + 20 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
a(0)=0, a(1)=22, a(2)=46, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, May 25 2014
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(21)/21 = A001008(21)/A102928(21) = 18858053/108636528, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/21 - 166770367/4888643760. (End)
From Stefano Spezia, Jan 30 2021: (Start)
O.g.f.: 2*x*(11 - 10*x)/(1 - x)^3.
E.g.f.: x*(22 + x)*exp(x). (End)

A132766 a(n) = n*(n+24).

Original entry on oeis.org

0, 25, 52, 81, 112, 145, 180, 217, 256, 297, 340, 385, 432, 481, 532, 585, 640, 697, 756, 817, 880, 945, 1012, 1081, 1152, 1225, 1300, 1377, 1456, 1537, 1620, 1705, 1792, 1881, 1972, 2065, 2160, 2257, 2356, 2457, 2560, 2665, 2772, 2881, 2992, 3105, 3220, 3337

Offset: 0

Omar E. Pol, Aug 28 2007

G. C. Greubel, Table of n, a(n) for n = 0..5000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566, A028569, A098603, A098849, A098850, A098847, A098848, A102928, A120071, A132759, A132760, A132761, A132762, A132763, A132764, A132765.

Table[n (n + 24), {n, 0, 50}] (* or *) LinearRecurrence[{3, -3, 1}, {0, 25, 52}, 50] (* Harvey P. Dale, Feb 11 2016 *)

a(n)=n*(n+24) \\ Charles R Greathouse IV, Jun 17 2017

[n*(n+24) for n in (0..50)] # G. C. Greubel, Mar 14 2022

a(n) = n*(n + 24).
a(n) = 2*n + a(n-1) + 23 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
a(0)=0, a(1)=25, a(2)=52; for n>2, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, Feb 11 2016
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(24)/24 = A001008(24)/A102928(24) = 1347822955/8566766208, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 3602044091/128501493120. (End)
From G. C. Greubel, Mar 14 2022: (Start)
G.f.: 2*x*(13 - 12*x)/(1-x)^3.
E.g.f.: x*(26 + x)*exp(x). (End)

A132767 a(n) = n*(n + 25).

Original entry on oeis.org

0, 26, 54, 84, 116, 150, 186, 224, 264, 306, 350, 396, 444, 494, 546, 600, 656, 714, 774, 836, 900, 966, 1034, 1104, 1176, 1250, 1326, 1404, 1484, 1566, 1650, 1736, 1824, 1914, 2006, 2100, 2196, 2294, 2394, 2496, 2600, 2706, 2814, 2924, 3036, 3150, 3266, 3384

Offset: 0

Omar E. Pol, Aug 28 2007

a(n) is the Zagreb 1 index of the Mycielskian of the cycle graph C[n]. See p. 205 of the D. B. West reference. - Emeric Deutsch, Nov 04 2016

Douglas B. West, Introduction to Graph Theory, 2nd ed., Prentice-Hall, NJ, 2001.

G. C. Greubel, Table of n, a(n) for n = 0..5000
Felix P. Muga II, Extending the Golden Ratio and the Binet-de Moivre Formula, Preprint on ResearchGate, March 2014.
Wikipedia, Mycielskian.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A001008, A002378, A005563, A028347, A028552, A028557, A028560, A028563, A028566, A028569, A098603, A098847, A098848, A098849, A098850, A102928, A120071, A132759, A132760, A132761, A132762, A132763, A132764, A132765, A132766.

Table[n (n + 25), {n, 0, 50}] (* Bruno Berselli, Aug 22 2018 *)
LinearRecurrence[{3,-3,1},{0,26,54},60] (* Harvey P. Dale, Feb 20 2023 *)

a(n)=n*(n+25) \\ Charles R Greathouse IV, Jun 17 2017

[n*(n+25) for n in (0..50)] # G. C. Greubel, Mar 13 2022

a(n) = 2*n + a(n-1) + 24 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
a(n) = n^2 + 25*n. - Omar E. Pol, Nov 04 2016
From Chai Wah Wu, Dec 17 2016: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2.
G.f.: 2*x*(13 - 12*x)/(1-x)^3. (End)
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(25)/25 = A001008(25)/A102928(25) = 34052522467/223092870000, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/25 - 19081066231/669278610000. (End)
E.g.f.: x*(26 + x)*exp(x). - G. C. Greubel, Mar 13 2022

cp's OEIS Frontend

A097480 Positive integers n such that 2n-15 is prime.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A098605 Positive integers n such that 2n - 17 is prime.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A158405 Triangle T(n,m) = 1+2*m of odd numbers read along rows, 0<=m

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

Extensions

A120071 a(n) = n*(n+20).

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A132765 a(n) = n*(n + 23).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Sage

Formula

A132762 a(n) = n*(n + 19).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Sage

Formula

A132764 a(n) = n*(n+22).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Sage

Formula