555.txt

高斯列主元消去法,求解其次线性方程组

//3.高斯列主元消去法,求解其次线性方程组

//C/C++ code
#include<stdio.h>
#include <math.h>
#define N 20
int main()
{
int n,i,j,k;
int mi,tmp,mx;
float a[N][N],b[N],x[N];
printf("\nInput n:");
scanf("%d",&n);
if(n>N)
{ printf("The input n should in(0,N)!\n");
getch();
return 1;
}
if(n<=0)
{ printf("The input n should in(0,N)!\n");
getch();
return 1;
}
printf("Now input a(i,j),i,j=0...%d:\n",n-1);
for(i=0;i<n;i++)
{ 
   for(j=0;j<n;j++)
   scanf("%f",&a[i][j]);
}
printf("Now input b(i),i,j=0...%d:\n",n-1);
for(i=0;i<n;i++)
   scanf("%f",&b[i]);
for(i=0;i<n-2;i++)
{ 
   for(j=i+1,mi=i,mx=fabs(a[i][j]);j<n-1;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];
}
for(i=0;i<n;i++)
   printf("Answer:\n x[%d]=%f\n",i,x[i]);
getch();
return 0;
}