线性方程组数值解法-列主元高斯消去法

 

//列主元高斯消去法
/*  2(x1)+4(x2)+(x3)=4
    2(x1)+6(x2)-(x3)=10
	(x1)+5(x2)+2(x3)=2
*/
#include<stdio.h>
#include<math.h>
void main(){
	void ColPivot(float *,int,float[]);
	int i;
	float x[3];
	float c[3][4]={
		   2,4,1,4,
           2,6,-1,10,
           1,5,2,2};
	ColPivot(c[0],3,x);
	for(i=0;i<=2;i++)
   printf("x[%d]=%f\n",i,x[i]);
}
void ColPivot(float *c,int n,float x[])
{
	int i,j,k,t;
	float p;
	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!=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[j]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
	}
}