Advertisement
Guest User

partial sums of iterated sine function

a guest
Jan 12th, 2013
360
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 5.01 KB | None | 0 0
  1. (*program start*)
  2. Clear[x, y, n, a, xmin, xmax, ymin, ymax, c]
  3. xmin = -1*Pi;
  4. xmax = +3*Pi;
  5. x = N[Range[xmin, xmax, 1/(50*Pi)], 90];
  6. y1 = Sin[x];
  7. y2 = Sin[y1];
  8. y3 = Sin[y1 + y2];
  9. y4 = Sin[y1 + y2 + y3];
  10. y5 = Sin[y1 + y2 + y3 + y4];
  11. y6 = Sin[y1 + y2 + y3 + y4 + y5];
  12. y7 = Sin[y1 + y2 + y3 + y4 + y5 + y6];
  13. y8 = Sin[y1 + y2 + y3 + y4 + y5 + y6 + y7];
  14. y9 = Sin[y1 + y2 + y3 + y4 + y5 + y6 + y7 + y8];
  15. y10 = Sin[y1 + y2 + y3 + y4 + y5 + y6 + y7 + y8 + y9];
  16. y11 = Sin[y1 + y2 + y3 + y4 + y5 + y6 + y7 + y8 + y9 + y10];
  17. y12 = Sin[y1 + y2 + y3 + y4 + y5 + y6 + y7 + y8 + y9 + y10 + y11];
  18. y = y1 + y2 + y3 + y4 + y5 + y6 + y7 + y8 + y9 + y10 + y11 + y12;
  19. Max[y]
  20. a = Table[{x[[n]], y[[n]]}, {n, 1, Length[x]}];
  21. a1 = Table[{x[[n]], y1[[n]]}, {n, 1, Length[x]}];
  22. a2 = Table[{x[[n]], y2[[n]] + y1[[n]]}, {n, 1, Length[x]}];
  23. a3 = Table[{x[[n]], y3[[n]] + y2[[n]] + y1[[n]]}, {n, 1, Length[x]}];
  24. a4 = Table[{x[[n]], y4[[n]] + y3[[n]] + y2[[n]] + y1[[n]]}, {n, 1,
  25. Length[x]}];
  26. a5 = Table[{x[[n]],
  27. y5[[n]] + y4[[n]] + y3[[n]] + y2[[n]] + y1[[n]]}, {n, 1,
  28. Length[x]}];
  29. a6 = Table[{x[[n]],
  30. y6[[n]] + y5[[n]] + y4[[n]] + y3[[n]] + y2[[n]] + y1[[n]]}, {n, 1,
  31. Length[x]}];
  32. a7 = Table[{x[[n]],
  33. y7[[n]] + y6[[n]] + y5[[n]] + y4[[n]] + y3[[n]] + y2[[n]] +
  34. y1[[n]]}, {n, 1, Length[x]}];
  35. a8 = Table[{x[[n]],
  36. y8[[n]] + y7[[n]] + y6[[n]] + y5[[n]] + y4[[n]] + y3[[n]] +
  37. y2[[n]] + y1[[n]]}, {n, 1, Length[x]}];
  38. a9 = Table[{x[[n]],
  39. y9[[n]] + y8[[n]] + y7[[n]] + y6[[n]] + y5[[n]] + y4[[n]] +
  40. y3[[n]] + y2[[n]] + y1[[n]]}, {n, 1, Length[x]}];
  41. a10 = Table[{x[[n]],
  42. y10[[n]] + y9[[n]] + y8[[n]] + y7[[n]] + y6[[n]] + y5[[n]] +
  43. y4[[n]] + y3[[n]] + y2[[n]] + y1[[n]]}, {n, 1, Length[x]}];
  44. a11 = Table[{x[[n]],
  45. y11[[n]] + y10[[n]] + y9[[n]] + y8[[n]] + y7[[n]] + y6[[n]] +
  46. y5[[n]] + y4[[n]] + y3[[n]] + y2[[n]] + y1[[n]]}, {n, 1,
  47. Length[x]}];
  48. a12 = Table[{x[[n]],
  49. y12[[n]] + y11[[n]] + y10[[n]] + y9[[n]] + y8[[n]] + y7[[n]] +
  50. y6[[n]] + y5[[n]] + y4[[n]] + y3[[n]] + y2[[n]] + y1[[n]]}, {n,
  51. 1, Length[x]}];
  52. ymin = Min[y]/Pi;
  53. ymax = Max[y]/Pi;
  54. ymin = -10/2*Pi;
  55. ymax = 10/2*Pi;
  56.  
  57. g = ListLinePlot[a,
  58. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  59. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  60. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2},
  61. PlotStyle -> RGBColor[1, 0, 1]];
  62. g1 = ListLinePlot[a1,
  63. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  64. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  65. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  66. g2 = ListLinePlot[a2,
  67. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  68. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  69. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  70. g3 = ListLinePlot[a3,
  71. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  72. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  73. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  74. g4 = ListLinePlot[a4,
  75. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  76. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  77. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  78. g5 = ListLinePlot[a5,
  79. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  80. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  81. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  82. g6 = ListLinePlot[a6,
  83. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  84. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  85. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  86. g7 = ListLinePlot[a7,
  87. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  88. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  89. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  90.  
  91. g8 = ListLinePlot[a8,
  92. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  93. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  94. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  95.  
  96. g9 = ListLinePlot[a9,
  97. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  98. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  99. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  100.  
  101. g10 = ListLinePlot[a10,
  102. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  103. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  104. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  105.  
  106. g11 = ListLinePlot[a11,
  107. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  108. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  109. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  110.  
  111. g12 = ListLinePlot[a12,
  112. Ticks -> {Flatten[{xmin, xmax, Table[n, {n, xmin, xmax, Pi}]}],
  113. Flatten[{ymin, ymax, Table[k, {k, ymin, ymax, Pi}]}]},
  114. ImageSize -> Full, PlotRange -> {-Pi - 1/2, Pi + 1/2}];
  115.  
  116. Show[g, g1, g2, g3, g4, g5, g6, g7, g8, g9, g10, g11, g12]
  117. Show[g5]
  118. (*program end*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement