Solution (source code)

= Solution

The second <partial derivatives> are $f_{xx}=4$, $f_{xy}=f_{yx}=3+2y$ and $f_{yy}=2+2x+30y$. Evaluating at the origin gives the <Hessian matrix>
$$
\boxed{H_f(0,0)=\begin{pmatrix}4&3\\3&2\end{pmatrix}.}
$$