function a = InterpV(x,y)
%  a = InterpV(x,y)
%  Computes coefficients of interpolating polynomial
%  in power form, using the Vandermonde matrix.
%
%  Input:  x - vector of distinct points
%          y - vector of data at the x-points
%          
%  Output: a - coefficients of interpolating
%              polynomial, a(1) is leading coefficient.

V = vander(x);
a = V \ y;