Function to calculate the Jacobian matrix for a set of nonlinear algebraic equations. A simple central differencing scheme is used.

• 用来计算一套非线性代数方程的雅可比矩阵。使用一个简单微分方法。