易截截图软件、单文件、免安装、纯绿色、仅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


相关问答:

jsp中 <c:if test> 相关问题

在查询后将查询出来的值赋给各输入框
<c:if test="${not empty dataValue}">
      fm.SAMPLING_DATE.value=" <c:out value='${dataValue.SAMPLING_DATE}'/ ......

mysql C API求助

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

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

【★C/C++奖励基金-3期★】

C/C++  2009-9专家榜
名次  专家名称  专家分  个人描述 
1  mstlq (面色铁路桥)  5843  因为菜,所以努力学习中 
2  whg01 (及时结贴是美德)&n ......
© 2009 ej38.com All Rights Reserved. 关于E健网联系我们 | 站点地图 | 赣ICP备09004571号