harmonic.py

# https://www.youtube.com/watch?v=d3mHfqd0VZY
def harmonic(x, y, n=6, epsilon=1e-3):
    r = sqrt(x/y)
    for a in range(1, n):
        for b in range(1, n):
            if abs(r - a/b) < epsilon:
                return True
    return False

for x in np.arange(1.6, 6.4, 0.2):
    for y in np.arange(1.6, 6.4, 0.2):
        status = '✔️' if harmonic(x, y) else '❌'
        #print("{:.2f}{}".format(sqrt(x)/sqrt(y), status), end=' ')
        print("{}".format(status), end=' ')
    print()