N=30;
for n = 1:N;
   m = 2 * n - 1;
   x = linspace(-1,1,n);
   V = vander(x); 
   y = linspace(-1,1,m);
   W = vander(y);
   W = W(:,m-n+1:m);      % W is m by n; keep only powers y^0,...,y^{n-1}
   A = W / V;             % = W V^(-1)
   f(n) = norm(A,'inf');
end
np=1:N;
semilogy(np,f,'ko','markersize',2)
hold on
semilogy(np,2.^np./(e*np.*log(np)),'k')
hold off
grid on