function y=expchirp(fs,f0,f1,t)
% expchirp - exponential chirp
%
% expchirp ::
%    F:real   ~'sampling frequency',
%    real     ~'initial frequency',
%    real     ~'final frequency',
%    D:nonneg ~'duration'
% -> [[1,D*F]].

y=chirp(0:(1/fs):t,f0,t,f1,'logarithmic',-90);
