The function $x^3 - 3x^2 - 6xy + 2y^2+ 8y + 8$ yields
$f_x = 3x^2 - 6x - 6y$
$f_y = 4y + 8 - 6x$
When we equate both of these to $0, y = \frac{x^2}{2} - x$ and $y = \frac{3}{2}x - 2$
Plotting these on a graph we can see the points of interaction are $(1, -0.5)$ and $(4, 4)$ and although that is the answer, I was wondering how to get those results using simultaneous equations, as the $x^2$ means that Gaussian elimination is inadequate, but I can't seem to figure it out using simultaneous equations.
Also, is the equation of a tangent plane the same formula as linear approximation? The tangent plane equation I found when at $(1, -2)$ is 13 + 9x - 6y