# include<stdio.h>
# include<math.h>
# define delta 1e-6
#define N 100
int main()
{
int i,j,t,r,n,u,c=0;
float p,L,max,s;
float X[N];
float a[N][N+1];
printf("请输入系数矩阵的阶数:\n");
scanf("%d",&n);
printf("输入的增广矩阵系数,中间用空格隔开:\n");
for(i=0;i<n;i++)
for(j=0;j<n+1;j++)
scanf("%f",&a[i][j]);
printf("你输入的增广矩阵为:\n");
for(i=0;i<n;i++)
{
for(j=0;j<n+1;j++)
printf("%f\t",a[i][j]);
printf("\n");
}
for(j=0;j<n-1;j++)
{
{
max=fabs(a[j][j]);
r=j;
}
for(i=j+1;i<n;i++)
if(fabs(a[i][j])>max)
{
max=a[i][j];
r=i;
}
if(fabs(a[i][j])<delta)
printf("矩阵奇异");
for(t=j;t<n+1;t++)
{
p=a[j][t];
a[j][t]=a[r][t];
a[r][t]=p;
}
for(i=j+1;i<N;i++)
{
L=a[i][j]/a[j][j];
for(t=j;t<n+1;t++)
a[i][t]=a[i][t]-L*a[j][t];
}
}
printf("输出原方程的解为:\n");
X[n-1]=a[n-1][n]/a[n-1][n-1];
for(i=n-2;i>=0;i--)
{
s=a[i][n];
for(j=i+1;j<n;j++)
s=s-a[i][j]*X[j];
X[i]=s/a[i][i];
}
for(u=0;u<n;u++)
for(j=0;j<n+1;j++)
{
printf("%12f",a[u][j]);
c++;
if(c%(n+1)==0)
printf("\n");
}
for(i=0;i<n;i++)
{
printf("x(%d)=%.4f\n",i+1,X[i]);
if(i==n-1)
printf("\n");
}
system("pause");
return 0;
}
高斯消去法
最新推荐文章于 2022-05-15 17:47:40 发布