数值分析实验报告冶金工业.docxVIP

数值分析实验报告冶金工业 数值分析实验报告 课题一：解线性方程组的直接方法 实验目的： 1、 通过该课题的实验，体会模块化结构程序设计方法 的优点； 2、 运用所学的计算方法，解决各类线性方程组的直接 算法； 3、 提高分析和解决问题的能力，做到学以致用； 4、 通过三对角形线性方程组的解法，体会稀疏线性方 程组解法的特点。 实验过程： 实验代码： #include/z/z #i nclude ttinclude usingn amespacest d; //Gauss 法 voidlzy(double**a, double*b, i ntn) inti, j, k; double 1, x[10], te mp; for(k=0 ;k{ for (j=k , i=k; j{ if(j==k) temp=fabs (a[j] [k]); elseif (temp{ temp =fabs (a[j] [k]); ? ? l=j; } } if (temp=0) { coutret urn; } else { for (j=k;j { temp=a[k] [j]; a[k] [j]=a[i] [j]; a [i] [j]=tem p; } temp=b [k]; b[k]=b[i]; b[i]=te mp; for (i=k+l;i { l=a[i][k]/a[k][k]; for