import numpy as np
R = np.array([100000.0, 1538.8, 72.0])
T = np.array([298.15, 423.15, 573.15])
Y = 1.0 / T
L = np.log(R)

X = np.column_stack((np.ones(3), L, L**3))
ABC = np.linalg.solve(X, Y)

A, B, C = ABC
print(f"A = {A:.10e}")
print(f"B = {B:.10e}")
print(f"C = {C:.10e}")