t = (2005:2013)';x = [130756;131448;132129;132246;133474;134100;134735;135404;136072];p=fittype('poly1') f=fit(t,x,p) plot(f,t,x);
