高斯-赛德尔迭代法的算法及程序设计.docxVIP

1 1 1 1 高斯-赛德尔迭代法的算法及程序设计 设方程组Ax =b的系数矩阵的对角线元素(i =1,2,|山n) , M为迭代次数容许的最大 值，；为容许误差。 1取初始向量令k=0. 2对i=1,2,…,n计算 3如果则输出结束；否则执行4 4如果则不收敛，终止程序；否则，转 2 源程序： #i nclude stdio.h #in elude math.h #defi ne N 600 void main() { int i; double x[4]; double c[4][5]={10,-1,2,0,-11,0,8,-1,3,-11,2,-1,10,0,6,-1,3,-1,11,25}; void GaussSeidel(double *,i nt,double[]); GaussSeidel(c[0],4,x); for(i=0;i=3;i++) prin tf(x[%d]=%f\n,i,x[i]);} void GaussSeidel(double *a,i nt n, double x[]) { int i,j,k=1; double d,dx,eps; for(i=0;i=3;i++) while(1) {eps=0; for(i=0;i=3;i++) { d=0