import numpy as np
from time import sleep

a = 123
b = 431
c = a*b
#a²-d*a - c = 0 where d = a - b

#a² - (a-b)*a - c
#b*a-c

#a² - d*a - c där d_a = 1

b = 100
a = 100
err = 1
while True:
    err,a,b = np.sort(abs(np.roots([err,1,c//(a//1)-c//(b//1),-c])))
    #b += 0.01 * err
    print(np.sum(err**2))
    sleep(0.001)