[inwave, fs] = wavread('800Hz.wav'); % Import wave file '800Hz.wav'
Modfreq = 6; % Vibrato function parameter, rate of modulating LFO in Hz
Width = 0.005; % Vibrato function parameter, basic delay in seconds of inwave

outwave=vibrato(inwave,fs,Modfreq,Width); % Call vibrato function

wavwrite(outwave, fs, '800Hz_with_vibrato.wav'); % Create wave file from function output