Guass列主元消去法及其应用

《Guass列主元消去法及其应用》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、说明：深圳大学实验报告实验课程名称：计算机数值方法实验项目名称：Guass列主元消去法及其应用学院：计算机软件学院专业：计算机科学与技术报告人：学号：班级：04指导教师：实验时间：2010.4.20实验报告提交时间：2010.5.3一、实验目的与要求了解Guass列主元消去法及其应用，掌握利用Guass列主元消去法解方程，解决实际问题二、模型求解0.x+0.4y+0.3z+60000=x0.3x+0.1y+0.2z+30000=y0.2x+0.1y+0.2z+50000=z三、源代码：#include#includ

2、e#include#defineN20#defineEPS1.0e-4doublea[N][N],b[N],c[N][N],d[N][N];voidswap(double&a,double&b)//换行{doubletemp;temp=a;a=b;a=temp;}voidgauss(doublea[][N],doubleb[N],intn)//gauss求解{inti,j,k,l;doubletemp,m;for(k=0;k<=n-2;k++){temp=fabs(a[k][k]);//选主元l=

3、k;for(i=k+1;itemp))//fabs为取绝对值函数。包含在头文件math.h中{temp=fabs(a[i][k]);l=i;}}if(fabs(a[l][k])

4、r(j=k;j=0;i--){for(j=n-1;j>i;j--)b[i]-=a[i][j]*b[j];b[i]/=a[i][i];}}voidmain(){inti,j,n;doubletemp=0.0;printf("请输入要求解的方程数:");//矩阵维数为nsc

5、anf("%d",&n);printf("请输入矩阵a[%d][%d]:",n,n);for(i=0;i

6、]:",n,n);for(i=0;i

7、1、报告内的项目或内容设置，可根据实际情况加以调整和补充。2、教师批改学生实验报告时间应在学生提交实验报告时间后10日内。