Fixed numerical issue for points exactly on the ellipse axes (square root of difference of nearly equal numbers). Improved error detection and handling. Increased robustness by increasing number of allowable iteratons and increasing the number of initial values tested.