import numpy as np
from time import sleep
"""
c = a*b
a = c/b
b = c/a
a-b = c/b - c/a
a-b = c*a-c*b/(a*b)
a*b*(a-b) = c*a - c*b

d = a-b
b = a-d

a = c/b = c/(a-d)
a²-d*a - c = 0

(a - d/2) - d²/4 = c

d = c/()

11*13
"""
a = 123
b = 431
c = a*b
d = 1#a-b
#a²-d*a - c = 0
x = np.linspace(0,1,100)
while True:
    p = np.roots((1,-d,-c))
    err = (x-p[0])*(x-p[1])
    d += 0.01 * (err[1] - err[1]//1)
    print(np.sum(err**2))
    sleep(0.001)