t=100+[.01 .02 .04]';
H=[10.10272 10.09367 10.09141]';
A=[ones(size(t)) 100.01-t (100.01-t).*(100.02-t)];
c=A\H;
H10003=c(1)+c(2)*(100.01-100.03)+c(3)*(100.01-100.03)*(100.02-100.03)