guss_sediel.cpp

这是一个计算机专业计算方法的课程试验

#include<iostream.h>
double G_S[3][4],x[4]={0,0,0,0};
void main()
{
	int i ,j;
	for(i=0;i<3;i++)
		for(j=0;j<4;j++)
			cin>>G_S[i][j];
		for(i=0;i<3;i++)
			for(j=0;j<4;j++)
			{
			/*	if(G_S[i][j]==0)
				{
					for(int k=i+1;k<3;k++)
					{
						if(G_S[k][j]!=0)
							for(int l=0;l<4;l++)
							{x[l]=G_S[i][j];
							G_S[i][j]=G_S[k][l];
							G_S[k][l]=x[l];
							} 
					}
				}
				*/

				cout<<"G_S["<<i<<"]["<<j<<"]="<<G_S[i][j]<<"  ";
				if(j==3)cout<<'\n';
			}
		for(i=0;i<3;i++)
		{
			x[i]=G_S[i][i];
			for(j=0;j<4;j++)
				G_S[i][j]/=x[i];
		}
			for(i=0;i<3;i++)
				for(j=0;j<4;j++)
				{
					cout<<"G_S["<<i<<"]["<<j<<"]="<<G_S[i][j]<<"  ";
					if(j==3)cout<<'\n';
				}
	/********************************************/
			
			for(i=0;i<8;i++)//要求做的次数
			{
				for(int I=0;I<3;I++)
					for(j=0;j<3;j++)
						if(I!=j)
						{
							//x[I]=G_S[I][j];
							x[I]-=G_S[I][j]*x[j];
						}
						x[I]+=G_S[I][3];
						cout<<x[I]<<" ";
			}
}