顺序高斯消元算法

#include<iostream>

#include<math.h>

using namespace std;

#define MAX 100

double A[MAX+1][MAX+1];

double B[MAX+1];

double X[MAX+1];

double e;

int n;

void OrderGaussianElimination()

{

int i,j,k;

double T;

//消元

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

{

if( fabs(A[k][k])<=e ) { printf("Error!\n"); return ;}

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

{

T=A[i][k]/A[k][k];

B[i]-=T*B[k];

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

A[i][j]-=T*A[k][j];

}

//回代

if( fabs(A[n][n])<=e ) { printf("Error!\n"); return ;}

X[n]=B[n]/A[n][n];

double S_Aij_Xj;

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

{

S_Aij_Xj=0;

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

S_Aij_Xj+=A[i][j]*X[j];

X[i]=(B[i]-S_Aij_Xj)/A[i][i];

}

void print(double X[])

{

int i;

printf("X[]:\n");

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

printf("%f\n",X[i]);

}

int main()

{

int i,j;

printf("n:");

scanf("%d",&n);

printf("A[][]:\n");

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

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

scanf("%lf",&A[i][j]);

printf("B[]:\n");

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

scanf("%lf",&B[i]);

printf("e:");

scanf("%lf",&e);

OrderGaussianElimination();

print(X);

return 0;

}

posted on 2009-05-16 16:32 wyiu 阅读(226) 评论(0) 编辑收藏引用所属分类: 数值分析

只有注册用户登录后才能发表评论。
【推荐】100%开源！大型工业跨平台软件C++源码提供，建模，组态！

相关文章: 全主元高斯消元法列主元高斯消元算法顺序高斯消元算法

网站导航: 博客园 IT新闻 BlogJava 知识库博问管理

常用链接

留言簿(1)

随笔分类(84)

随笔档案(100)

向高手学习

搜索

积分与排名

最新评论

阅读排行榜

评论排行榜