c语言描述超松弛算法的源代码

对超松弛算法用c语言进行了详细的描述，可在vc++6.0上运行

#include 
#include
#include 
#include
using namespace std;
void function（char *file）
{

	ifstream FileIn;
	FileIn.open（file）;
	/////////////////////////////////////////////////////
	//动态给二位数组分配存储空间
	int line;
	int row;
	FileIn>>line;
	FileIn>>row;
	double **coefficient_matrix;
	coefficient_matrix = new double *[line];
	int loop;
	for（loop = 0;loop <  line;loop++）
		coefficient_matrix[loop] = new double [row];
	//由用户自行输入系数矩阵A
	cout.setf（ios_base::left ios_base::adjustfield ）;
	int factor;
	for （loop = 0;loop< line;loop++）
		for（factor = 0;factor <  row;factor++）
		{
			//cout <<“ 请输入第“<			FileIn >>   coefficient_matrix[loop][factor];
		}
		///////////////////////////////////////
		double *constant_Matrix ;
		constant_Matrix = new double [line];
		//cout << “\n请输入常系数矩阵b：“<		for（loop = 0;loop < line;loop++）
		{
			//cout <<“请输入第“<			FileIn >>constant_Matrix[loop];
		}
		/////////////////////////////////////////
		double *First_iteratorValue; 
		First_iteratorValue = new double  [line];//给行分配大小
		//////////////////////////////////////////
		//输入数据
		//cout<<“\n请输入近似解的初始矩阵值：“<		for（loop = 0;loop		{
			//cout <<“请输入第“<			FileIn>> First_iteratorValue[loop];
		}
		///////////////////////////////////////
		//声明x（k+1）
		double  *Next_iterator = new double [line];
		for（loop = 0;loop < line;loop++）
			Next_iterator[loop] = 0;
		int count = 0;
		double MaxPrecision = 0.001;
		cout<<“\t\t\tGuass 算法“<		cout<<“\n______________________________________________________________“<		cout.width（15）;
		cout<<“k“;

		for（loop = 0;loop < row;loop++）
		{
			cout.width（15）;
			cout<<“X“;
		}

		cout<<“\n______________________________________________________________“<		int index;
		for（index = 0;index < line;index++）
		{
			if（index == 0）
			{
				cout.width（15）;
				cout<			}

			cout.width （15）;
			cout<		}
		cout<		cout.clear （）;
		int w = 1.2;
		do
		{
			double sum = 0;
			double total = 0;
			int index_i= 0;
			int index_j = 0;
			index_j = 0;
			do
			{

				if（index_j < index_i-1）
					do{

						sum += coefficient_matrix[index_i][index_j]*Next_iterator[index_j];
						index_j++;

					}while（index_j < index_i -1）;

				index_j = index_i+1;
				for（;index_j < row;index_j++）
					total += coefficient_matrix[index_i][index_j]*First_iteratorValue[index_j];

				Next_iterator[index_i] = First_iteratorValue[index_i]*（1-w）+w*（constant_Matrix[index_i]-total-sum）/coefficient_matrix[index_i][index_i];
				index_i++;
			}while（index_i < row）;

			int FOR_loop = 0;
			for（;FOR_loop < line;FOR_loop++）
			{  

				First_iteratorValue [FOR_loop] = Next_iterator[FOR_loop];			
			}



			for（index = 0;index < line;index++）
			{
				if（index == 0）
				{
					cout

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2010-03-11 16:38  超松弛\
     文件          46  2010-01-15 07:19  超松弛\data.txt
     文件          43  2010-01-15 08:17  超松弛\data2.txt
     文件          49  2010-01-15 08:18  超松弛\data3.txt
     文件          44  2010-01-15 08:19  超松弛\data4.txt
     文件          26  2010-01-15 08:19  超松弛\data5.txt
     目录           0  2010-03-11 16:38  超松弛\Debug\
     文件        8828  2010-01-15 08:34  超松弛\Debug\BuildLog.htm
     文件       47507  2010-01-15 08:34  超松弛\Debug\Function.obj
     文件        5076  2010-01-15 08:34  超松弛\Debug\main.obj
     文件          67  2010-01-15 08:34  超松弛\Debug\mt.dep
     文件      175104  2010-01-15 08:34  超松弛\Debug\vc90.idb
     文件      225280  2010-01-15 08:34  超松弛\Debug\vc90.pdb
     文件       44032  2010-01-15 08:34  超松弛\Debug\超松弛.exe
     文件         663  2010-01-15 08:34  超松弛\Debug\超松弛.exe.embed.manifest
     文件         728  2010-01-15 08:34  超松弛\Debug\超松弛.exe.embed.manifest.res
     文件         621  2010-01-15 08:34  超松弛\Debug\超松弛.exe.intermediate.manifest
     文件      402592  2010-01-15 08:34  超松弛\Debug\超松弛.ilk
     文件      535552  2010-01-15 08:34  超松弛\Debug\超松弛.pdb
     文件        3457  2010-01-15 08:34  超松弛\Function.cpp
     文件         262  2010-01-15 08:34  超松弛\main.cpp
     文件     1223680  2010-01-15 08:45  超松弛\超松弛.ncb
     文件         883  2010-01-15 08:28  超松弛\超松弛.sln
     文件        8192  2010-01-15 08:45  超松弛\超松弛.suo
     文件        3978  2010-01-15 08:34  超松弛\超松弛.vcproj
     文件        1427  2010-01-15 08:45  超松弛\超松弛.vcproj.PC-20100113JOHO.Administrator.user