code: def mktk03(terms, seed, G2, G3):
      
  # random.seed(n)
      p = 0
      n = seed
      gh = []
      g10, sfact, afact = -18e3, 3.8, 2.4
      g20 = G2 * g10
      g30 = G3 * g10
      alpha = g10/afact
      s1 = s_l(1, alpha)
      s10 = sfact * s1
      gnew = random.normal(g10, s10)
      if p == 1:
          print(1, 0, gnew, 0)
      gh.append(gnew)
      gh.append(random.normal(0, s1))
      gnew = gh[-1]
      gh.append(random.normal(0, s1))
      hnew = gh[-1]
      if p == 1:
          print(1, 1, gnew, hnew)
      for l in range(2, terms + 1):
          for m in range(l + 1):
              OFF = 0.0
              if l == 2 and m == 0:
                  OFF = g20
              if l == 3 and m == 0:
                  OFF = g30
              s = s_l(l, alpha)
              j = (l - m) % 2
              if j == 1:
                  s = s * sfact
              gh.append(random.normal(OFF, s))
              gnew = gh[-1]
              if m == 0:
                  hnew = 0
              else:
                  gh.append(random.normal(0, s))
                  hnew = gh[-1]
              if p == 1:
                  print(l, m, gnew, hnew)
      return gh
docs: generates a list of gauss coefficients drawn from the TK03 distribution