用C语言求解N阶线性矩阵方程Ax=b的简单解法

1. #include 2. #include 3. #include 4. #define dim 10 //定义最大的维数10，为防止计算值溢出 5. double a[dim+1][dim+1],b[dim+1],x[dim+1]; //定义双精度数组 6. double temp; 7. double getarray(int n); //定义输入矩阵元素的函数 8. double showarray(int n); //定义输出化简系数矩阵过程的函数 9. int n,i,j,k,p,q; 10. double main() 11. { 12. 13. printf("请输入系数矩阵的阶数n(ndim) 17. { 18. printf("错误：元数超过初设定的值%d，请重启程序重新输入\n",dim); 19. exit(0); 20. } 21. 22. /*输入系数矩阵和常数矩阵（即增广矩阵）的元素*/ 23. getarray(n); 24. 25. /*使对角线上的主元素不为0*/ 26. for(j=1;j<=n-1;j++) 27. { 28. if(a[j][j]==0) 29. for(i=j+1;i<=n;i++) 30. { 31. if(a[i][j]!=0) 32. { 33. /*交换增广矩阵的第i行与第j行的所有元素*/