function [y] = vander(v);
% ======================================================================
% Vandermonde-Matrix
% ----------------------------------------------------------------------
global nfc 
nfc = nfc+1;
xa=0;
xb=1;

n=length(v);
x=linspace(xa,xb,n)';
e=ones(n,1);

a(1,1)=n;
for i=2:2*n-1
    a(i,1)=e'*x.^(i-1);
end;

for i=1:n
    y(i,1)=a(i:n+i-1,1)'*v;
end;

return;
