C++实现水准网与导线网平差程序

#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];   
       }   
}   
//************************************************************************************************   
   
void AXB(double A[][MAX],double B[][1],double C[][1],int m,int n)   
{   
for(int i=0;i<m;i++)   
       for(int j=0;j<1;j++)   
       {