(交换第k行和第r行的所有元素)Р for(i=0;i<n+1;i++)Р {Р temp=c[k][i];Р c[k][i]=c[r][i];Р c[r][i]=temp;Р }Р Р //高斯消元Р for(i=k+1;i<n;i++)Р {Р l=c[i][k]/c[k][k];Р for(j=k;j<n+1;j++)Р c[i][j]=c[i][j]-l*c[k][j];Р Р }Р?Р cout<<"输出增广矩阵:"<<endl;Р?Р for(i=0;i<n;i++)Р {Р for(j=0;j<n+1;j++)Р cout<<c[i][j]<<" ";Р cout<<endl;Р }Р?Р?}Р?Р?//输出上三角矩阵Р?cout<<"输出上三角矩阵:"<<endl;Р?for(i=0;i<n;i++)Р?{Р for(j=0;j<n+1;j++)Р {Р cout<<c[i][j]<<" ";Р }Р cout<<endl;Р?}Р?//回代过程Р?x[n-1]=c[n-1][n]/c[n-1][n-1];Р?for(k=n-2;k>=0;k--)Р?{ s=0;Р for(j=n-1;j>k;j--)Р s=s+c[k][j]*x[j];Р Р x[k]=(c[k][n]-s)/c[k][k];Р?}Р Р //输出方程组的解Р?cout<<"方程组的解为:"<<endl;Р for(i=0;i<n;i++)Р cout<<x[i]<<" ";Р?return 0;}Р?Р3 牛顿法解非线性方程组Р3.1牛顿法解非线性方程组算法说明Р 设已知。Р 第1步:计算函数Р Р 第2步:计算雅可比矩阵Р Р 第3步:求线性方程组Р 的解。Р 第4步:计算下一点Р Р 重复上述过程。Р3.2 牛顿法解非线性方程组算法流程图Р图3-1 算法流程图Р图3-1 算法流程图
全文预览
