【文章內容簡介】 G1[i][m])fabs(GtG1[i+1][m])) ik=i+1。 else 。 } if(ik!=m) { for(i=m。i=k。i++) { temp=GtG1[m][i]。 GtG1[m][i]=GtG1[ik][i]。 GtG1[ik][i]=temp。 } temp=Gt[m][l]。 Gt[m][l]=Gt[ik][l]。 Gt[ik][l]=temp。 } for(i=m+1。i=k。i++) { q=GtG1[i][m]/GtG1[m][m]。 for(j=m。j=k。j++) { GtG1[i][j]=GtG1[i][j]q*GtG1[m][j]。 } Gt[i][l]=Gt[i][l]q*Gt[m][l]。 } } D[k][l]=Gt[k][l]/GtG1[k][k]。 for(m=k1。m=0。m) { temp=0。 for(i=m+1。i=k。i++) { temp+=D[i][l]*GtG1[m][i]。 } D[m][l]=(Gt[m][l]temp)/GtG1[m][m]。 } } for(i=0。i=k。i++) { for(j=0。j=k。j++) { temp=0。 for(m=0。m=20。m++) { temp+=A[i][m]*D[j][m]。 } C[i][j]=temp。 } } sigma=0。/*歸零，開始計算sigma值*/ for(i=0。i=10。i++) { for(j=0。j=20。j++) { p[i][j]=0。 for(i1=0。i1=k。i1++) { for(j1=0。j1=k。j1++) { p[i][j]+=C[i1][j1]*pow(x[i],i1)*pow(y[j],j1)。 } } sigma+=pow(p[i][j]z[i][j],2)。 } } printf(k=%d sigma=%.12e\n,k,sigma)。 k=k+1。 } k。 printf(達到精度要求時的k和sigma值：\n)。 printf(k=%d sigma=%.12e\n,k,sigma)。 kz=k。 /*定義為全局變量，便于duibi()調用*/}void duibi() /*對比??和p逼近效果的子程序*/{ int i,j,i1,j1。 double p[8][5]。 for(i=0。i=7。i++) { for(j=0。j=4。j++) { x[i]=*(i+1)。 y[j]=+*(j+1)。}} /*重新輸入節(jié)點*/ Newton(x,y)。 fpeccz(t,u)。 /*求解??（x*，y*）*/ for(i=0。i=7。i++) /*求解p（x*，y*）*/ { for(j=0。j=4。j++) { p[i][j]=0。 for(i1=0。i1=kz。i1++) { for(j1=0。j1=kz。j1++) { p[i][j]+=C[i1][j1]*pow(x[i],i1)*pow(y[j],j1)。 } } printf(x[%d]=%.6f y[%d]=%.6f\n,i+1,x[i],j+1,y[j])。 printf(f(x[%d],y[%d])=%.12e\n,i+1,j+1,z[i][j])。 printf(p(x[%d],y[%d])=%.12e\n,i+1,j+1,p[i][j])。 printf(σ=%.12e\n,z[i][j]p[i][j])。 }/*數表xi*,yi*,?? (xi*,yi*),p（xi*,yi*）*/} }void main(){ int i,j。 for(i=0。i=10。i++) { for(j=0。j=20。j++) { x[i]=*i。 y[j]=+*j。}} /*輸入節(jié)點*/ Newton(x,y)。 fpeccz(t,u)。 for(i=0。i=10。i++) { for(j=0。j=20。j++) { printf(x[%d]=%.6f y[%d]=%.6f z[%d][%d]=%.12e \n,i,x[i],j,y[j],i,j,z[i][j])。 }} /*數表：xi,yi, ?? (xi,yi)*/ qmnh()。 for(i=0。i=kz。i++) { for(j=0。j=kz。j++) { printf(C[%d][%d]=%.12e \n,i,j,C[i][j])。 }}/*數表：C[r][s]*/ duibi()。 }三、運行結果輸出數表數表：xi,yi,??(xi,yi)x[0]= y[0]= z[0][0]=x[0]= y[1]= z[0][1]=x[0]= y[2]= z[0][2]=x[0]= y[3]= z[0][3]=x[0]= y[4]= z[0][4]=x[0]= y[5]= z[0][5]=x[0]= y[6]= z[0][6]=x[0]= y[7]= z[0][7]=x[0]= y[8]= z[0][8]=x[0]= y[9]= z[0][9]=x[0]= y[10]= z[0][10]=x[0]= y[11]= z[0][11]=x[0]= y[12]= z[0][12]=x[0]= y[13]= z[0][13]=x[0]= y[14]= z[0][14]=x[0]= y[15]= z[0][15]=x[0]= y[16]= z[0][16]=x[0]= y[17]= z[0][17]=x[0]= y[18]= z[0][18]=x[0]= y[19]= z[0][19]=x[0]= y[20]= z[0][20]=x[1]= y[0]= z[1][0]=x[1]= y[1]= z[1][1]=x[1]= y[2]= z[1][2]=x[1]= y[3]= z[1][3]=x[1]= y[4]= z[1][4]=x[1]= y[5]= z[1][5]=x[1]= y[6]= z[1][6]=x[1]= y[7]= z[1][7]=x[1]= y[8]= z[1][8]=x[1]= y[9]= z[1][9]=x[1]= y[10]= z[1][10]=x[1]= y[11]= z[1][11]=x[1]= y[12]= z[1][12]=x[1]= y[13]= z[1][13]=x[1]= y[14]= z[1][14]=x[1]= y[15]= z[1][15]=x[1]= y[16]= z[1][16]=x[1]= y[17]= z[1][17]=x[1]= y[