雅可比,高斯-赛德尔以及sor迭代法

《雅可比,高斯-赛德尔以及sor迭代法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、昨天闲来无事,就把雅可比,高斯-赛德尔和SOR迭代法解方程组的程序写了.把主要程序发了,可能有谁用的着./**********雅可比迭代法**********/voidJacobi(floatxk[]){ inti,j; floatt=0.0; floattt=0.0; float*xl; xl=(float*)malloc(sizeof(float)*(n+1)); for(i=1;i

2、 xl[i]=(b[i]-t)/a[i][i]; } t=0.0; for(i=1;ic) {   if(t<=p)   printf("Reachthegivenprecision!");  else     printf("overthemaximalcount!");  printf("Countnumberis%d",k)

3、; } else  if(t>p)  {   k++;   Jacobi(xk);  }  else  {    printf("Reachthegivenprecision!");   printf("Countnumberis%d",k);  }}/**********高斯-赛德尔迭代法**********/voidGauss_Siedel(floatxk[]){ inti,j; floatt=0.0; floattt=0.0; float*xl; xl=(float*)malloc(sizeof(float)*(n+1)); for(i=1;i

4、i++) {  t=0.0;  tt=0.0;  for(j=1;jc) {   if(t<=p)   printf

5、("Reachthegivenprecision!");  else     printf("overthemaximalcount!");  printf("Countnumberis%d",k); } else  if(t>p)  {   k++;   Gauss_Siedel(xk);  }  else  {    printf("Reachthegivenprecision!");   printf("Countnumberis%d",k);  }}/**********SOR迭代法**********/voidSOR(floa

6、txk[]){ inti,j; floatt=0.0; floattt=0.0; float*xl; xl=(float*)malloc(sizeof(float)*(n+1)); for(i=1;i

7、i]-xk[i]);   tt=tt*tt;   t+=tt;  } t=sqrt(t); for(i=1;ic) {   if(t<=p)   printf("Reachthegivenprecision!");  else     printf("overthemaximalcount!");  printf("Countnumberis%d",k); } else  if(t>p)  {    k++;   S