欢迎来到天天文库
浏览记录
ID:39500951
大小:48.00 KB
页数:3页
时间:2019-07-04
《矩阵三角分解法》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、班级:研1104学号:2011020109姓名:宁建业#include#defineN4main(){doublea[N][N]={111,-1,0,0,-1,112,-2,0,0,-2,113,-2,0,0,-3,114};doublel[N][N],u[N][N],f[N]={3,1,0,-5};doublex[N],y[N];inti,j;l[0][0]=a[0][0];u[0][0]=1;for(i=0;ii)l[i][j]=0;if(j2、j)u[i][j]=1;if(j-i>=2)u[i][j]=0;if(i-j==1)l[i][j]=a[i][j];if(i-j>=2)l[i][j]=0;}}for(i=0;i3、;}printf("");}printf("");printf("u矩阵:");for(i=0;i=0;i--){if4、(i==N-1)x[N-1]=y[N-1];elsex[i]=y[i]-u[i][i+1]*x[i+1];printf("x[%d]=%.4f",i+1,x[i]);}printf("");}
2、j)u[i][j]=1;if(j-i>=2)u[i][j]=0;if(i-j==1)l[i][j]=a[i][j];if(i-j>=2)l[i][j]=0;}}for(i=0;i3、;}printf("");}printf("");printf("u矩阵:");for(i=0;i=0;i--){if4、(i==N-1)x[N-1]=y[N-1];elsex[i]=y[i]-u[i][i+1]*x[i+1];printf("x[%d]=%.4f",i+1,x[i]);}printf("");}
3、;}printf("");}printf("");printf("u矩阵:");for(i=0;i=0;i--){if
4、(i==N-1)x[N-1]=y[N-1];elsex[i]=y[i]-u[i][i+1]*x[i+1];printf("x[%d]=%.4f",i+1,x[i]);}printf("");}
此文档下载收益归作者所有