As presented at RUMA on 2022-10-12HEADv1.0master

1 files changed, 24 insertions, 0 deletions

diff --git a/figs/atoms.fig/crystal.py b/figs/atoms.fig/crystal.py
new file mode 100644
index 0000000..4b1ba3c
--- /dev/null
+++ b/figs/atoms.fig/crystal.py
@@ -0,0 +1,24 @@
+#!/usr/bin/env python3
+
+from math import *
+import random
+import sys
+
+# size of lattice
+N=5
+
+# configuration
+config=[]
+for i in range(N):
+    for j in range(N):
+        for k in range(N):
+            config.append([2*i+((j+k)%2),sqrt(3)*(j+(k%2)/3),2*sqrt(6)/3*k])
+
+
+for i in range(len(config)):
+    print(str(config[i][0])+"+cos(u)*sin(v)", end=",")
+    print(str(config[i][1])+"+sin(u)*sin(v)", end=",")
+    print(str(config[i][2])+"+cos(v)", end=" ")
+    print("with pm3d", end="")
+    if i<len(config)-1:
+        print(", \\")