单向后方交会实验报告c++

?{ for(intj=0;j<6;j++) { D[i][j]=m0*Asum[i][j]; cout<<setw(10)<<D[i][j]; f6<<setw(10)<<D[i][j]; } cout<<endl; f6<<endl;?}?cout<<"所得中误差为"<<endl;?for(i=0;i<6;i++) cout<<sqrt(D[i][i])<<endl;?f6.close();getchar();?return0;}voidinverse(doublec[n][n]){?inti,j,h,k;?doublep;?doubleq[n][12];?for(i=0;i<n;i++)?for(j=0;j<n;j++) q[i][j]=c[i][j];?for(i=0;i<n;i++)?for(j=n;j<12;j++)?{ if(i+6==j) q[i][j]=1; else q[i][j]=0;?}?for(h=k=0;k<n-1;k++,h++)?for(i=k+1;i<n;i++)?{ if(q[i][h]==0) continue; p=q[k][h]/q[i][h]; for(j=0;j<12;j++) { q[i][j]*=p; q[i][j]-=q[k][j]; }?}?for(h=k=n-1;k>0;k--,h--)?for(i=k-1;i>=0;i--)?{ if(q[i][h]==0) continue; p=q[k][h]/q[i][h]; for(j=0;j<12;j++) { q[i][j]*=p; q[i][j]-=q[k][j]; }?}?for(i=0;i<n;i++)?{ p=1.0/q[i][i]; for(j=0;j<12;j++) q[i][j]*=p;?}?for(i=0;i<n;i++)?for(j=0;j<n;j++)