Guass列主元消去法及其应用.doc

PAGE 实验目的与要求 了解Guass列主元消去法及其应用，掌握利用Guass列主元消去法解方程，解决实际问题 模型求解 0.x+0.4y+0.3z+60000=x 0.3x+0.1y+0.2z+30000=y 0.2x+0.1y+0.2z+50000=z 源代码： #include stdio.h #include stdlib.h #include math.h #define N 20 #define EPS 1.0e-4 double a[N][N],b[N],c[N][N],d[N][N]; void swap(double a,double b)//换行 { double temp; temp=a; a=b; a=temp; } void gauss(double a[][N],double b[N],int n)//gauss求解 { int i,j,k,l; double temp,m; for(k=0;k=n-2;k++) { temp=fabs(a[k][k]); //选主元 l=k; for(i=k+1;in;i++) { if(fabs(a[i][k]temp))//fabs为取绝对值函数。包含在头文件math.h中 { temp=fabs(a[i][k]); l=i; }