主元素约当法程序的设计.pptVIP

主元素约当法程序设计 例:2x1-x2+3x3=1 4x1+2x2+5X3=4 ●X1+2x2=7 #define n 3 #define plus 4 #includesstdio. h #includemath.h float aa[n][n={{2,-1,3},{4,2,5} 1,2,0}; float bb[n]=(1, 4,7 float a[n]plus],t,r,e int p, q, L,j, k for (k=0; kn; k++) [p=k; q=l; e=a[k][1; for(j=0; jn: j++) If(fabs(allnfabs(e)) I e=allii p=l; q=j forg=l: j=n+1; j++) °{t=ak][j a[klj=alp[j; april=t; For(=0: In; I++) if((=k)(a[kq=o)) I r=allllqlalkllql for(j=k+1j=n+1:j++) a[j]=a[[]-r*ak]j]; °m[k]=q for(I=0;In;++) {q=m[1];