function y = EvalFourier(a,b,x)
%
% Evaluate Fourier Series at given points.
%
n = length(a)-1;
a = a(:);
b = b(:);
x = x(:);
y = a(1) + cos(x*(1:n))*a(2:n+1)+sin(x*(1:n-1))*b(1:n-1);



