function v = erfintegrand(t);
% ERFINTEGRAND  Computes the coefficients of the 8th degree polynomial
%     interpolant of f(t) = e^{-t^2}.  Input t = (interpolation points).
% Example, to solve 5(a):
%   >> erfintegrand(0.0:0.1:0.8)

A = fliplr(vander(t));  % Vandermonde matrix with preferred order
y = exp(-t.^2);
v = (A \ y')';

