--- a/scripts/image/hot.m	Sun Mar 24 13:22:06 2013 -0400
+++ b/scripts/image/hot.m	Mon Mar 25 14:32:49 2013 -0700
@@ -46,12 +46,27 @@
   if (n == 1)
     map = [1, 1, 1];
   elseif (n > 1)
-    x = linspace (0, 1, n)';
-    r = (x < 2/5) .* (5/2 * x) ...
-      + (x >= 2/5);
-    g = (x >= 2/5 & x < 4/5) .* (5/2 * x - 1) ...
-      + (x >= 4/5);
-    b = (x >= 4/5) .* (5 * x - 4);
+    idx = floor (3/8 * n);
+    nel = idx;
+
+    r = ones (n, 1);
+    if (nel > 0) 
+      r(1:idx, 1) = [1:nel]' / nel;
+    endif
+
+    g = zeros (n, 1);
+    g(idx+1:2*idx, 1) = r(1:idx);
+    g(2*idx+1:end, 1) = 1;
+
+    idx = floor (3/4 * n);
+    if (any (mod (n, 8) == [0, 1, 3, 6]))
+      idx++;
+    endif
+    nel = n - idx + 1;
+
+    b = zeros (n, 1);
+    b(idx:end, 1) = [1:nel]' / nel;
+
     map = [r, g, b];
   else
     map = zeros (0, 3);