用DSolve解三维波动方程的柯西问题

Question

如上图所示的方程用DSolve求解....

DSolve[{D[u[x, y, z, t], {t, 2}] == 
   a^2 (D[u[x, y, z, t], {x, 2}] + D[u[x, y, z, t], {y, 2}] + 
      D[u[x, y, z, t], {z, 2}]),
  u[x, y, z, 0] == 0,
  (D[u[x, y, z, t], t] /. t -> 0) == x^2 + y z},
 u[x, y, z, t], {x, y, z, t}]

上面这段程序解不出来，请帮忙看看哪里弄错了⁄(⁄ ⁄•⁄ω⁄•⁄ ⁄)⁄