% Example2_6
%
Re = [ 0.1 1 10 10^2 10^3 10^4 10^5 ]';
Nu = [ 0.45 0.84 1.83 5.1 15.7 56.5 245 ]';
n = length(Re);
p = 3;   % number of parameters
y = log10( Nu);
t = log10(Re);
X = [ones(n,1) t t.^2 ];
b_LS = X \ y;
fprintf('   i  ')
fprintf(' %6d',[1:n]);
fprintf('\n')
fprintf('  t_i ')
fprintf(' %6.3f',t);
fprintf('\n')
fprintf('  y_i ')
fprintf(' %6.3f',y);
fprintf('\n')
fprintf('X-tranpose:\n')
for ir = 1:p
    fprintf(' %6.3f',X(:,ir));
    fprintf('\n')
end
fprintf(' b_LS:')
fprintf(' %6.4f',b_LS);
fprintf('\n')



