function [fs,f]=myspectrum(s,dt)
% MYSPECTRUM makes amplitude spectrum of time series
% [fs,f]=myspectrum(s,dt)
% s seismogram
% dt time step
% f frequency vector

disp(' Warning! Taper set to 1000! ')
s=taper(s,-1000);
s=taper(s,1000);

fs=fft(s);
nt=length(fs);
fny=1./(2*dt);
finc=fny/(nt/2);
f(1:nt/2)=(1:nt/2)*finc;

% end