高斯-主销元.txt

高斯主元,雅可比迭代以及牛顿算法 对于数值分析有很大的帮助

#include<iostream.h>
#include<math.h>

void main()
{float a[10][10],b[10],m[10][10],x[10],sum;
 int i,j,k,n;
cout<<"请输入行列式 a[n][n] 的 n 的值:"<<endl;
 cin>>n;
cout<<endl;
cout<<"请输入行列式矩阵 a[n][n]"<<endl;
 for(i=0;i<n;i++)
  for(j=0;j<n;j++)
   cin>>a[i][j] ;
  cout<<"请输入常数矩阵 b[i]"<<endl;
    for(i=0;i<n;i++)
     cin>>b[i];
	cout<<"这道题目的结果是"<<endl;
      for(k=0;k<n-1;k++)
      {if(a[k][k]==0)
       cout<<"error"<<endl;
	else for(i=k+1;i<n;i++)
	     {m[i][k]=a[i][k]/a[k][k];
	      a[i][k]=m[i][k];
	      b[i]=b[i]-m[i][k]*b[k];
	 for(j=k+1;j<n;j++)
	  a[i][j]=a[i][j]-m[i][k]*a[k][j];
      }}
 if(a[n-1][n-1]==0)
 cout<<"error"<<endl;
 else x[n-1]=b[n-1]/a[n-1][n-1];
 b[n-1]=x[n-1];
 for(i=n-2;i>=0;i--)
 {sum=0;
  for(j=i+1;j<n;j++)
   {sum+=a[i][j]*x[j];}
    x[i]=(b[i]-sum)/a[i][i];
    b[i]=x[i];
   }
 for(i=0;i<n;i++)
 cout<<x[i]<<endl;
}