易截截图软件、单文件、免安装、纯绿色、仅160KB

关于列主高斯消元法的C程序实现的问题

下面是我写的关于高斯列主消元法的C程序实现方法。但是求解结果始终不正确啊!检查很久了,还是检查不出来,麻烦哪位大哥给指点下吧!先O(∩_∩)O谢谢了!
程序代码:

C/C++ code:
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define N 4
#define EPS 1e-10

void main()
{
int A[N][N+1]={{5,2,3,2,-1},{2,4,1,-2,5},{1,-3,4,3,4},{3,2,2,8,-6}},x[N];
int i,j,k,Pr,t,total;
for(k=1;k<N;k++)
{
Pr=k-1;
for(i=k;i<N;i++)
{
if(abs(A[i][k-1])>abs(A[Pr][k-1]))
{
Pr=i;
}
}//End For i
if(abs(A[i][k-1])<EPS)
{
printf("主元接近于0,方法失效\n");
exit(0);
}
if(Pr>k-1)
{
for(j=k-1;j<N+1;j++)
{
t=A[k-1][j];
A[k-1][j]=A[Pr][j];
A[Pr][j]=t;
}//End For j
} //End if
for(i=k;i<N;i++)
{
t=A[i][k-1]/A[k-1][k-1];
for(j=k-1;j<N+1;j++)
{
A[i][j]=A[i][j]-t*A[k-1][j];
} //End For j
} //End For i
} //End For k
if(abs(A[N-1][N-1])<EPS)
{
printf("主元接近于0,方法失效\n");
exit(0);
} //End If
x[N


相关问答:

mysql C API求助

写了一段程序
C/C++ code:

..
MYSQL* mysql;
MYSQL_RES* res;
//我想实现一个查看有没有 mysql结构,如果没有就连接。
void connect()
{
mysql_init(mysql);
mysql_realconnect(. ......

在c版块的你们都从事什么工作?

我看学生比较多...
额是业余的……

偶是业余学生

我本来想学c侵入式学习
可觉得学java容易些
所以改为j2ee方向

看来我的估计是准确的, 因为大部分工作的人都用c++或者java, .NET之类...

c/php/p ......

c陷阱与缺陷,说真的看过这本的觉得这书怎么样?

我个人感觉有的地方没必要看,繁琐得狠,用自己总结的方法更好,大家说说自己的意见
我感觉便捷计算与不对称边界完全根据自己的经验来,用不着看哦,大家说说

你看不看那是你的事

其实都是基础的东西啦

......
© 2009 ej38.com All Rights Reserved. 关于E健网联系我们 | 站点地图 | 赣ICP备09004571号