Advertisement
MatsGranvik

Generalization of rule 30

Dec 7th, 2019
287
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 3.71 KB | None | 0 0
  1. (*triangle formatting*)
  2. (*start*)
  3. Clear[t, n, k]; t[n_, 1] = 0; t[n_, 2] = 0;
  4. t[n_, 3] = 1;
  5. t[n_, k_] :=
  6. t[n, k] =
  7. If[2*n + 1 >=
  8. k, ((1 - (t[n - 1, k - 2]*t[n - 1, k - 1]*t[n - 1, k]))*(1 -
  9. t[n - 1, k - 2]*
  10. t[n - 1,
  11. k - 1]*(1 - t[n - 1, k]))*(1 - (t[n - 1,
  12. k - 2]*(1 - t[n - 1, k - 1])*
  13. t[n - 1,
  14. k]))*(1 - ((1 - t[n - 1, k - 2])*(1 - t[n - 1, k - 1])*(1 -
  15. t[n - 1, k])))) + ((t[n - 1,
  16. k - 2]*(1 - t[n - 1, k - 1])*(1 - t[n - 1, k]))*((1 -
  17. t[n - 1, k - 2])*t[n - 1, k - 1]*
  18. t[n - 1, k])*((1 - t[n - 1, k - 2])*
  19. t[n - 1,
  20. k - 1]*(1 - t[n - 1, k]))*((1 - t[n - 1, k - 2])*(1 -
  21. t[n - 1, k - 1])*t[n - 1, k])), 0]; TableForm[
  22. ArrayPlot[Table[Table[t[n, k], {k, 3, 2*n + 1}], {n, 1,
  23. 100}]]] (*Mats Granvik, Dec 06 2019*)
  24. (*end*)
  25.  
  26. (*pyramid formatting*)
  27. (*start*)Clear[t, n, k];
  28. nn = 100;
  29. t[1, k_] := t[1, k] = If[k == nn, 1, 0];
  30. t[n_, k_] :=
  31. t[n, k] = ((1 - (t[n - 1, k - 1]*t[n - 1, k - 0]*
  32. t[n - 1, k + 1]))*(1 -
  33. t[n - 1, k - 1]*
  34. t[n - 1,
  35. k - 0]*(1 - t[n - 1, k + 1]))*(1 - (t[n - 1,
  36. k - 1]*(1 - t[n - 1, k - 0])*
  37. t[n - 1,
  38. k + 1]))*(1 - ((1 - t[n - 1, k - 1])*(1 -
  39. t[n - 1, k - 0])*(1 - t[n - 1, k + 1])))) + ((t[n - 1,
  40. k - 1]*(1 - t[n - 1, k - 0])*(1 - t[n - 1, k + 1]))*((1 -
  41. t[n - 1, k - 1])*t[n - 1, k - 0]*
  42. t[n - 1, k + 1])*((1 - t[n - 1, k - 1])*
  43. t[n - 1,
  44. k - 0]*(1 - t[n - 1, k + 1]))*((1 - t[n - 1, k - 1])*(1 -
  45. t[n - 1, k - 0])*t[n - 1, k + 1]));
  46. ArrayPlot[Table[Table[t[n, k], {k, 1, 2*nn}], {n, 1, nn}]]
  47. (*Mats Granvik,Dec 07 2019*)
  48. (*end*)
  49.  
  50. (*entries as polynomials*)
  51. (*start*)
  52. Clear[t, n, k];
  53. nn = 4;
  54. t[1, k_] := t[1, k] = If[k == nn, x, 0];
  55. t[n_, k_] :=
  56. t[n, k] = ((1 - (t[n - 1, k - 1]*t[n - 1, k - 0]*
  57. t[n - 1, k + 1]))*(1 -
  58. t[n - 1, k - 1]*
  59. t[n - 1,
  60. k - 0]*(1 - t[n - 1, k + 1]))*(1 - (t[n - 1,
  61. k - 1]*(1 - t[n - 1, k - 0])*
  62. t[n - 1,
  63. k + 1]))*(1 - ((1 - t[n - 1, k - 1])*(1 -
  64. t[n - 1, k - 0])*(1 - t[n - 1, k + 1])))) + ((t[n - 1,
  65. k - 1]*(1 - t[n - 1, k - 0])*(1 - t[n - 1, k + 1]))*((1 -
  66. t[n - 1, k - 1])*t[n - 1, k - 0]*
  67. t[n - 1, k + 1])*((1 - t[n - 1, k - 1])*
  68. t[n - 1,
  69. k - 0]*(1 - t[n - 1, k + 1]))*((1 - t[n - 1, k - 1])*(1 -
  70. t[n - 1, k - 0])*t[n - 1, k + 1]));
  71. TableForm[Table[Table[Expand[t[n, k]], {k, 1, 2*nn}], {n, 1, nn}]]
  72. (*Mats Granvik,Dec 07 2019*)
  73. (*end*)
  74.  
  75.  
  76. (*start*)
  77. Clear[t, n, k, x];
  78. nn = 9;
  79. t[1, k_] := t[1, k] = If[k == nn, x, 0];
  80. t[n_, k_] :=
  81. t[n, k] =
  82. t[-1 + n, -1 + k] +
  83. t[-1 + n, 0 + k]*y + (1 + t[-1 + n, 0 + k]) t[-1 + n, +1 + k];
  84. TableForm[Table[Expand[t[n, nn]], {n, 1, nn}]];
  85. TableForm[Table[CoefficientList[Expand[t[n, nn]]/x, x], {n, 1, nn}]]
  86. Table[Sum[
  87. Binomial[n, k]*Binomial[n - k, k]*Binomial[n - 2*k, 2*k], {k, 0,
  88. n}], {n, 0, 12}]
  89. (*end*)
  90.  
  91. From Wikipedia 1.1.2020:
  92.  
  93. (*start*)
  94. Clear[t, n, k, x];
  95. nn = 200;
  96. t[1, k_] := t[1, k] = If[k == nn, 0, 1];
  97. t[n_, k_] :=
  98. t[n, k] =
  99. Mod[t[-1 + n, -1 + k] + 1 + t[-1 + n, 0 + k] t[-1 + n, +1 + k], 2];
  100. ArrayPlot[Table[Table[1 - t[n, k], {k, 1, 2*nn}], {n, 1, nn}]]
  101. (*end*)
  102.  
  103. Binomials found in terms of:
  104. (*start*)
  105. (*Some of the binomial coefficients in front of terms*)
  106. Clear[t, n, k, x];
  107. nn = 16;
  108.  
  109. t[1, k_] := t[1, k] = If[k == nn, x, 0];
  110. t[n_, k_] :=
  111. t[n, k] =
  112. t[-1 + n, -1 + k] +
  113. t[-1 + n, 0 + k] + (1 + t[-1 + n, 0 + k]) t[-1 + n, +1 + k];
  114.  
  115. Do[
  116. Print[TableForm[Table[Part[t[n, nn], i], {n, 1, nn}]]], {i, 1, nn}]
  117. (*end*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement