资源描述:
《算法报告实例》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、《数值分析》实验报告成绩:班级:信息12・2小组成员:方MM(121246402XX)卜MM(121246402XX)谭GG(121246402XX)2014年12月25日广东石油化工学院理学院信息系目录一.课题:线性方程组的矩阵分解法1.运行代码32.实验题目43.运行结果5二、课题:线性方程组的迭代法71.运行代码72•实验题目10103・运行结果.三、课题:拉格朗日插值法151.运行代码152.实验题目16173.运行结果.心得20、课题:线性方程组的矩阵分解法1、运行代码:#include<
2、stdio.h>#include#defineMAX_N20intmain(intargc,char*argv[]){intn;//未知数个数inti,j,k;staticdoublea[MAX_N][MAX_N],b[MAX_N],x[MAX_N],y[MAX_N];staticdoublel[MAX_N][MAX_N],u[MAX_N][MAX_N];printf("lnputnvalue(dimofAx=b):");scanf("%dH,&n);if(n>MAX_N){pr
3、intf("TheinputnislargerthanMAX_N,pleaseredefinetheMAX_N.rT);return1;}if(n<=0){printf("Pleaseinputanumberbetween1and%d.",MAX_N);return1;}//{{程序输入printf("Nowinputthematrixa(i,j),i,j=0,%d:",n-1);for(i=0;i4、j]);pr
5、intf("Nowinputthematrixb(i),i=0,%d:",n-1);for(i=0;ivn;i++)scanf(n%ir,&b[i]);//}}程序输入for(i=0;ivn;i++)u[i][i]=l;//for(k=0;k6、k]U]=a[k][j];for(i=0;i<=k-l;i++)u[k]U]-=(l[k][i]*u[i]Ul);u[k]U]/=l[k][k];}}for(i=0;i=0;i-)//解UX=Y{x[i]=y[i];for(j=i+l;j7、x_i=n);//输出结果for(i=0;i4)、41014202400-1(4,321)6012_4_2-830-2兀211-101305兀3=1710-492£-1022-117_玉._5_3)、(1"丄丄1)9-4(1丄0)32111314232-7-6012).(100,1)_7200027300037-105)、00-17-3000-36_0000-20__11_0兀2250兀3230兀41()-2
8、兀568兀638(1,2,3,4,5,6)3、运行结果1)、Inputnvalue(dimofAx=b):4Nowinputthematrixa(i>j),i,jx0,3:4100tZ2°324-1)0-14^owinputthematrixb(i),i=0,•・・,3:19口20152Solve・・・x_i=4.0993792.6024842・745342L1S63352)、口EACYuY*nCYu¥dnbinwwtemp.exeInputnvalue(dimofAx=b):4Nowinpu
9、tthematrixa(i,j),i,j=0,・・・,3:1113-6142032-71-60114Nowinputthematrixb(i),i=0,・・・,3:8Solve・.・x_i=1.0000000.000000-0.0000001.0000003)、4)、'E:CYuYanCYuYanbinwwtemp.exeInputnvalue(dimofAx=b):3Nowinputthematrixa(i,j),i,j=0,•…2:9-41-432128N