def calc_pep_index(A,E):
  R = E
  for i in xrange(1000):
    R1 = .975*(A*R)+.025*E
    if norm(R1-R) < some_tolerance:
      break
    R = R1
  else:
    raise Exception("no convergence")
  return R