毕业论文论文范文课程设计实践报告法律论文英语论文教学论文医学论文农学论文艺术论文行政论文管理论文计算机安全
您现在的位置: 毕业论文 >> 报告 >> 正文

计算方法列主元高斯消去法求解方程组 第2页

更新时间:2007-9-26:  来源:毕业论文

 if(det!=0)

 {

    for(i=1;i<=n;i++)

           printf("[%d]=%f\n",i,x[i]);

    printf("\n");

 }

}

 int solve(a,x)

 double a[n+1][n+2],x[n+1];

 {

  int i,j,k,r;

  double c;

  for(k=1;k<=n;k++)

  {

   r=k;

   for(i=k;i<=n;i++)

          if(fabs(a[i][k])>fabs(a[r][k]))

                 r=i;

          if(fabs(a[r][k])+1.0==1.0)

          {

                 printf("\ndet=0 fail!");

                 return 0;

          }

   if(r!=k)

          for(j=k;j<=n+1;j++)

          {

           c=a[k][j]; a[k][j]=a[r][j]; a[r][j]=c;

          }

   for(i=k+1;i<=n;i++)

   {

    c=a[i][k]/a[k][k];

       for(j=k+1;j<=n+1;j++)

              a[i][j]=a[i][j]-c*a[k][j];

   }

  }

  if(fabs(a[n][n])+1.0==1.0)

  {

   printf("\n det=0. fail!");

   return (0);

  }

  for(k=n;k>=1;k--)

  {

         x[k]=a[k][n+1];

         for(j=k+1;j<=n;j++)

                x[k]=x[k]-a[k][j]*x[j];

         x[k]=x[k]/a[k][k];

  }

  return (1);

 }

上一页  [1] [2] [3] 下一页

计算方法列主元高斯消去法求解方程组 第2页下载如图片无法显示或论文不完整,请联系qq752018766
设为首页 | 联系站长 | 友情链接 | 网站地图 |

copyright©751com.cn 辣文论文网 严禁转载
如果本毕业论文网损害了您的利益或者侵犯了您的权利,请及时联系,我们一定会及时改正。