牛顿迭代法 c++程序设计求解 {0=x＊x-2＊x-y+0.5； 0=x*x+4 ＊ y*y—4;}的方程＃include＃include#define N 2 // 非线性方程组中方程个数、未知量个数 #define Epsilon 0.0001 // 差向量 1 范数的上限#define Max 100 //最大迭代次数using namespace std；const int N2=2＊N;int main（){void ff（float xx［N］,float yy［N］）； //计算向量函数的因变量向量 yy[N］void ffjacobian(float xx[N]，float yy[N］[N]）；/ /计算雅克比矩阵 yy［N]［N]void inv_jacobian(float yy[N]［N］,float inv[N]［N］）； //计算雅克比矩阵的逆矩阵 invvoid newdundiedai(float x0［N］， float inv［N］[N]，float y0[N],float x1［N])； //由近似解向量 x0 计算近似解向量 x1float x0［N］={2.0,0.25}，y0［N]，jacobian［N］［N］，invjacobian[N][N］，x1［N],errornorm；int i,j，iter=0； //假如取消对 x0 的初始化，撤销下面两行的注释符， 就可以由键盘向 x0 读入初始近似解向量for( i=0；ix0［i］;cout<〈”初始近似解向量："<〈endl;for (i=0；i
