Best way to write this function in matlab while using eps in the denominator to avoid zero

Here's what I had z = (xx.*yy.*(xx.^2-yy.^2))/(eps(xx.^2+yy.^2)); but I kept getting a rank deficient error

[KSSV]

Your denominator is x^2+y^2, it will be zero for x = y = 0. You have to skip these values. Post your full code.

Answer 1

[Walter Roberson]

z = (xx.*yy.*(xx.^2-yy.^2)) ./ (eps+(xx.^2+yy.^2));