水波

/* h = A*sin(2*pi*W*t)*/
  int dx, dy; /* Disturbance Point */
  float s, v, W, t;
  int i, j;

  dx = WSIZE >> 1;
  dy = WSIZE >> 1;
  W = 0.3;
  v = -4; /* Wave speed */

  for (i=0; i<WSIZE; i++)
    for (j=0; j<WSIZE; j++)
    {
      s = sqrt((j-dx)*(j-dx) + (i-dy)*(i-dy));//距离中心点越远s越大
      wt[i][j] += 0.1;
      t = s/v;
      water[i][j] = A*sin(2*PI*W*(wt[i][j]+t));//距离中心点相同的点具有相同的t值，也就具有了相同的振幅。也就是整个面呈现正形波状态。
            //针对同一个点，由于t值一直相同，故增加一个值，这样同一个点的波形在变化，所有的点增加
            //相同的值。
    }