Hi,
Should the argument of the arccos not be (E_0-E)/2t? Or am I missing something?
you are correct, thank you for reporting - fixed now