c语言求解矩阵方程,用C语言求解N阶线性矩阵方程Ax=b的简单解法（附源码）

#include

#include

#include

#define dim 10                              //定义最大的维数10，为防止计算值溢出

double a[dim+1][dim+1],b[dim+1],x[dim+1];   //定义双精度数组

double temp;

double getarray(int n);                     //定义输入矩阵元素的函数

double showarray(int n);                    //定义输出化简系数矩阵过程的函数

int n,i,j,k,p,q;

double main()

{

printf("请输入系数矩阵的阶数n(n<10):");

scanf("%d",&n);

/*判断矩阵阶数是否超过界定值*/

if(n>dim)

{

printf("错误：元数超过初设定的值%d，请重启程序重新输入\n",dim);

exit(0);

}

/*输入系数矩阵和常数矩阵(即增广矩阵)的元素*/

getarray(n);

/*使对角线上的主元素不为0*/

for(j=1;j<=n-1;j++)

{

if(a[j][j]==0)

for(i=j+1;i<=n;i++)

{

if(a[j]!=