(** nessa questao a função vale uma expressao quando x difere de 0 e vale 1 quando x=0 . Eu tentei diferenciar esses dominios separando a função em 3 mas nao funcionou, o resultado esperado é esse e ficou bem diferente:

Plot3D[Sin[x y]/(x y), {x, -5, 5}, {y, -5, 5},PlotRange->All]