雅可比迭代和高斯赛德尔迭代

 

 

 

 

     算法伪代码如下，有错误的话请给我指出来：

 

//假设所有a[i][i]都不是0，为0的话还不知道怎么办
		输入系数于数组a[][];
		输入常数于数组b[];
		
		//计算新系数于数组c[][]
		for(int i=0; ij)
					c[i][j] = -a[i][j]/a[i][i];
				else
					continue;
			}
		}	
		//迭代100次，自己加上精度判断吧
		//从向量形式可以看出，初试解为全0的话第一次解为c[i][n-1]，如果为了省着一次时间的话可以先把
		//这个复制给解数组 root[],再仔细看看会发现每次迭代的话都需要加上c[i][n-1]
		while(T--) {
			//k行
			for(int j=0; jk)
						temp += c[j][k]*root[k];
					else if(j