function z = P(x,y,r)

z =x*y-r*((1/(3-(x^2+y)))+(1/x^2)+(1/y^2));