本文介绍了利用C语言编写程序来解决N阶线性矩阵方程Ax=b的一种简便算法和实现方式，旨在为编程爱好者与工程技术人员提供参考。
1. 包含头文件
2. 包含头文件
3. 包含头文件
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. 输出提示信息请输入系数矩阵的阶数n(n<10):
14. 读取用户输入的整数值，存入变量n中
15. //判断矩阵阶数是否超过设定值
16. 如果 n > dim
17. {
18. 输出错误消息：元数超过初设定的值%d，请重启程序重新输入\n ，其中dim为定义的最大维数
19. 退出程序
20. }
21.
22. //调用函数，输入系数矩阵和常数矩阵（即增广矩阵）的元素
23. 调用getarray(n)
24.
25. //确保主对角线上的主元不为零
26. 对于j从1到n-1
27. {
28. 如果 a[j][j] == 0
29. 则对于i从 j+1 到 n
30. {
31. 如果 a[i][j] != 0
32. 则交换增广矩阵的第 i 行与第 j 行的所有元素
5星
