from scipy import stats
y = [1,3,5,7,12,17,27,40,55,70,85,117,116,]
x = [2009,2010,2011,2012,2013,2014,2015,2016,2017,2018,2019,2020,2021]
slope,intercept,r,p,std_err = stats.linregress(x,y)
def myfunc(x):
return slope * x + intercept
print(myfunc(2018))
ZnJvbSBzY2lweSBpbXBvcnQgc3RhdHMKCnkgPSBbMSwzLDUsNywxMiwxNywyNyw0MCw1NSw3MCw4NSwxMTcsMTE2LF0KeCA9IFsyMDA5LDIwMTAsMjAxMSwyMDEyLDIwMTMsMjAxNCwyMDE1LDIwMTYsMjAxNywyMDE4LDIwMTksMjAyMCwyMDIxXQoKc2xvcGUsaW50ZXJjZXB0LHIscCxzdGRfZXJyID0gc3RhdHMubGlucmVncmVzcyh4LHkpCgpkZWYgbXlmdW5jKHgpOgogICByZXR1cm4gc2xvcGUgKiB4ICsgaW50ZXJjZXB0CgpwcmludChteWZ1bmMoMjAxOCkp