%
%  w03cond.m -- djm -- 23 sept 2010
%

N = 21;

%  polynomial fitting matrix
M = vander((1:N)/N);

%  random unit vectors
reset(RandStream.getDefaultStream);
aj = 2*rand(N,1)-1;
%aj = floor(1e5*aj)*1e-5;

yj = M*aj;

%  vector of coefficients
Cj = M\yj;

[aj Cj]

