function v=vanderint(x,y,u)

tic

A=vander(x);
c=A\y.';
v=polyval(c,u);

toc