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
	//將數(shù)組A分解成數(shù)組u和數(shù)組l
	//r代表的是行數(shù)i代表列數(shù)
	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求得數(shù)組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求得數(shù)組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個文件信息