实验2 LU分解求线性方程组

一．实验名称：LU分解法解线性方程组

二、实验目的：

（1）学会对矩阵进行LU分解。

（2）学会用LU分解法解线性代数方程组。

三、实验要求

（1）按照题目要求完成实验内容

（2）写出相应的C语言程序

（3）给出实验结果

（4）写出相应的实验报告

四、实验题目

1.  LU分解(Doolittle分解)

编写对矩阵A进行LU分解的程序，其中

2.在以上程序的基础上进一步编写出利用LU分解的方法求解线性方程组的程序。

五、实验结果

输出L、U及方程的解x。

#include

#include

int main(){

int n=4,i,k,j,m=n+1;

float a[m][m],u[m][m],l[m][m],t,t1,x[m],y[m],b[m];

for(i=1;i<=n;i++)

for(j=1;j<=n;j++)

{

}

for(i=1;i<=n;i++)

 l[i][i]=1;

b[1]=2;

b[2]=8;

b[3]=4;

b[4]=9;

printf("请输入%d阶矩阵:\\n",n);

for(i=1;i<=n;i++)

for(i=1;i<=n;i++){

y[1]=b[1]; 

x[n]=y[n]/u[n][n]; 

for(i=n-1;i>=1;i--)

for(i=1;i<=n;i++){

}

printf("\\n");

for(i=1;i<=n;i++){

}

for(i=1;i<=n;i++){

}

}