If we don't know what to do, we can always try to use the partial derivative test for
looking for extrema of two variables function
(look also @Christian Blatter answer)
$$\frac{\partial}{\partial x}f(x,y) = 9 x^2+2 x (4 y-7)+y^2-7 y+3 = 0$$
$$\frac{\partial}{\partial y}f(x,y) = 4 x^2+2xy -7x+1 = 0$$
From the second equation we have $ y = \frac{7x-1-4x^2}{2x}$, so if we put it in the first one, we will get
$$ -3x^2 +14x - \frac{45}{4} + \frac{1}{4x^2} = 0 $$
which is equivalent to
$$ 12x^4 - 56x^3 + 45x^2 -1 = 0 $$
We already know one of the roots of this equation ($x_{1}=1$), thus by dividing polynomial above by $(x-1)$, we get
$$ 12x^3 - 44x^2 +x +1 =0$$
Fortunately, this equation have rational root, which we seek manually by testing $$ -1,1,-\frac{1}{2},\frac{1}{2},-\frac{1}{3},\frac{1}{3},-\frac{1}{4},\frac{1}{4},-\frac{1}{6},\frac{1}{6},-\frac{1}{12},\frac{1}{12}$$
Because we only looks for $x>0$, thus we ignore negative ones, and finally get the root for $ x_{2} = \frac{1}{6} $. The other two roots we get from
$$ (12x^3 - 44x^2 +x +1)/(x-\frac{1}{6}) = 6 (2 x^2-7 x-1) = 0 $$
which are $x_{3} = \frac{1}{4} (7-\sqrt{57}) < 0 $ and $x_{4} = \frac{1}{4} (7+\sqrt{57})$
Now, if we put $x_{1}, x_{2} \text{ and } x_{4}$ into $\frac{\partial}{\partial y}f(x,y) =0$, we get
$$ y_{1} = \frac{7x_{1}-1-4x_{1}^2}{2x_{1}} = 1 $$
$$ y_{2} = \frac{7x_{2}-1-4x_{2}^2}{2x_{2}} = \frac{1}{6} $$
$$ y_{3} = \frac{7x_{4}-1-4x_{4}^2}{2x_{4}} = \frac{1}{4} (7-3 \sqrt{57}) < 0 $$
So we can have extrema (for $x,y > 0$) in $(1,1)$ and $(\frac{1}{6},\frac{1}{6})$.
Now we seek for second partial derivatives:
$$ \frac{\partial^2 f(x,y)}{\partial x^2} = 18x +8y - 14 $$
$$ \frac{\partial^2 f(x,y)}{\partial x \partial y} = 8x +2y -7 $$
$$ \frac{\partial^2 f(x,y)}{\partial y^2} = 2x $$
Now, the determinant of Hessian matrix $D(x,y)$, gives us
$ D(1,1) = 15 > 0 $, $ \frac{\partial^2 f(1,1)}{\partial x^2} = 12 > 0 $, thus minimum $f(1,1) = 0$;
$ D(\frac{1}{6},\frac{1}{6}) = -\frac{285}{9} < 0 $, $ \frac{\partial^2 f(\frac{1}{6},\frac{1}{6})}{\partial x^2} = \frac{-29}{3} < 0 $, $ \frac{\partial^2 f(\frac{1}{6},\frac{1}{6})}{\partial y^2} = \frac{1}{3} > 0 $, thus saddle point,
$f(\frac{1}{6},\frac{1}{6}) = 9 - \frac{8664}{1296} = \frac{3000}{1296} > 0$
Also, for $x,y>0$ we have
$$ f(x,0) = (3x-7)(x^2+1)+9 > 0 $$
$$ f(0,y) = y+2 > 0 $$
and I think it should be enough now. (If not, there would be extrema somewhere between, where function would need to take negative value, but except $(1,1)$ and $(\frac{1}{6},\frac{1}{6})$ there are no such points).
Sorry for so long and hope there is no formal gap above (tell me if there is, then I will edit it)