
N=16;
f=inline('exp(-i*x)+exp(7*i*x)+sin(x)');
Y=fft_coeffs(f,N);

x=(0:0.01:2*pi)';
fx=f(x);
tx=fft_eval(Y,x);

err=norm(fx-tx,inf)

fprintf('\n');
