To solve a system of nonlinear equations which are equal to zero, ${f_i}(X) = 0\forall i \in 1:m,X = \left[ {{x_1},{x_2}, \cdots ,{x_n}} \right]$, we use Jacobian matrix and Newton algorithm. Hence, at each iteration, updated $X$ can be find by the following equation
${X^{(k + 1)}} = {X^{(k)}} + {J^{ - 1}}\left[ {\begin{array}{*{20}{c}} {{f_i}({X^{(k)}})} \\ {{f_2}({X^{(k)}})} \\ \vdots \\ {{f_m}({X^{(k)}})} \end{array}} \right]k = 0:j$
Where $J = \left[ {\begin{array}{*{20}{c}} {\frac{{\partial {f_1}({X^k})}}{{\partial {x_1}}}}&{\frac{{\partial {f_1}({X^k})}}{{\partial {x_2}}}}& \cdots &{\frac{{\partial {f_1}({X^k})}}{{\partial {x_n}}}}\\ {\frac{{\partial {f_2}({X^k})}}{{\partial {x_1}}}}&{\frac{{\partial {f_2}({X^k})}}{{\partial {x_2}}}}& \cdots &{\frac{{\partial {f_2}({X^k})}}{{\partial {x_n}}}}\\ \vdots & \vdots & \vdots & \vdots \\ {\frac{{\partial {f_m}({X^k})}}{{\partial {x_1}}}}&{\frac{{\partial {f_m}({X^k})}}{{\partial {x_2}}}}& \cdots &{\frac{{\partial {f_m}({X^k})}}{{\partial {x_n}}}} \end{array}} \right]$
The question is how to represent the above in terms of Jacobian matrix and first differential equations of ${f_i}(X)$ or $X$.