gaosixiaoyuan.rar_gauss_数值算法_线性方程组_高斯消元_高斯消元法

#include<iostream.h> #include<math.h> void main () { int i,j,n=3,k; float array[3][3]={1,2,1,2,2,3,-1,-3,8},b[3]={0,3,10},x[3],l[3][3],s=0; // cout<< "输入方程组中每一个方程的系数:"<<endl; // for(i=0;i<n;i++) // for(j=0;j<n;j++) // cin>>array[i][j]; for(i=0;i<n;i++) { for(j=0;j<n;j++) cout<<array[i][j]<<" "; cout<<endl; } // cout<<"输入b[i]:"<<endl; // for(i=0;i<n;i++) // cin>>b[i]; for(i=0;i<n;i++) cout<<b[i]<<endl; for(k=0;k<n-1;k++) for(i=k+1;i<n;i++) { if(array[k][k]!=0) { l[i][k]=array[i][k]/array[k][k]; b[i]-=l[i][k]*b[k]; } for(j=k+1;j<n;j++) array[i][j]-=l[i][k]*array[k][j]; } x[n-1]=b[n-1]/array[n-1][n-1]; cout<<"x["<<n-1<<"]="<<x[n-1]<<endl; for(k=n-2;k>=0;k--) { { for(j=n-1;j>=k+1;j--) s+=array[k][j]*x[j]; } x[k]=b[k]-s; x[k]/=array[k][k]; s=0; cout<<"x["<<k<<"]="<<x[k]<<endl; } }