code: def edf_mtotdev(N, m, alpha):
      
      assert(alpha in [2, 1, 0, -1, -2])
      NIST_SP1065_table8 = [(1.90, 2.1), (1.20, 1.40), (1.10, 1.2), (0.85, 0.50), (0.75, 0.31)]
      #(b, c) = NIST_SP1065_table8[ abs(alpha-2) ]
      (b, c) = NIST_SP1065_table8[abs(alpha-2)]
      edf = b*(float(N)/float(m))-c
      print("mtotdev b,c= ", (b, c), " edf=", edf)
      return edf
docs: Equivalent degrees of freedom for Modified Total Deviation
      
          NIST SP1065 page 41, Table 8