通过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 (n < 10):); 14. scanf(%d, &n); 15. /*检查矩阵的阶数是否超过预设的限制*/ 16. if (n > dim) 17. { 18. printf(错误：矩阵的维度超过了最初设定的值 %d，请重新启动程序并重新输入。\n, dim); 19. exit(0); 20. } 21. 22. /*输入系数矩阵和常数矩阵（即增广矩阵）的元素*/ 23. getarray(n); 24. 25. /*确保主对角线上的主元不为零*/ 26. for (j = 1; j <= n - 1; j++) { 27 . if (a[j][j] == 0) { //如果主元为零... 继续处理... } //如果主元为零... 继续处理... } //如果主元为零... 继续处理... } //如果主元为零... 继续处理... } //如果主元为零... 继续处理... } //如果主元为零... 继续处理... } //如果主元为零... 继续处理... } //如果主元为零... 继续处理... } //如果主元为零... 继续处理... } //如果主元为零... 继续处理... } //如果主元为零，则交换行元素以消除奇异性。