【正文】 接近的特征值λik。 else return a。) { k=500j。 /* u[501]為初始迭代向量*/ double a,b,c=0。j=min2(i+2,500)。i=500。 } array[jk+2][k]=array[jk+2][k]/array[2][k]。i=500。 } u[i]=u[i]/array[2][i]。 chuzhi(u)。 d=fmifa(u,MatrixC,a)。i++) { chuzhi2(u,i)。若為0，則以此類推。 printf(改變迭代初始向量對λmax計(jì)算結(jié)果的測試如下：\n)。i=39。 printf( SY1103120 朱舜杰\n)。i=0。i=500。j=min2((k+2),500)。i=500。 } //求uk for(i=0。) { k=j。 j++。②矩陣A的行列式可先對矩陣A進(jìn)行LU分解后，detA等于U所有對角線上元素的乘積。二． 源程序includeincludeinclude include includeincludeincludedefine E /*定義全局變量相對誤差限*/int max2(int a,int b) /*求兩個整型數(shù)最大值的子程序*/{ if(ab) return a。 } i++。 j++。i=500。i++) { a[i]=(*rand()/RAND_MAX)。j++) { for(t=max3(0,k2,j2)。i++) { a=a+u[i]*u[i]。i) { for(j=i+1。//將矩陣A轉(zhuǎn)存為MatrixCassignment(MatrixC)。i++) { a=d1+(i*(d501d1))/40。 assignment(MatrixC)。i=500。 chuzhi(u)。 d=mifa(u,MatrixC,0)。j++) { u[i]=array[ij+2][j]*u[j]。 //求yk1 for(i=0。t++) { array[jk+2][k]=array[jt+2][t]*array[tk+2][k]。 for(i=0。 for(j=max2(i2,0