tds=[9205,974,274];	
str=[12126,799,272];
tdsb=mean(tds);
strb=mean(str);
tdsc=tds-tdsb;
strc=str-mean(str);
a=sum(tdsc.*strc)/sum(tdsc.^2)
b=strb-a*tdsb