Gauss消去法、列主元素消去法、Doolittle三角分解法c++程序

Gauss消去法、列主元素消去法、Doolittle三角分解法c++程序，编写Gauss消去法类、列主元素消去法类、Doolittle三角分解法类，再通过run.cpp主程序调用，并且把求解的每一步都打印出来了。

class Doolittle: public IOput
{
public:
	double l[N][N]u[N][N]y[N];
	Doolittle（） {}
	void FenJie（）;
	void QiuJie（）;
	~Doolittle（） {}
};
void Doolittle::FenJie（）
{
	int ijkr;
	for（i=0;i	{
		x[i]=y[i]=0;
		for（j=0;j		l[i][j]=u[i][j]=0;
	}
	double de;
	for（i=0;i
	//将数组A分解成数组u和数组l
	//r代表的是行数i代表列数
	for（r=0;r		if（r==0）
		{
			for（i=0;i			for（i=1;i		}
		else 
		{
			for（i=r;i			{
				d=0;
				for（k=0;k<=r-1;k++） 
					d=d+l[r][k]*u[k][i];
				u[r][i]=A[r][i]-d;
			}//公式3.6			
			for（i=r+1;i				if（r!=m-1）
				{
					e=0;
					for（k=0;k<=r-1;k++）
						e=e+l[i][k]*u[k][r];
						l[i][r]=（A[i][r]-e）/u[r][r];
				}//公式3.7
         }
	cout<<“输出矩阵l:“<	for（i=0;i		for（j=0;j		{
			cout<			if（j==m-1）
				cout<		}
		cout<<“输出矩阵u:“<	for（i=0;i		for（j=0;j		{
			cout<			if（j==m-1）
				cout<		}
}
void Doolittle::QiuJie（）
{
	int ir;
	double fg;
	for（r=0;r<=m-1;r++）//lY=B求得数组Y
		if（r==0）
			y[r]=b[r];//公式3.10
		else 
		{
			f=0;
			for（i=0;i				f=f+l[r][i]*y[i];
		    y[r]=b[r]-f;
		}//公式3.10
  
	for（r=m-2;r>=0;r--）//uX=Y求得数组X
		if（r==m-1） x[r]=y[r]/u[m-1][m-1];//公式3.11
		else 
		{
			g=0;
			for（i=r+1;i				g=g+u[r][i]*x[i];
			x[r]=（y[r]-g）/u[r][r];
		}//公式3.11
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      21784  2010-03-04 15:12  Gauss_Doolittle\Debug\Gauss&Doolittle.obj

     文件     548939  2010-03-05 18:41  Gauss_Doolittle\Debug\Gauss.exe

     文件     792688  2010-03-05 18:41  Gauss_Doolittle\Debug\Gauss.ilk

     文件    2017284  2010-03-05 18:54  Gauss_Doolittle\Debug\Gauss.pch

     文件    1123328  2010-03-05 18:41  Gauss_Doolittle\Debug\Gauss.pdb

     文件        203  2010-03-05 14:59  Gauss_Doolittle\Debug\Gauss.sln

    ..A..H.      7680  2010-03-05 14:59  Gauss_Doolittle\Debug\Gauss.suo

     文件     557132  2010-03-04 13:23  Gauss_Doolittle\Debug\Gauss1.exe

     文件     804676  2010-03-04 13:23  Gauss_Doolittle\Debug\Gauss1.ilk

     文件     297369  2010-03-04 13:23  Gauss_Doolittle\Debug\Gauss1.obj

     文件    2272008  2010-03-04 09:19  Gauss_Doolittle\Debug\Gauss1.pch

     文件    1123328  2010-03-04 13:23  Gauss_Doolittle\Debug\Gauss1.pdb

     文件     553045  2010-03-02 21:23  Gauss_Doolittle\Debug\Gauss_Doolittle.exe

     文件     798952  2010-03-02 21:23  Gauss_Doolittle\Debug\Gauss_Doolittle.ilk

     文件    2247428  2010-03-02 21:03  Gauss_Doolittle\Debug\Gauss_Doolittle.pch

     文件    1115136  2010-03-02 21:23  Gauss_Doolittle\Debug\Gauss_Doolittle.pdb

     文件     557137  2010-03-04 09:02  Gauss_Doolittle\Debug\Gauss消去法.exe

     文件     802712  2010-03-04 09:02  Gauss_Doolittle\Debug\Gauss消去法.ilk

     文件     294464  2010-03-04 09:08  Gauss_Doolittle\Debug\Gauss消去法.obj

     文件    2272008  2010-03-03 19:17  Gauss_Doolittle\Debug\Gauss消去法.pch

     文件    1123328  2010-03-04 09:02  Gauss_Doolittle\Debug\Gauss消去法.pdb

     文件     224873  2010-03-03 19:41  Gauss_Doolittle\Debug\InputData.obj

     文件     557142  2010-03-06 20:27  Gauss_Doolittle\Debug\Run.exe

     文件     800240  2010-03-06 20:27  Gauss_Doolittle\Debug\Run.ilk

     文件     281682  2010-03-06 20:27  Gauss_Doolittle\Debug\Run.obj

     文件    2034292  2010-03-06 20:04  Gauss_Doolittle\Debug\Run.pch

     文件    1115136  2010-03-06 20:27  Gauss_Doolittle\Debug\Run.pdb

     文件        203  2010-03-06 14:43  Gauss_Doolittle\Debug\Run.sln

    ..A..H.      8192  2010-03-06 14:43  Gauss_Doolittle\Debug\Run.suo

     文件      91136  2010-03-06 20:27  Gauss_Doolittle\Debug\vc60.idb

............此处省略16个文件信息