C语言实现单纯形算法详解与应用

/*本程序是单纯形法，参考书目:河海大学<吴凤平>运筹学方法与应用*/   
  #include<stdio.h>   
  #include<math.h>   
  #define   X   5   
  #define   Y   7   
  void   xi_max(int   *m2,int   *mn1,float   *c,int   *is,int   *ir,int   *j0,float   (*a)[X][Y])   
  {   
      int   j;   
      *c=0;   
      for(j=1;j<=*is;j++)     
          {   
              if((*a)[*ir][j]-*c>0)   
          {   
              *c=(*a)[*ir][j];   
              *j0=j;   
          }   
            }   
  }   
    
              /*****************   参数说明   **********************/   
              /*       m_约束方程个数(基变量个数),n_非基变量个数       */   
              /*       m2-m+2整个变量,(*a)[X][Y]存放初始数据               */   
              /*       (*k)[]存放基变量脚标,(*x)存放基变量最优值       */   
              /*********   http://happyyangxu.home.sunbo.net*******/     
  int   xi_sm(int   m,int   n,int   m2,int   mn1,int   l1,float   (*a)[X][Y],   
                int(*k)[],float(*x)[])   
  {   
      int   i,m1,mn,j0,i0,j;   
      float   c,g;   
      m1=m+1;