% seismic.m loads and graphs the seismic data file aBGSOZIP1999.09.21.14.

load aBGSOZIP1999.09.21.18 -ascii  % load data
s = aBGSOZIP1999;                  % give the data a shorter name
t=(1:length(s))/600;

figure(1)
plot(t,s)
xlabel('Time in minutes')
title('Seismic data, 2:00 PM to 3:00 PM on Tuesday, September 21, 1999')
zoom on                            % allow zooming by clicking on the graph

figure(2)
h=[1 1]/sqrt(2);            % multiresolution coefficients for Haar wavelet
                            % put the h you want to use second
showdwt(s',h,600);          % 600 means 600 data points per time unit (minutes)