Given an $n$-dimensional variable $\mathbf{x}\in\mathbb{R}^n$ and the functions $h_i: \mathbb{R}^n \rightarrow \mathbb{R}$, $i=1,\dots,l$, we would like to find a solution of the following equation:
$$ \mathbf{x} + \sum_{i=1}^{l} exp\{-h_i^2(\mathbf{x})\} \frac{\partial}{\partial\mathbf{x}}h_i(\mathbf{x}) + \mathbf{c} = \mathbf{0}, $$
where
$$ h_i(\mathbf{x}) = -\frac{\mathbf{x}^T\mathbf{a}_i+b-1}{\sqrt{2\mathbf{x}^T\Sigma_i\mathbf{x}}}, $$
$\mathbf{c}=(c_1,\cdots,c_n)^T, \mathbf{a}_i=(a_{i,1},\cdots,a_{i,n})^T \in \mathbb{R}^n$, $b\in\mathbb{R}$, and $\Sigma_i$ is an $n\times n$ positive definite symmetric matrix, i.e., $\Sigma_i \in \mathbb{S}_{++}^{n}$.
Thus, unless I have some mistake, it holds that
$$ \frac{\partial}{\partial\mathbf{x}}h_i(\mathbf{x}) = \frac{(\mathbf{x}^T\mathbf{a}_i+b-1)2\mathbf{x}^T\Sigma_i - (2\mathbf{x}^T\Sigma_i\mathbf{x})\mathbf{a}_i^T}{2\mathbf{x}^T\Sigma_i\mathbf{x}\sqrt{2\mathbf{x}^T\Sigma_i\mathbf{x}}} $$
Is that a feasible equation? Could I find a $\mathbf{x}\in\mathbb{R}^n$ that satisfies the above equation?
Thanks in advance!