guass.c

使用高斯列主消元法实现高次方程组的求解过程

#include<stdio.h>
#include<math.h>
#include<iostream.h>
#include<conio.h>

#define MAX_N 20

void main()
{
	int n;
	int i,j,k;
	int mi;
	double mx,tmp;
	static double a[MAX_N][MAX_N],b[MAX_N],x[MAX_N];
	printf("Input n value dim of Ax=b:");
	scanf("%d",&n);
	if (n > MAX_N )
	{
		printf("The input  ....");
		//return 1;

	}
	if (n <= 0)
	{
		printf("Please input....");
		//return 1;
	}
	printf("Now input matrix a(i,j)|b(i):\n");
	for(i=0;i<n;i++)
	{
		for(j=0;j<n;j++)
		{
			scanf("%lf",&a[i][j]);
	         //printf("%lf\n",a[i][j]);	
		}
		scanf("%lf",&b[i]);
	}

    

    for(i=0;i<n;i++)
	{
		for(j=0;j<n;j++)
		{
		     printf("%f",a[i][j]);
			 if(j==n-1)
				 printf("X%d = ",j+1);
			 else
				 printf("X%d + ",j+1);
		}
		printf("%lf\n",b[i]);
	}


	for(i=0;i<n-1;i++)
	{
		for(j=i+1,mi=i,mx=fabs(a[i][j]);j<n;j++)
		{
			if(fabs(a[j][i])>mx)
			{
				mi=j;
				mx=fabs(a[j][i]);
			}
		}


		if(i < mi)
		{
			tmp=b[i];
			b[i]=b[mi];
			b[mi]=tmp;
			for(j=i;j<n;j++)
			{
				tmp = a[i][j];
			    a[i][j] = a[mi][j];
			    a[mi][j] = tmp;
			}
		}

         for(j=i+1;j<n;j++)
		 {
			 tmp=-a[j][i]/a[i][i];
			 b[j]+=b[i]*tmp;
			 for(k=i;k<n;k++)
			 {
				 a[j][k]+=a[i][k]*tmp;
			 }
		 }
	}

	x[n-1]=b[n-1]/a[n-1][n-1];

	for(i=n-2;i>=0;i--)
	{
		x[i]=b[i];
		for(j=i+1;j<n;j++)
		{
			x[i]-=a[i][j]*x[j];
		}
		x[i]/=a[i][i];
	}
	printf("the roots of this Equation are  x_i:\n");
	for(i=0;i<=n-1;i++)
	{
		printf("X%d=%f\n",i+1,x[i]);
	}
	getch();
}