integer nT
      double precision beta_h(maxT),prob_limit
      common /free/ beta_h,prob_limit,nT