from moebius import *
-n_u = 20
+n_u = 40
m = Moebius(n_u)
-ts = np.linspace(0, 1, 20)
+ts = np.linspace(0, 1, 40)
-thetas = np.linspace(0, tau, 20)
+thetas = np.linspace(0, tau, 40)
c0 = curve( color=color.red, pos = [ m.edge.point(th) for th in thetas ] )
c1 = curve( color=color.blue, pos = [ m.midline.point(th) for th in thetas ] )
for ix_u in range(0, n_u):
c2 = curve( pos = [ m.point(ix_u, t) for t in ts ] )
+
+for t in ts:
+ c3 = curve( color=color.yellow, pos = [ m.point(ix_u, t)
+ for ix_u in range(0, n_u+1) ] )