#include<stdio.h>
#include<math.h>
#include<conio.h>
#define NUM 20
#define Enter 13 /*换行符*/
void exchRow(int r,int k); /*矩阵行交换函数*/
float maxCoef(int k); /*求最大系数*/
void message(); /*用于显示提示信息的函数*/
float a[NUM][NUM+1]; /*用于存放方程系数矩阵及向量*/
int maxI,n;
main()
{
float x[NUM]; /*用于存放方程解*/
float s;
int k,i,j;
printf("\n用高斯列主元消元法解线性方程组");
printf("\n请输入方程组的维数:n=");
scanf("%d",&n);
printf("\n输入系数矩阵A和向量b:");
for(i=1;i<=n;i++)
{
printf("\n输入a[%d][1]-a[%d][%d]系数及向量b[%d]:",i,i,n,i); /*输入每一行的系数和向量*/
for(j=1;j<=n+1;j++) /*将输入的数据存入数组*/
scanf("%f",&a[i][j]);
}
for(k=1;k<=n;k++)
{
if(maxCoef(k)==0) /*判断方程是否为线性方程*/
{
printf("\n该方程非线性!");
message();
}
else if(maxI!=k)
exchRow(maxI,k);
for(i=k+1;i<=n;i++)
for(j=k+1;j<=n+1;j++)
a[i][j]=a[i][j]-a[k][j]*a[i][k]/a[k][k];
}
x[n]=a[n][n+1]/a[n][n];
for(k=n-1;k>=1;k--)
{
s=0;
for(j=k+1;j<=n;j++)
s=s+a[k][j]*x[j];
x[k]=(a[k][n+1]-s)/a[k][k];
}
printf("\n求得方程解为:");
for(i=1;i<=n;i++)
printf("\nx%d=%f",i,x[i]);
message();
}
void exchRow(int r,int k) /*矩阵行交换函数*/
{
int i;
for(i=1;i<=n;i++)
a[0][i]=a[r][i];
for(i=1;i<=n+1;i++)
a[r][i]=a[k][i];
for(i=1;i<=n+1;i++)
a[k][i]=a[0][i];
}
float maxCoef(int k) /*求最大系数*/
{
int i;
float t=0;
for(i=k;i<=n;i++)
if(fabs(a[i][k]>t))
{
t=(fabs(a[i][k]));
maxI=i;
}
return t;
}
void message() /*用于显示提示信息的函数*/
{
printf("\n\n继续运算请按Enter,退出请按其他任意键!\n");
if(getch()==Enter)
main();
}
运行后这么多错误,怎么解决,万分感谢!!