from scipy import stats
x = [2019,2020,2021,2022,2023]
y = [161000,182000,257000,282000,307000]
slope,intercept,r,p,std_err = stats.linregress(x,y)
def myfunc(x):
return slope * x + intercept
print(f'{myfunc(2024):,}')
ZnJvbSBzY2lweSBpbXBvcnQgc3RhdHMKCnggPSBbMjAxOSwyMDIwLDIwMjEsMjAyMiwyMDIzXQp5ID0gWzE2MTAwMCwxODIwMDAsMjU3MDAwLDI4MjAwMCwzMDcwMDBdCgpzbG9wZSxpbnRlcmNlcHQscixwLHN0ZF9lcnIgPSBzdGF0cy5saW5yZWdyZXNzKHgseSkKCmRlZiBteWZ1bmMoeCk6CiAgIHJldHVybiBzbG9wZSAqIHggKyBpbnRlcmNlcHQKCnByaW50KGYne215ZnVuYygyMDI0KTosfScp