I need some help to compute the Jacobian and Hessian of a function $f : \mathbb{R}^n \rightarrow \mathbb{R}^n$ which takes as input a vector $x$ of length $n > 0$. The other symbols can be assumed to be constant.
\begin{equation} f(x) = \sum_{i,j=1}^{n}\rho_{ij}\sigma_i\sigma_j x_i x_j \end{equation}
For the Jacobian do I just make $x_i=0$ for $\frac{\partial}{\partial x_i}f$? How does the Hessian differ once the first derivative of $x_i$ is already $0$?
Let $A=(a_{ij})$ be the matrix with entries $a_{ij}=\rho_{ij}\sigma_i\sigma_j.$ Then \begin{equation}f(x)=\langle Ax,x\rangle, \end{equation}so \begin{equation}f(x+h)=\langle A(x+h),x+h\rangle=\langle Ax,x\rangle +\langle A ,h\rangle+ \langle Ah,x\rangle+\langle Ah,h\rangle=f(x)+\langle (A+A^t)x,h\rangle+\langle Ah,h\rangle \end{equation} and therefore the Hessian matrix is just $A$ itself.