octave> h=gaussian1d(64,32,3.0);
octave> drum(h*h',h*h',100,1,0.005,3); % Click on picture to see animation!

Frames 1-23*

Frames 24-47*

*Thanks to Anoop Menon.