计算方法-解线性方程组的直接法实验报告(2)

}

}

solve(l,u,b,x,n); printf(\解为:\\n\ for(i=1;i<=n;i++)

printf(\

（2）

（3）

#include

#include

#include

#define N 20

double a[N], b[N], c[N-1], f[N], r[N]; int n;

void LUDecompose();

// LU分解

void backSubs(); // 回代

void main() {

printf(\请输入方程的维数n＝\ scanf(\ getchar(); if(n>N||n<=0) { }

printf(\输入下三角元素\\n\ printf(\输入%d个a值: \ for (int i=1; i

scanf(\

printf(\由于该维数过于犀利, 导致程序退出!\ return;

}

getchar();

printf(\输入主对角线元素\\n\ printf(\输入%d个b值: \ for (i=0; i

scanf(\

getchar();

printf(\输入上三角元素\\n\ printf(\输入%d个c值: \ for (i=0; i

scanf(\

getchar();

printf(\输入%d个方程组右端项: \\n\ for (i=0; i

scanf(\

getchar(); LUDecompose(); backSubs();

printf(\线性方程组的解为: \\n\ for (i=0; i

printf(\

void LUDecompose(){ }

void backSubs(){

f[0]=f[0]/b[0]; for(int i=1; i

f[i]=(f[i]-r[i]*f[i-1])/b[i]; c[0]=c[0]/b[0]; for(int i=1;i

b[i]=b[i]-r[i]*c[i-1]; r[i]=a[i];

b[i]=b[i]-r[i]*c[i-1]; c[i]=c[i]/b[i];

f[n-1]=f[n-1]; for(i=n-2;i>=0;i--)

f[i]=f[i]-c[i]*f[i+1];

}

7．【实验结果与分析总结（含运行结果截图）】