How to determine that a surface is symmetric Given a surface $f(x,y,z)=0$, how could you determine that it's symmetric about some plane, and, if so, how would you find this plane.
The special case where $f$ is a polynomial is of some interest.
The question is somewhat related to this one: How to identify surfaces of revolution
Here's a simple 2D example to experiment with:
$$
27 x^3 + 108 x^2 y + 144 x y^2 + 64 y^3
         - 80 x^2 + 120 x y  -   45 y^2 - 200 x   + 150 y - 125 = 0$$
This one is symmetric, as the following picture suggests:

In fact, if we do the translation/rotation described by the substitution $u = \tfrac15(3x+4y)$, $v = \tfrac15(4x−3y) + 1$, then the curve is just $u^3 - v^2 = 0$, which is obviously symmetric about the line $v=0$.
But how would you discover this translation/rotation if I didn't tell you, and how would you do similar things in the 3D surface case?
Added after a few days of thought:
We can consider the surface as an object made of thin sheet metal. As such, it has a center of mass, provided it's bounded, and any plane(s) of symmetry must pass though this center of mass. The plane then has only two remained degrees of freedom, so may be easier to find.
The same sort of reasoning applies in the case of a bounded 2D curve. Again, any line of symmetry must pass through the curve's centroid, so it has only one remaining degree of freedom, namely its slope/angle.
For curves and surfaces given by implicit equations, I don't really know how to calculate centroids, but I expect this can be done.
Fabricated from Comments Below:
Several people suggested looking at highest-degree terms only. So, in my example, we just look at the equation
$$
27 x^3 + 108 x^2 y + 144 x y^2 + 64 y^3 = 0
$$
Putting $w = y/x$, this is roughly equivalent to
$$
64w^3 + 144w^2 + 108 w + 27 = 0
$$
But the polynomial on the left is just $(4w + 3)^3$, so we have a root $w = -3/4$, with multiplicity three. Is the repeated root an accident that happens only in this case, or will it always happen?? Anyway, the vector $(-4,3)$ gives us the normal to the line of symmetry, and that surely can not be an accident.
I don't really understand why this magic process works, but it looks very promising for the 2D curve case.
I don't know how to generalize to the 3D surface case, or to non-polynomial cases.
 A: Here is an answer to this question in the case of compact surfaces (without  boundary); maybe these ideas can be used in the general case as well. Let $F$ be a compact surface in $R^3$, bounding a solid $S$. In the setting you are interested in, one will probably have $F=\{x: f(x)=c\}$ and $S=\{x: f(x)\le c\}$. I will also assume that $f$ is a polynomial (I do not think it is really necessary, but I use this assumption to detect surfaces of revolution). 
The solid $S$ then has a Loewner-Jones ellipsoid, which is the unique ellipsoid $E$ of the least volume containing $F$ (in fact, every compact in $R^n$ with not contained in a hyperplane has such ellipsoid). This ellipsoid is also known as Jones ellipsoid and least volume ellipsoid. It was discovered by Loewner and then rediscovered by Jones, if I remember correctly. People also use the largest volume ellipsoid contained in $S$, that one will also work. 
The key is that, in view of uniqueness of $E$, every symmetry of $F$ is also a symmetry of $E$. There is a vast literature in the computational math community describing various algorithms for finding $E$, here is just a random paper on this subject which I found by googling. Now, suppose you are lucky and your ellipsoid has distinct (beyond the margin of error in the algorithm you would be using) principal axes. Then $E$ has exactly 3 planes of symmetry which pass through its center $C$ and pairs of axes. Then you test if symmetries in these planes preserve $S$. 
Now, suppose you are unlucky and $E$ has a rotational symmetry in its axis $A$ (but is not the sphere). By intersecting $S$ with the plane $P$ orthogonal to $A$, we now get a 2-dimensional problem: Given a (closed and bounded) curve $\Gamma$ in the plane and a center $C$, determine if $\Gamma$ has a reflection symmetry in a line $L$ (through $C$). This is relatively easy since each point of intersection of $\Gamma$ with $L$ is a critical point of the distance function from $C$, so you can find such points using Lagrange multipliers and then check if any of the corresponding symmetries preserve $F$. This method will fail if $\Gamma$ is a circle. Instead of intersecting $F$ with planes passing through $C$, we can use other planes $P_k$ orthogonal to $L$. Now, if too many  (comparing to the degree of $f$) of these intersections with $F$ are circles, then $F$ is a surface of revolution with the axis $A$. (I think, it should follow from Bezout theorem, but this has to be checked.) In this case, of course, any plane through $A$ will work.    
The remaining case when $E$ is a sphere is similar to the circular case in one dimension lower, you would be then looking for critical points of the distance function from the origin. 
Remark: The argument with Lowener-Jones ellipsoid should also handle your second question, about surfaces of revolution, in the bounded case. 
