测量平差编程(7)

}

int main(int argc, char* argv[]) {

long n,m,k,i,j;

double a[14][15],c[14][13],b[15][13]; FILE *stream;

stream = fopen("矩阵输入.txt","r"); fscanf(stream,"%ld %ld",&n,&m); for (i=0;i<n;i++) { }

fscanf(stream,"%ld %ld",&m,&k); for(i=0;i<m;i++) { }

fclose(stream);

matrixMultiply(a,b,c,4,5,3);

stream = fopen("矩阵计算结果.txt","w"); for (i=0;i<=3;i++) {

for(j=0;j<=2;j++)

fprintf(stream,"%16.7e ",c[i][j]); fprintf(stream,"\n");

for(j=0;j<m;j++) { }

fscanf(stream,"%lf",&a[i][j]);

for(j=0;j<k;j++) { }

fscanf(stream,"%lf",&b[i][j]);