I wrote one of these in matlab for a partial differential equations class a year ago.
I had fun with that, it's fun to mix neumann with direchlet boundary conditions, you get some cool results.
Haha, I'd just stare at the output for long stretches of time. Good class. Waves are cool.