1;
## bug with constant y

pkg load statistics

rand("seed",1)
A = rand()
n = 5;

X = [ones(n,1),(1:n)'];
y = A * ones(n,1);

[b1, bint1, r1, rint1, stats1] = regress (y, X)



