[[a == 1000 * (r1 + sqrt(r1 ^ 2 + 2000 * r1 - 1000000)) / (r1 + sqrt(r1 ^ 2 +
    2000 * r1 - 1000000) + 1000), b == -1 / 2 * r1 - 1 / 2 * sqrt(r1 ^ 2 + 
    2000 * r1 - 1000000) + 500, c == r1], [a == 1000 * (r2 - sqrt(r2 ^ 2 + 
    2000 * r2 - 1000000)) / (r2 - sqrt(r2 ^ 2 + 2000 * r2 - 1000000) + 1000
    ), b == -1 / 2 * r2 + 1 / 2 * sqrt(r2 ^ 2 + 2000 * r2 - 1000000) + 500,
    c == r2]]
Simplest way to solve mathematical equations in Python