diff --git a/scripts/gradients.m b/scripts/gradients.m
new file mode 100644
index 0000000000000000000000000000000000000000..04ba93b441f392c0acf0888d8a9bdb55fff3a08a
--- /dev/null
+++ b/scripts/gradients.m
@@ -0,0 +1,16 @@
+clear all
+clf
+
+tx = ty = linspace (-2, 2, 41)';
+[xx, yy] = meshgrid (tx, ty);
+tz = xx.^2 + yy.^2;
+
+mesh(tx, ty, tz)
+
+hold on
+
+[x, y] = meshgrid (-2:0.2:2);
+h = quiver (x, y, 2*x, 2*y);
+set (h, "maxheadsize", 0.33);
+
+hold off