function  dopulse (tmax, tw, n, fmax, df)
%function  dopulse (tmax, tw, n, fmax, df);
%Generate time pulse, calculate Fourier Transform and plot;
%assumes an even time signal;
%tmax=max time, tw=pulse width from 0, n=number of time steps;
%fmax=max frequency, df=freq step;
%example:  dopulse (1,.1,100,50,1);

[t,g]=pulse (tmax,tw,n);
subplot(1,2,1);
plot(t,g);
[f,G]=evenft(t,g,fmax,df);
subplot(1,2,2);
plot(f,G);