LU分解法解线性方程组（C语言）

#include #include #define N 10 //矩阵大小范围 /* * 使用已经求出的x，向前计算x（供getx()调用） * float a[][] 矩阵U * float x[] 方程组解 * int i 解的序号（数组X元素序号） * int n 矩阵大小 * return 公式中需要的和 */ float getmx(float a[N][N], float x[N], int i, int n) { float mx = 0; int r; for(r=i+1; rN){ printf("矩阵过大！/n"); continue; } flag=0; } printf("请输入系数矩阵值：/n"); for(i=0; i=0; i--) { x[i] = getx(u,y,x,i,n); } /*显示数组X*/ printf("/n/n数组X：/n"); for(i=0; i