#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我是用高斯列主元素消去法解的，觉的没什么错，但就是不对，请帮助看看，谢谢

谢谢大家了，不用了
我知哪里不对了