文档介绍：数值分析实验报告篇一: 1 2 3 4 5篇二: 课题一:解线性方程组的直接方法 : 1、通过该课题的实验,体会模块化结构程序设计方法的优点; 2、运用所学的计算方法,解决各类线性方程组的直接算法; 3、提高分析和解决问题的能力,做到学以致用; 4、通过三对角形线性方程组的解法,体会稀疏线性方程组解法的特点。 : 实验代码: #include"" #include"" #include usingnamespacestd; //Gauss法 voidlzy { inti,j,k; doublel,x[10],temp; for { for { if temp=fabs; elseif) { temp=fabs; i=j; } } if { cout return; } else { for { temp=a[k][j]; a[k][j]=a[i][j]; a[i][j]=temp; } temp=b[k]; b[k]=b[i]; b[i]=temp; } for { l=a[i][k]/a[k][k]; for a[i][j]=a[i][j]-l*a[k][j]; b[i]=b[i]-l*b[k]; } } if { cout return; } x[n-1]=b[n-1]/a[n-1][n-1]; for { temp=0; for temp=temp+a[i][j]*x[j]; x[i]=/a[i][i]; } for { printf; printf; } } //平方根法 voidpfg { inti,k,m; doublex[8],y[8],temp; for { temp=0; for temp=temp+pow; if return; a[k][k]=pow,); for { temp=0; for temp=temp+a[i][m]*a[k][m];a[i][k]=/a[k][k]; } temp=0; for temp=temp+a[k][m]*y[m]; y[k]=/a[k][k]; } x[n-1]=y[n-1]/a[n-1][n-1]; for { temp=0; for temp=temp+a[m][k]*x[m]; x[k]=/a[k][k]; } for { printf,x[i]); printf; } } //追赶法 voidzgf { inti; doublea0[10],c[10],d[10],a1[10],b1[10],x[10],y[1