05.cpp

高思列主元消去法是计算方法的一种方法..这是一个要做的实验....希望可以帮到大家好好理解它

#include <stdio.h>
#include <math.h>
#define MAX_N 20

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("实验五  列主元高斯消元法解 Ax=b\n\n");
	printf("\nInput n value(dim of Ax=b):");
	scanf("%d",&n);
	if (n>MAX_N)
	{
		printf("The input n is larger than MAX_N,please redefine the MAX_N.\n");
		return 1;
	}
	if (n<=0)
	{
		printf("Please input a number between 1 and %d.\n",MAX_N);
		return 1;
	}
	printf("Now input the matrix a(i,j),i,j=0,...,%d: \n",n-1);
	for (i=0;i<n;i++)
		for(j=0;j<n;j++)
		    scanf("%lf",&a[i][j]);
	printf("Now input the matrix b(i),i=0,...m%d: \n",n-1);
	for(i=0;i<n;i++) scanf("%lf",&b[i]);
	for(i=0;i<n-1;i++)
	{
		for(j=i+1,mi=i,mx=fabs(a[i][i]);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[j]+b[i]*tmp;
			for(k=i;k<n;k++) a[j][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]=x[i]-a[i][j]*x[j];
	    x[i]=x[i]/a[i][i];
		
	}
	printf("Solve...x_i= \n");
	for (i=0;i<n;i++) printf("%f\n",x[i]);
	return 0;
	scanf("%d",&mi);
}