fmt=:[: ;:inv DP&$: : (4 :0)&.>
  x{.deb (x*2j1)":y
)

root=: ln@] exp@% [

epsilon=: 1r9^DP