高斯消去法求方程组.cpp

利用高斯消去法求线性方程组 自己写的 简陋了点

#include<iostream.h>
#include<math.h>
 void main()
 { int n,i,k,z,y,m;
   double a[100][100],l[100][100],x[100],b[100],c[100];
   cout<<"n="<<endl;
   cin>>n;
   for(i=1;i<=n;i++)
   {for (k=1;k<=n;k++)
   {cout<<"a"<<i<<k<<"=";
	   cin>>z;
   a[i][k]=z;
   }
   }
   for(i=1;i<=n;i++)
   {cout<<"b"<<i<<"=";
   cin>>y;
   b[i]=y;
   }
  for(i=1;i<=n;i++)
   {for(k=1;k<=n+1;k++)
   {c[i]+=a[i][k]*x[k];
   }
   }
 for(k=1;k<=n;k++)
   {for(i=n;i>k;i--)
   {l[i][k]=a[i][k]/a[k][k];
   for( m=1;m<=n;m++)
	   a[i][m]=a[i][m]-l[i][k]*a[k][m];
       b[i]=b[i]-l[i][k]*b[k];
	   cout<<"l"<<i<<k<<"="<<l[i][k]<<"\n";
	   cout<<"b"<<i<<"="<<b[i]<<"\n";
   }

 }
 if(k=n)
	 x[n]=b[n]/a[n][n];
     cout<<"x"<<n<<"="<<x[n]<<"\n";
   for(i=n-1;i>0;i--)
   {for(m=n;m>i;m--)
   {   b[i]=b[i]-a[i][m]*x[m];}

	   x[i]=b[i]/a[i][i];
	   cout<<"x"<<i<<"="<<x[i]<<"\n";
   }
    gh
 }