def solve_hyperb(xs,ys):
x1,x2,x3 = xs
y1,y2,y3 = ys
a = x1 - x3
b = x1 - x2
c = x1 * y1 - x2 * y2
d = x1 * y1 - x3 * y3
e = y1 - y2
f = y1 - y3
x0 = (a * c - b * d) / ( a * e - b * f )
y0 = (c - x0 * e) / b
a = (y1 - y0) * (x1 - x0)
return (a,x0,y0)
xs = [4,5,6]
ys = [2.6667,2.2500,2.0000]
print(solve_hyperb(xs,ys))
CmRlZiBzb2x2ZV9oeXBlcmIoeHMseXMpOgogICAgeDEseDIseDMgPSB4cwogICAgeTEseTIseTMgPSB5cwogICAgYSA9IHgxIC0geDMKICAgIGIgPSB4MSAtIHgyCiAgICBjID0geDEgKiB5MSAtIHgyICogeTIKICAgIGQgPSB4MSAqIHkxIC0geDMgKiB5MwogICAgZSA9IHkxIC0geTIKICAgIGYgPSB5MSAtIHkzCiAgICB4MCA9IChhICogYyAtIGIgKiBkKSAvICggYSAqIGUgLSBiICogZiApCiAgICB5MCA9IChjIC0geDAgKiBlKSAvIGIKICAgIGEgPSAoeTEgLSB5MCkgKiAoeDEgLSB4MCkKICAgIHJldHVybiAoYSx4MCx5MCkKCnhzID0gWzQsNSw2XQp5cyA9IFsyLjY2NjcsMi4yNTAwLDIuMDAwMF0KCnByaW50KHNvbHZlX2h5cGVyYih4cyx5cykp