from scipy import stats
x = [1,2,3,4]
y = [1500,3000,4500,5500]
slope,intercept,r,p,std_err = stats.linregress(x,y)
def myfunc(x):
return slope * x + intercept
result = myfunc(572)
print(f'{result:,.0f}')
ZnJvbSBzY2lweSBpbXBvcnQgc3RhdHMKCnggPSBbMSwyLDMsNF0KeSA9IFsxNTAwLDMwMDAsNDUwMCw1NTAwXQoKc2xvcGUsaW50ZXJjZXB0LHIscCxzdGRfZXJyID0gc3RhdHMubGlucmVncmVzcyh4LHkpCgpkZWYgbXlmdW5jKHgpOgogICByZXR1cm4gc2xvcGUgKiB4ICsgaW50ZXJjZXB0CgpyZXN1bHQgPSBteWZ1bmMoNTcyKQpwcmludChmJ3tyZXN1bHQ6LC4wZn0nKQ==