dist2_diff_solve(664515.9375 + 0j, 3998106.0, 0.0536749947614598, 
    -71340.5618328239, 23709057427.2661)
[664515.7418921513 + 3.552713678800501e-15j, 664600.9266076663 + 
    5.329070518200751e-15j, 664564.8069210749 - 1.4210854715202004e-14j]
sympy lambdify RuntimeWarning: invalid value encountered in double_scalars