D[omega J1(omega sqrt(x^2+y^2))/(2pi sqrt(x^2+y^2)), {x, 3}]