【文章內(nèi)容簡(jiǎn)介】 resh()。 sleep(2)。 endwin()。 exit(EXIT_SUCCESS)。 } 二、 C 語言程序清單 1 void grns(u,g,r,n,a) int n。 double u,g,*r,a[]。 { int i,k,m。 double s,w,v,t。 s=。 w=。 v=。 for (k=0。 k=n1。 k++) { t=。 for (i=1。 i=12。 i++) {*r=(*r)*w+v。 m=(int)(*r/s)。 *r=*rm*s。 t=t+(*r)/s。 } 16 a[k]=u+g*()。 } return。 } 2 include include main() { int i,j。 double u,g,r,a[50]。 r=。 u=。 g=。 printf(\n)。 grns(u,g,amp。r,50,a)。 for (i=0。 i+9。 i++) { for (j=0。 j+4。 j++) printf(% ,a[i*5+j])。 printf(\n)。 } printf(\n)。 } 3 include include int lluu(a,n,l,u) int n。 double a[],l[],u[]。 {int i,j,k,w,v,ll。 for (k=0。k=n2。k++) {ll=k*n+k。 if (fabs(a[ll])+==) {printf(fail\n)。 return(0)。} for (i=k+1。i=n1。i++) {w=i*n+k。a[w]=a[w]/a[ll]。} for (i=k+1。i=n1。i++) {w=i*n+k。 for (j=k+1。j=n1。j++) {v=i*n+j。 a[v]=a[v]a[w]*a[k*n+j]。 } } } for (i=0。i=n1。i++) {for (j=0。j=i。j++) 17 {w=i*n+j。l[w]=a[w]。u[w]=。} w=i*n+i。 l[w]=。u[w]=a[w]。 for(j=i+1。j=n1。j++) {w=i*n+j。l[w]=。u[w]=a[w]。} } return(1)。 } 4 include include main() {int i,j。 static double l[4][4],u[4][4]。 static double a[4][4]={{,},{,},{,},{,}}。 i=lluu(a,4,l,u)。 if (i!=0) {printf(MAT L IS:\n)。 for(i=0。i=3。i++) {for(j=0。j=3。j++) printf(%,l[i][j])。 printf(\n)。 } printf(\n)。 printf(MAT U IS:\n)。 for(i=0。i=3。i++) {for(j=0。j=3。j++) printf(%,u[i][j])。 printf(\n)。 } printf(\n)。 } } 5 include include int chlk(a,n,m,d) int n,m。 double a[],d[]。 {int i,j,k,u,v。 if ((a[0]+==)||(a[0])) { printf(fail\n)。 return(2)。 } a[0]=sqrt(a[0])。 18 for (j=1。j=n1。j++) a[j]=a[j]/a[0]。 for (i=1。i=n1。i++) { u=i*n+i。 for (j=1。j=i。j++) {v=(j1)*n+i。 a[u]=a[u]a[v]*a[v]。 } if ((a[u]+==)||(a[u])) {printf(fail\n)。 return(2)。 } a[u]=sqrt(a[u])。 if (i!=(n1)) {for (j=i+1。j=n1。j++) {v=i*n+j。 for (k=1。k=i。k++) a[v]=a[v]a[(k1)*n+i]*a[(k1)*n+j]。 a[v]=a[v]/a[u]。 } } } for (j=0。j=m1。j++) { d[j]=d[j]/a[0]。 for (i=0。i=n1。i++) {u=i*n+i。v=i*m+j。 for (k=1。k=i。k++) d[v]=d[v]a[(k1)*n+i]*d[(k1)*m+j]。 d[v]=d[v]/a[u]。 } } for (j=0。j=m1。j++) {u=(n1)*m+j。 d[u]=d[u]/a[n*n1]。 for (k=n1。k=1。k) { u=(k1)*m+j。 for (i=k。i=n1。i++) {v=(k1)*n+i。 d[u]=d[u]a[v]*d[i*m+j]。 } v=(k1)*n+k1。 d[u]=d[u]/a[v]。 } } return(2)。 } 6 19 include include main() {int i。 static double a[4][4]={{,},{,},{,},{,}}。 static double d[4][2]={{,},{,},{,},{,}}。 if (chlk(a,4,2,d)0) for (i=0。i=3。i++) printf(x(%d)=%,%\n,i,d[i][0],d[i][1])。 } 7 include include int gsdl (a,b,n,x,eps) int n。 double a[],b[],x[],eps。 { int i,j,u,v。 double p,t,s,q。 for (i=0。i=n1。i++) {u=i*n+1。 p=。 x[i]=。 for (j=0。j=n1。j++) if (i!=j) {v=i*n+j。 p=p+fabs(a[v])。} if (p=fabs(a[u])) { printf(fail\n)。 return(1)。} } p=eps+。 while (p=eps) {p=。 for (i=0。 i=n1。 j++) { t=x[i]。 s=。 for (j=0。 j=n1。 j++) if (j!=i) s=s+a[i*n+j]*x[j]。 x[i]=(b[i]s)/a[i*n*i]。 q=fabs(x[i]t)/(+fabs(x[i]))。 if (qp) p=q。 } } return(1)。 } 8 include 20 include main() {int i。 double eps。 static double a[4][4]={{,},{,},