Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- n = 100;
- \[Theta]0 = RandomReal[{0, 2 Pi}, n];
- \[Theta]p = RandomReal[{0, 0.6}, n];
- T = 20 Pi;
- \[Lambda] = 0.3;
- soln = NDSolve[
- Table[{\[Theta][i]'[
- t] == \[Theta]p[[i]] - \[Lambda] Sin[(\[Theta][i][t] -
- Sum[\[Theta][i][t], {i, 1, n}]/n)], \[Theta][i][
- 0] == \[Theta]0[[i]]}, {i, n}],
- Table[\[Theta][i], {i, n}], {t, 0, T}];
- Manipulate[Graphics[
- {
- Block[{\[Mu]},
- \[Mu] =
- If[norm, First[Sum[\[Theta][i][t], {i, 1, n}]/n /. soln], 0];
- Table[
- {
- Opacity[0.5],
- ColorData["DarkRainbow"][i/n],
- Disk[{Cos[\[Theta][i][t] - \[Mu]], Sin[\[Theta][i][t] - \[Mu]]},
- 0.1] /. soln}, {i, n}]
- ]
- }
- , PlotRange -> 1.2]
- ,
- {t, 0, T},
- {norm, {True, False}}]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement