高斯法解线性方程组.cpp

高斯法解线性方程组

#include"linequ.h"

void main()
{
	cout<<"                       欢迎进入解线性方程组程序";
begin:
	int n;
	cout<<endl<<endl<<endl<<"输入阶数(n>=2):";
	cin>>n;
	double *a,*b;
	a=new double[n*n];
	b=new double[n];
	for(int i=0;i<n*n;i++){
		int c,d;c=i%n;d=(i-c)/n;
		cout<<"输入a["<<d+1<<"]["<<c+1<<"]:";
		cin>>a[i];
	}
	for(int j=0;j<n;j++){
		cout<<"输入b["<<j+1<<"]:";
		cin>>b[j];
	}
	Linequ equ1(n);
	equ1.setLinequ(a,b);
	equ1.printL();
	if(equ1.Solve())
		equ1.showX();
	else
		cout<<"Fail"<<endl;
	cout<<endl<<endl<<"是否解新的方程?(y,n):";
	char answer;
	cin>>answer;
	if(answer=='y')
		goto begin;
}