def falseMatchDendritevsMPatterns(n_=1000):
  a_ = 30
  theta_ = 10
  s_ = 20

  # Arrays used for plotting
  MList = []
  errorList = []

  print "\n\nn=%d, a=%d, theta=%d, s=%d" % (n_,a_,theta_,s_)

  error = 0.0
  for M_ in range(1, 40):

    # Number of bits active in At
    atAfterUnion = numOnBits.subs(n, n_).subs(s, a_).subs(M,M_).evalf()
    if error >= 0.9999999999:
      error = 1.0
    else:
      if M_ <= 8:
        eq3 = subsampledFpFSlow.subs(n, n_).subs(a, atAfterUnion).subs(theta, theta_)
      else:
        eq3 = subsampledFpF.subs(n, n_).subs(a, atAfterUnion).subs(theta, theta_)
      error = eq3.subs(s,s_).evalf()

    print M_,atAfterUnion,error

    MList.append(M_)
    errorList.append(error)

  print MList
  print errorList
  return MList,errorList