Advertisement
Guest User

Untitled

a guest
Mar 23rd, 2018
1,671
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.27 KB | None | 0 0
  1. n = 10;
  2. dy = 10;
  3. ss = Table[s, {s, 0.1, 2, 0.1}];
  4. d[s_, t_,
  5. e_] := (s EllipticE[
  6. ArcTan[s Tan[2 Pi t]] + Floor[4 t] Pi - Floor[2 t] Pi, e]);
  7. frame[t_] := (Show[{
  8. Table[
  9. Module[{l, e},
  10. e = 1 - 1/s^2;
  11. l = 4 s EllipticE[e];
  12. ParametricPlot[
  13. {
  14. {d[s, t, e] \[Theta]/(2 Pi), dy s},
  15. {Cos[\[Theta]] Sin[2 \[Pi] t] +
  16. s Cos[2 \[Pi] t] Sin[\[Theta]],
  17. Cos[2 \[Pi] t] Cos[\[Theta]] -
  18. s Sin[2 \[Pi] t] Sin[\[Theta]]} + {
  19. d[s, t, e] -
  20. ((s^2 - 1) Sin[2 \[Pi] t] Sign[Cos[2 Pi t]])/Sqrt[
  21. 1 + s^2 Tan[2 \[Pi] t]^2],
  22. Abs[Cos[2 \[Pi] t]] Sqrt[1 + s^2 Tan[2 \[Pi] t]^2] + dy s}
  23. },
  24. {\[Theta], 0, 2 Pi}, Axes -> None,
  25. PlotStyle -> If[s == 1, Red]]
  26. ]
  27. , {s, ss}],
  28. Graphics[{
  29. Line[{{0, 0}, {0, 22}}],
  30. {
  31. Module[{ps},
  32. ps = Table[{d[s, t, 1 - 1/s^2], dy s}, {s, ss}];
  33. {Point[ps],
  34. Darker@Red, Thick, Line[ps]}
  35. ],
  36. Opacity[0.4],
  37. Table[{
  38. Line[{{0, dy s}, {7, dy s}}]}
  39. , {s, ss}]
  40. }}]
  41. }
  42. , PlotRange -> {{-2.1, 12}, {0, 24}}, ImageSize -> 300]
  43. );
  44. Manipulate[
  45. frame[t],
  46. {t, 0, 1}]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement