plots

abstable.py

@@ -5,8 +5,8 @@ from idmlin import IDMLin
 
 
 t = 1.5
-amin, amax, acount = 0.2, 0.6, 15
-bmin, bmax, bcount = 2, 4, 20
+amin, amax, acount = 0.1, 0.6, 50
+bmin, bmax, bcount = 2, 4, 100
 
 arng = linspace(amin, amax, acount)
 brng = linspace(bmin, bmax, bcount)
@@ -14,10 +14,10 @@ s0 = 2
 v0 = 30
 re = 0.5
 
-ns = (12, 15, 20, 50)
+ns = (10, 25, 50)
 
-#f, axarr = plt.subplots(2, len(ns)+1, sharex=True, sharey=True)
-f, axarr = plt.subplots(2, 2, sharex=True, sharey=True)
+f, axarr = plt.subplots(2, len(ns)+1, sharex=True, sharey=True)
+#f, axarr = plt.subplots(2, 2, sharex=True, sharey=True)
 
 allstab = zeros((len(arng), len(brng)))
 alloss = zeros((len(arng), len(brng)))
@@ -39,12 +39,12 @@ for i, n in enumerate(ns):
             stab[ai, bi], oss[ai, bi] = isstable(a,b)
 
 
-    axarr[0, i*0].imshow(stab, cmap=plt.cm.Reds, 
+    axarr[0, i*1].imshow(stab, cmap=plt.cm.Reds, 
                     interpolation='none', 
                     extent=[bmin, bmax, amax, amin], 
                     aspect='auto',
                     )
-    axarr[1, i*0].imshow(oss, cmap=plt.cm.Reds, 
+    axarr[1, i*1].imshow(oss, cmap=plt.cm.Reds, 
                     interpolation='none', 
                     extent=[bmin, bmax, amax, amin], 
                     aspect='auto',