def falseMatchvsM(n_=1000):
  a_ = 200
  theta_ = 15
  s_ = 25

  # 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, 20):

    # Need this otherwise calculation goes to 0
    numSynapsesAfterUnion = numOnBits.subs(n, n_).subs(s, s_).subs(M,
                                                                   M_).evalf()
    if error >= 0.9999999999:
      error = 1.0
    else:
      if M_ <= 2:
        eq3 = subsampledFpFSlow.subs(n, n_).subs(a, a_).subs(theta, theta_)
      else:
        eq3 = subsampledFpF.subs(n, n_).subs(a, a_).subs(theta, theta_)
      error = eq3.subs(s,round(numSynapsesAfterUnion)).evalf()

    print M_,numSynapsesAfterUnion,error

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

  print MList
  print errorList
  return MList,errorList