C++实现的GPS网络平差程序设计

#include<iostream.h>    
#include<fstream.h>    
#define MAX 100    
#define PI 3.14159265358979312    
#define rou (180.0*60*60/PI)    
#include<math.h>    
#include<string.h>    
   
//***************************矩阵求逆函数说明**************************************    
int inverse(double C[][MAX],double B[][MAX],int n);   
//define the function of inverseing;    
//return 0 means mat can't be inversed or;    
//return 1 means mat inversed correctly, and B is inversed mat;    
   
//***************************矩阵乘积函数说明**************************************    
void AXB(double A[][MAX],double B[][MAX],double C[][MAX],int m,int n,int k);   
void AXB(double a,double A[][MAX], double aA[][MAX],int m,int n);    
void AXB(double A[][MAX],double B[][1],double C[][1],int m,int n);   
//define the time fuction of mats or mat and a number    
   
//***************************矩阵转置函数说明**************************************    
void AT(double A[][MAX],double AH[][MAX],int m,int n);   
void AT(double A[][1],double AH[][MAX],int m);   
void AT(double A[][MAX],double AH[][1],int m);   
//define the fuction to turn-over a mat     
   
//***************************平差计算相关函数**************************************    
void ATPA(double A[][MAX],double P[][MAX],double ATPA[][MAX],int m,int n);   
void ATPL(double A[][MAX],double P[][MAX],double L[][1],double ATPL[][1],int m,int n);   
double VPV(double V[][1],double P[][MAX],int m);   

//define initial fuctions of doadj()    
   
//**************************矩阵显示***********************************************    
void matdis(double A[][MAX],int n,int m);   
void matdis(double A[][1],int n);   
//define two fuctions to output a mat to screen      
   
//**********************通用平差相关结构与函数*************************************    
struct adj;   
void ksetadj(adj &a);     //input data from keyboard     
int fsetadj(adj &aa,char name[20]);  //input data from a file     
int  doadj(adj &a);       //adjusting a adj     
int rubust(adj &a);       //抗差估计    
void adjdis(adj &aa);     //output a adj results to screen     
int foutadj(adj &aa, char name[20]); //outputdate to a file    
//define adj and it related fuctions    
   
//************************************************************************************************    
void AXB(double A[][MAX],double B[][MAX],double C[][MAX],int m,int n,int k)   
{   
   for(int i=0;i<m;i++)   
       for(int j=0;j<k;j++)   
       {   
           C[i][j]=0;   
          for(int l=0;l<n;l++)   
          C[i][j]+=A[i][l]*B[l][j];   
       }   
}   
//************************************************************************************************