为什么手算的和编程算,不一样了,请帮助看看,程序哪里错了,谢谢
#include "stdio.h"#include "math.h"
#include <malloc.h>
float *gaoshi(float*c,int n)
{
int i,j,k,t;
float *x,p;
x=(float*)malloc(n*sizeof(float));
for (i=0;i<=n-2;i++)
{
k=i;
for (j=i+1;j<=n-1;j++)
if(fabs(*(c+j*(n+1)+i))>fabs(*(c+k*(n+1)+i)))
{ k=j;
if(k!=j)
for (j=i;j<=n;j++)
{
p=*(c+i*(n+1)+j);
*(c+i*(n+1)+j)=*(c+k*(n+1)+j);
*(c+k*(n+1)+j)=p;
}
}
else
{
for (j=i+1;j<=n-1;j++)
{
p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
for(t=i;t<=n;t++)
*(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
}
}
}
for (i=n-1;i>=0;i--)
{
for(j=n-1;j>=i+1;j--)
{ (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));}
x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)*i));
}
return x;
}
void main()
{
float a[4][5]={1,1,0,3,4,
2,1,-1,1,1,
3,-1,-1,3,-3,
-1,2,3,-1,4};
float *x;
//x=a[0];
int i;
x=gaoshi(a[0],4);
for (i=0;i<=3;i++)
printf(" x[%d]=%f\n",i+1,x[i]);
}我要算的是个方程
x1+x2+3*x4=4
2*x1+x2-x3+x4=1
3*x1-x2-x3+3*x4=-3
-x1+2*x2+3*x3-x4=4我是用高斯列主元素消去法解的,觉的没什么错,但就是不对,请帮助看看,谢谢