Newton interpolation coefficient

This function finds the coefficients for Newton interpolation polynomial

% [b0 b1 b2...] = NewtonInterpolation(x,y) where vector [b0 b1 b2...] is
% the coefficients of Newton interpolation polynomial:
%
% N(x) = b0+b1(x-x0)+b2(x-x0)(x-x1)+...
%
% Note that both x, y should be 1-D vector with the same dimension.