% Set wavelet name. 
    wname = 'db2';

% Plot wavelet function. 
    [phi,psi,xval] = wavefun(wname,7);
    subplot(211); plot(xval,psi); title('Wavelet'); 

% Compute and plot wavelet integrals approximations 
% on a dyadic grid. 
    [integ,xval] = intwave(wname,7); 
    subplot(212); plot(xval,integ); 
    title(['Wavelet integrals over [-Inf x] ' ... 
        'for each value of xval']);