function [area] = inter(psd, dt)

% determine the area under the curve of a psd (integral) using the psd curve 
% and freqency resolution.
% Meg Palmsten, USF/USGS, St. Pete, 5/99.

area1 = psd * dt;
line = length(psd);
area = sum(area1(2:line));

% done