Improving this implementation of Bancroft's method

Question

Background

I\'ve written an algorithm (with much assistance from SO) to solve the 3-dimensional Time Difference of Arrival (TDoA) multi-