x = rand(100,1);
gamma = randn(100,1) * 0.5;
y = 2 + 3 * x + gamma;

A = ones(100,2);
A(:,2) = x;

[U W V] = svd(A,0);
coef = V*inv(W)*U'*y;

a = coef(1)
b = coef(2)

X = [0:0.1:1];
Y = a + b*X;

plot(x,y,'o',X,Y);

w = diag(W);
ae = (V(1,1)^2)/w(1) + (V(1,2)^2)/w(2)
be = (V(2,1)^2)/w(1) + (V(2,2)^2)/w(2)

[bootstat,bootsam] = bootstrp(100,'corrcoef',x,y);