from scipy import stats
x = [1,3,5,7,12,17,27,40,55,70,85,117,116,]
y = [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(2008))
ZnJvbSBzY2lweSBpbXBvcnQgc3RhdHMKCnggPSBbMSwzLDUsNywxMiwxNywyNyw0MCw1NSw3MCw4NSwxMTcsMTE2LF0KeSA9IFsyMDA5LDIwMTAsMjAxMSwyMDEyLDIwMTMsMjAxNCwyMDE1LDIwMTYsMjAxNywyMDE4LDIwMTksMjAyMCwyMDIxXQoKc2xvcGUsaW50ZXJjZXB0LHIscCxzdGRfZXJyID0gc3RhdHMubGlucmVncmVzcyh4LHkpCgpkZWYgbXlmdW5jKHgpOgogICByZXR1cm4gc2xvcGUgKiB4ICsgaW50ZXJjZXB0CgpwcmludChteWZ1bmMoMjAwOCkp