function [pt,pf,f,V1]=plotpair(t,v1,fdisplay);
npts=length(v1);
pt=figure;plot(t,v1);grid on;
moose=axis;axis([moose(1:2),moose(3)-(moose(4)-moose(3))*0.1,...
		  moose(4)+(moose(4)-moose(3))*0.1]);
xlabel('t, Time, sec');ylabel('Signal Volts');
%  Fourier transform
V1=fftshift(fft(v1)/npts*max(t));
f=fftaxisshift(fftaxis(t)); % needs 2 functions from Chuck DiMarzio
% plot the spectra
test=find(abs(f)<fdisplay);
pf=figure;plot(f(test),real(V1(test)),'b-o',...
	    f(test),imag(V1(test)),'r-x',...
	    f(test),abs(V1(test)),'g-');grid on;
moose=axis;axis([moose(1:2),moose(3)-(moose(4)-moose(3))*0.1,...
		  moose(4)+(moose(4)-moose(3))*0.1]);
xlabel('F, Frequency, Hz');ylabel('V, Volts/Hz');
