clear 

n=10;

t=zeros(n,1);



t=[1:10]';



x=[ ones(n,1) t];



k=size(x,2);



y=[4 5 6 7 2 3 3 1 8 5]';



b=inv(x'*x)*x'*y

e=y-x*b;

ESS=e'*e;

TSS=(y-mean(y))'*(y-mean(y));

R2=1-ESS/TSS

modR2=1-(ESS/(n-k))/(TSS/(n-1))

s2=ESS/(n-k);



s2b=diag(inv(x'*x))*s2;

sb=s2b.^0.5



tval=b./s2b