标题:解线性方程组用的顺序Gauss消元法的编程,麻烦高手指点一下哪里有问题?
只看楼主
ainiwcc1
Rank: 1
等 级:新手上路
帖 子:7
专家分:0
注 册:2009-10-10
结帖率:66.67%
 问题点数:0 回复次数:0 
解线性方程组用的顺序Gauss消元法的编程,麻烦高手指点一下哪里有问题?
#include"iostream.h"
 void main()
 {
     int n,i,j,k;
     int s=0;
     cin>>n;
     float a[100][100];
     float b[100];
     float x[100];
     for(i=0;i<n;i++)
         for(j=0;j<n;j++)
         
             cin>>a[i][j];
            
   
         for(i=0;i<n;i++)
             cin>>b[i];//输入数组;
         for(k=0;k<n-1;k++)
            for(i=k+1;i<n;i++)
               for(j=k+1;j<n;j++)
               { a[i][j]=a[i][j]-a[i][k]*a[k][j]/a[k][k];//对a[i][j]进行计算;
                   b[i]=b[i]-a[i][k]*b[k]/a[k][k];//对b[i]进行计算;
               }
      x[n-1]=b[n-1]/a[n-1][n-1];
      for(k=n-2;k>=0;k--)
      for(j=k+1;j<n;j++)
      {s=s+a[k][j]*x[j];
         x[k]=(b[k]-s)/a[k][k];}//回代;
 
         cout<<"方程组的解为:"<<endl;
         for(i=0;i<n;i++)
             cout<<x[i]<<endl;
}
搜索更多相关主题的帖子: 麻烦 顺序 线性方程 Gauss 
2009-10-10 12:53



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-288043-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.055056 second(s), 8 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved