c求解非线性方程组的牛顿顿迭代法
牛顿迭代法c++程序设计求解 0=x*x-2*x-y+0.5; 0=x*x+4*y*y-4;的方程#include<iostream>#include<cmath>#define N 2
学习必备欢迎下载 牛顿迭代法c++程序设计 求解 0=x*x-2*x-y+0.5; 的方程 0=x*x+4*y*y-4; #include<iostream> #include<cmath> #defineN2// 非线性方程组中方程个数、未知量个数 #defineEpsilon0.0001//1 差向量范数的上限 #defineMax100// 最大迭代次数 usingnamespacestd; constintN2=2*N; intmain() { voidff(floatxx[N],floatyy[N]);//yy[N] 计算向量函数的因变量向量 voidffjacobian(floatxx[N],floatyy[N][N]);//yy[N][N] 计算雅克比矩阵 voidinv_jacobian(floatyy[N][N],floatinv[N][N]);//inv 计算雅克比矩阵的逆矩阵 voidnewdundiedai(floatx0[N],floatinv[N][N],floaty0[N],floatx1[N]);//x0 由近似解向量 x1 计算近似解向量 floatx0[N]={2.0,0.25},y0[N],jacobian[N][N],invjacobian[N][N],x1[N],errornorm; inti,j,iter=0; //x0 如果取消对的初始化，撤销下面两行的注释 x0 符，就可以由键盘向读入初始近似解向量 for(i=0;i<N;i++) cin>>x0[i]; cout<<""<<endl; 初始近似解向量： for(i=0;i<N;i++) cout<<x0[i]<<""; cout<<endl;cout<<endl; do { iter=iter+1; cout<<""<<iter<<""<<endl;//y0 第次迭代开始计算向量函数的因变量向量 ff(x0,y0);//jacobian 计算雅克比矩阵
