datain = 'acoustic.wav';
% read in wav sample
[ x, fs, N ] = wavread(datain);

% damping factor
% lower the damping factor the smaller the pass band
damp = 0.05;
% min and max centre cutoff frequency of variable bandpass filter
minf=500;
maxf=3000;
% wah frequency, how many Hz per second are cycled through
fw = 2000; 

% Function call
dataout = mywahwahfunction_2(x, fs, damp, minf, maxf, fw);

% Output as wave file
wavwrite(dataout, fs, N, 'out_wah.wav');

figure(1)
hold on
plot(x,'r');
plot(yb,'b');
title('Wah-wah and original Signal');