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
print(myfunc(200))
ZnJvbSBzY2lweSBpbXBvcnQgc3RhdHMKCnggPSBbMSwyLDMsNF0KeSA9IFsxNTAwLDMwMDAsNDUwMCw1NTAwXQoKc2xvcGUsaW50ZXJjZXB0LHIscCxzdGRfZXJyID0gc3RhdHMubGlucmVncmVzcyh4LHkpCgpkZWYgbXlmdW5jKHgpOgogICByZXR1cm4gc2xvcGUgKiB4ICsgaW50ZXJjZXB0CgpwcmludChteWZ1bmMoMjAwKSk=