%% Ricker Wavelet function
function [ws] = ricker_wave( t,nu0 )

sigma = 1/(pi*nu0*sqrt(2));
t0 = 6*sigma;
ws = 200*( 1 - 2*pi^2*nu0^2*(t-t0).^2 ) .* exp( -pi^2 * nu0^2 * (t-t0).^2 );

return;

