高斯方法解线性方程组c程序.docVIP

高斯消去法和高斯主元消去法解线性方程组 高斯消元法： #includestdio.h #includemath.h main() {int gauss(int n,double a[],double b[]); int i; double a[3][3]={{3,-1,4},{-1,2,-2},{2,-3,-2}}; double b[3]={7,-1,0}; if(gauss(3,a[0][0],b)!=0) for(i=0;i=2;i++) printf(\nx[%d]=%f\n,i,b[i]);} int gauss(int n,double a[],double b[]) {int i,k,j,p,q; double d,t; for(k=0;k=n-2;k++) {d=a[k*n+k]; if(d==0) return(0); for(j=k+1;j=n-1;j++) {p=k*n+j;a[p]=a[p]/d;} b[k]=b[k]/d; for(i=k+1;i=n-1;i++) {for(j=k+1;j=n-1;j++) {p=i*n+j; a[p]=a[p]-a[i*n+k]*a[k*n+j];} b[i]=b[i]-a[i*n+k]*b[k];}} d=a[(n-1)*n+n-1]; if(fabs(d)+1.0==1.0) {printf(fail\n); ret