What are the extrema of $f(x, y, z) = xyz$ on the intersection of the cylinder $x^2 + y^2 = r^2$ and the plane through the origin $ax + by + cz = 0$ in $\mathbf{R^3}$? It is sensible to assume that $c \neq 0$, so that the plane and the cylinder intersect to form an ellipse. Note that this function has extrema in this case, since the intersection is compact.
I tried to use Lagrange multipliers, which leads to the system of equations
\begin{align*} yz &= 2 \lambda_1 x + \lambda_2 a \\ xz &= 2 \lambda_1 y + \lambda_2 b \\ xy &= \lambda_2 c \\ ax + by + cz &= 0 \\ x^2 + y^2 &= r^2 \end{align*} where $\lambda_1$ and $\lambda_2$ are the Lagrange multipliers. Giving these to a computer algebra system gives extraordinarily complicated results.
Is there a better way to approach this problem?
HINT
Let parametrize as
then
$$f(x,y,z)=g(\theta)=\frac{r^3}{c}\cos\theta\sin \theta (-a\cos \theta-b\sin \theta)$$