Если достаточно приближенной модели, можно сделать так:
r1=sqrt(x^2+y^2)
a=A*p/sqrt(p^2+x^2+y^2)
z=a*sin(k*r1-w*t)
где p - размер области в окрестности нуля с насильственно ограниченной амплитудой,
A - максимальная амплитуда (в нуле)
k - пространственная частота (т.е. 2*pi/k - расстояние между волнами)
w - временнАя частота (2*pi/w - период).
w/k - скорость распространения.