Erro na função PERSP

Olá, pessoal.

Estou com um problema tentando plotar a função f(x,y)=1/(x+y). Eis o meu código:

Y <- seq(-5, 5, length= 50)
X <- seq(-5, 5, length= 50)

Z = function(X, Y) +1/(X+Y)

Z2 <- outer(X, Y, Z)

persp(X, Y, Z2, theta = 30, phi = 30, expand = 0.5, col = "brown", ltheta = 120, shade = 0.75, ticktype = "detailed", xlab = "X", ylab = "Y", zlab = "Z", main="Z=X+Y")

O que eu percebi é que quando eu chamo Z2, algumas entradas da matriz estão definidas como Inf (óbvio, devem ser as combinações que fazem Z tender ao infinito).

Alguma ideia de como posso plotar este gráfico?