A Theorem Regarding Triangles
Start with any triangle at all. You can construct an equilateral triangle on each of its three sides. These three equilateral triangles have centroids, which you can connect to form another triangle. Surprisingly, this last triangle is always equilateral, regardless of the starting triangle!
Here’s a video my friend Dan over at http://danielwalsh.tumblr.com/ made with Mathematica - it demonstrates (with wacky sound effects for good measure) that no matter how you move the vertices of the starting triangle, the resultant triangle of this construction is always equilateral.
I came across this pretty geometric result by modifying a construction found in Tristan Needham’s outstanding book Visual Complex Analysis. In that book, one learns that if you construct squares on the four sides of any quadrilateral, and then connect opposite square centers, the two resulting segments are congruent and perpendicular. I had to see what happens when you do the analogous construction with triangles, and sure enough the result is as nice as anyone could hope for. Let’s now prove this fact.
Taking a cue from Visual Complex Analysis, we represent the arbitrary starting triangle using some nonzero complex numbers z1, z2 and z3. One vertex will be z1, the second will be z1 + z2, and the third will be z1 + z2 + z3, and without loss of generality, the third vertex can be the origin of the complex plane, so that z1 + z2 + z3 = 0 is assumed to hold. Now to represent the centroid of the first equilateral triangle, we go halfway to vertex z1, turn left 90 degrees, and go a distance of Sqrt[3]/6. This is because the height of an equilateral triangle is Sqrt[3]/2 times the length of one of its sides, and the centroid is one-third up the height from any side. The left turn in the complex plane is achieved via multiplication by i, which rotates the plane 90 degrees counterclockwise. Therefore the first centroid is given by (z1)/2 + i(z1)Sqrt[3]/6. See the picture:

The other two centroids are obtained from similar operations involving z2 and z3. Here I define the centroids in Mathematica:

Actually we could have turned right instead of left when finding the centroids of the three equilateral triangles of the construction - it makes no difference if these triangles are constructed ‘outward’ or ‘inward’ relative to the starting triangle. This is reflected by the fact that i and -i are algebraic mirror images of each other.
Now the first side of our centroid-triangle will be the difference (as vectors in the complex plane) of c2 and c1; similarly for the other two sides. In Mathematica I define:

We need to let Mathematica know that z1 + z2 + z3 = 0 holds for all the expressions input above - otherwise our triangle may actually be a degenerate line segment. This can be accomplished with the slash-dot & right-arrow command, which replaces certain strings with other strings. In this case we use /. z3 -> -z2 -z1 because z3 = -z2 - z1 directly follows from z1 + z2 + z3 = 0.

Now, a slick way to show that the triangle with sides s1, s2 and s3 is equilateral is to show that any of those sides can be obtained by performing a rigid-body motion on one of the others. Witness:

What this says is that s1 times (-1/2 + i Sqrt[3]/2) equals s3, and also that s2 & s1 and s3 & s2 stand in this same relation. What’s the significance of (-1/2 + i Sqrt[3]/2)? It’s the complex number that rotates the plane 120 degrees counterclockwise - in particular, multiplication by (-1/2 + i Sqrt[3]/2) is a rigid-body motion. So, this settles the theorem! But here’s another, more direct way to see it: by taking norms.

So s1, s2 and s3 all have the same absolute value or length. Without Mathematica, it’s much easier to work out the rotation approach on paper than crunching through the absolute values. But I don’t think any conceptual understanding is lost in employing the software for the grunt work.
I wonder what we can say about constructing regular pentagons on the five sides of an arbitrary pentagon? Or regular n-gons on the n sides of an arbitrary n-gon? Or what happens in higher dimensions?
And isn’t it awesome that imaginary numbers help us prove this purely geometric theorem about triangles?