x, y = var('x y')
plot3d(sin(pi*(x^2+y^2))/sqrt(x^2+y^2),
       (x,-1.5,1.5),(y,-1.5,1.5),
       mesh=True,color='yellow')