>19
スマソ。でももう一度。
俺的改良版。理論的に正しいかは不明。。
float _RT2 = 1.0f / 1.41423563f;
for(int i = width + 1; i < p.length - width - 1; i++) {
 velocity[i] += WAVE * (-(_RT2 + 1) * 4 * p[i]
  + p[i - width] + p[i - 1] + p[i + 1] + p[i + width]
  + _RT2 * (p[i - width - 1] + p[i - width + 1] + p[i + width + 1] + p[i + width - 1]));
 if(velocity[i] < VMIN)
  velocity[i] = VMIN;
}
速度の上限決めてやると、もっさりとしていい感じ。