mc;

% Compute time-series
Nsecs = 1e3; % Number of seconds desired
fs    = 10;

s = ao(plist('tsfcn', 'sin(2*pi*1.137*t)', 'fs', fs, 'Nsecs', Nsecs));
n = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'Nsecs', Nsecs));
b = n + s;

% Make FFT
tic
bxx = fft(b);
toc

% Make DFT
tic
dxx = dft(b);
toc

% Make DFT at a particular frequency
dsx = dft(b, plist('f', 1.137));

% Plot
iplot(bxx, dxx, dsx, plist('Markers', {'none', 'none', 'o'}));

iplot(bxx./dxx)