m = Moebius(10)
-ts = arange(0, 1, 0.1)
+ts = np.linspace(0, 1, 20)
-thetas = arange(0, tau, 20)
+thetas = np.linspace(0, tau, 20)
-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 ] )
+c0 = curve( color=color.blue, pos = [ m.edge.point(th) for th in thetas ] )
+c1 = curve( color=color.red, 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 ] )