def est_e(limit = 20, F=None, E=None):
  x = 1.391925581569705

  if F == None:
    _f = 596.25 #596.25
  else:
    _f = F

  if E == None:
    _e = (((1/(((((1/log(x, pi-0.5))-1))-1)*(x-1)))-(1/((((1/_f)*17)+17)**2))))
  else:
    _e = E

  i = 0
  while i < limit:
    v = log(1/(log(_e)-1))
    _e = ((((1/(((((1/log(x, pi-0.5))-1))-1)*(x-1)))-(1/((((1/_f)*17)+17)**2)))) - (1/(_e**v)))
    _f = 1/((((1/((1/(((((1/log(x, pi-0.5))-1))-1)*(x-1)))-_e))**0.5)/17)-1)
    i = i + 1

  print("done!")
  print(f"estimate e: {_e},  actual: {e},  ratio: {_e/e}, F: {F}")
  return _f, _e

'''
Editted slightly for those who want to do something like

F, E = est_e(1, 1)
and then

i = 0
somelimit = 20
while i < limit:
  F, E = est_e(1, F, E)
  i = i + 1
'''