用列主元消去法解方程组*/
（1）
#include <stdio.h>
#include <math.h>
#include <alloc.h>
void main()
{
int i;
float *x;
float c[4][5]={{1,1,0,3,4},
{2,1,-1,1,1},
{3,-1,-1,3,-3},
{-1,2,3,-1,4}
};
float *ColPivot(float *,int);
x=ColPivot(c[0],4);
clrscr();
for(i=0;i<=3;i++)
printf("x[%d]=%f\n",i+1,x
);
getch();
}
float *ColPivot(float *c,int n)
{
int i,j,t,k;
float *x,p;
x=(float *)malloc(n*sizeof(float));
for(i=0;i<=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!=i)
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;
}
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=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
}
return x;
}
