jiangli1.cpp

自己编的高斯消元法求线性方程组的根 可以供数值算法的朋友参考

#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; 
				} 
			
}