编程求解方程组.docVIP

用c语言编程求线性方程组的解 悬赏分：30 - 解决时间：2009-5-27 23:58 例如 A=｛1 1｝B=｛1｝ 求 AX=B的 x的解 1 2 2 在线等 问题补充： 答的对有加分。加分不低于30分哦 提问者： winwyj - 三级 最佳答案 #include stdlib.h #include math.h #include stdio.h int rgauss(n,a,b) int n; double a[],b[]; { int *js,l,k,i,j,is,p,q; double d,t; js=malloc(n*sizeof(int)); l=1; for (k=0;k=n-2;k++) { d=0.0; for (i=k;i=n-1;i++) for (j=k;j=n-1;j++) { t=fabs(a[i*n+j]); if (td) { d=t; js[k]=j; is=i;} } if (d+1.0==1.0) l=0; else { if (js[k]!=k) for (i=0;i=n-1;i++) { p=i*n+k; q=i*n+js[k]; t=a[p]; a[p]=a[q]; a[q]=t; } if (is!=k) { for (j=k;j=n-1;j++) { p=k*n+j; q=is*n+j; t=a[p]; a[p]=a[q]; a[q]=t; } t=b[k]; b[k]=b[is]; b[is]=t; } } if (l==0) { free(js); printf(fail\n); return(0); } d=a[k*n+k]; 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) { free(js); printf(fail\n); return(0); } b[n-1]=b[n-1]/d; for (i=n-2;i=0;i--) { t=0.0; for (j=i+1;j=n-1;j++) t=t+a[i*n+j]*b[j]; b[i]=b[i]-t; } js[n-1]=n-1; for (k=n-1;k=0;k--) if (js[k]!=k) { t=b[k]; b[k]=b[js[k]]; b[js[k]]=t;} free(js); return(1); } 徐世良老师的.原地工作全选主元的高斯消元.用的符号有点诡异,而且将矩阵以线性表输入的. 我也写过一个,不过估计你不太喜欢用我的,这个用来解b也是矩阵的. #include stdio.h #include stdlib.h #include math.h void swap(double *a,double *b){ double temp; temp=*a; *a=*b; *b=temp; } void printMatrix(int n,int m,double **a){ int i,j; for (i=0;in;i++){ for (j=0;jm;j++) printf(%e\t,a[i][j]); printf(\n); } } int rgauss(int n,int m,double **a,double **b){//,double *det){ int *colex,rowex,s,i,j,k; double d,t; // *det=1; s=0; colex=malloc(n*sizeof(int)); for (k=0;kn;k++){ d=0.0; for (i=k;in;i++){ for (j=k;jn;j++){ t=fabs(a[i][j]); if (td){ d=t; colex[k]=j; rowex=i; } } } if (d+1.0==1.0){ free(colex); printf(Singular Matrix!\n); return 0; } else{ if (colex[k]!=k){ f