Fix indexing in particle heading plots for Vicsek model validation

Author: Daniel Precioso, PhD

modules/collective-motion/vicsek-validation.qmd

@@ -37,8 +37,8 @@ for t in range(10):
 # Plot the headings of the particles
 plt.figure(figsize=(6, 6))
 plt.quiver(
-    xy[:, 0],
-    xy[:, 1],
+    xy[0],
+    xy[1],
     np.cos(theta),
     np.sin(theta),
     angles='xy',
@@ -98,8 +98,8 @@ for t in range(10):
 # Plot the headings of the particles
 plt.figure(figsize=(6, 6))
 plt.quiver(
-    xy[:, 0],
-    xy[:, 1],
+    xy[0],
+    xy[1],
     np.cos(theta),
     np.sin(theta),
     angles='xy',
@@ -158,8 +158,8 @@ for t in range(10):
 # Plot the headings of the particles
 plt.figure(figsize=(6, 6))
 plt.quiver(
-    xy[:, 0],
-    xy[:, 1],
+    xy[0],
+    xy[1],
     np.cos(theta),
     np.sin(theta),
     angles='xy',
@@ -218,8 +218,8 @@ for t in range(10):
 # Plot the headings of the particles
 plt.figure(figsize=(6, 6))
 plt.quiver(
-    xy[:, 0],
-    xy[:, 1],
+    xy[0],
+    xy[1],
     np.cos(theta),
     np.sin(theta),
     angles='xy',