t1 = (2*x1*x3 + 2*y1*y3 + 2*z1*z3 + sqrt(-(4*x3^2 + 4*y3^2 + 4*z3^2)*(x1^2 + y1^2 + z1^2 - v^2) + (-2*x1*x3 - 2*y1*y3 - 2*z1*z3)^2)) / (2*x3^2 + 2*y3^2 + 2*z3^2)
t2 = (2*x1*x3 + 2*y1*y3 + 2*z1*z3 - sqrt(-(4*x3^2 + 4*y3^2 + 4*z3^2)*(x1^2 + y1^2 + z1^2 - v^2) + (-2*x1*x3 - 2*y1*y3 - 2*z1*z3)^2)) / (2*x3^2 + 2*y3^2 + 2*z3^2)