高斯_塞德尔迭代求解方程组

#include "iostream.h"
void main()
{  
int i,j,n;
    float **a,e;
void gaussSeidel(float **a,float e,int n);
cout<<"请输入方程次数和要求精度:";
cin>>n>>e;
a=new float*[n];
for(i=0;i<n;i++)
  a=new float[n+2];
cout<<"请输入增广矩阵:"<<endl;
for(i=0;i<n;i++)
  for(j=0;j<n+1;j++)
   cin>>a[j];
cout<<"请输入迭代初值:"<<endl;
for(i=0;i<n;i++)
  cin>>a[n+1];
gaussSeidel(a,e,n);
cout<<"方程的解为："<<endl;
for(i=0;i<n;i++)
  cout<<"X"<<i+1<<"="<<a[n+1]<<endl; 
}
void gaussSeidel(float **a,float e,int n)
{
int i,j,flag;
float k,t;
for(i=0;i<n-1;i++) [url=file://选/]file://选[/url]主元素
{
  for(flag=i,j=i+1;j<n;j++)
   if(a[j]>a[flag])
    flag=j;
  if(flag!=i)
   for(j=0;j<n+1;j++)
   {
    t=a[flag][j];
    a[flag][j]=a[j];
    a[j]=t;
   }
}
for(i=0;i<n;i++)    [url=file://建/]file://建[/url]立迭代矩阵
  for(k=a,j=0;j<n+1;j++)
   a[j]/=k;
while(flag)   [url=file://迭/]file://迭[/url]代过程
{
  flag=0;
  for(i=0;i<n;i++)
  {
   for(k=a[n],j=0;j<n;j++)
    if(j!=i)
     k-=a[j]*a[j][n+1];
   if(k-a[n+1]>e)
    flag=1;
   a[n+1]=k;
  } 
}  
}