f := proc(x,y) if x = 0 and y = 0 then RETURN( 0 ) else RETURN( x*y*(x^2-y^2)/(x^2+y^2) ) fi; end; assume(a>0, b >0): f(a,b); fx := unapply(diff(f(x,y),x),x,y); fx(0,y); fy := unapply(diff(f(x,y),y),x,y); fy(x,0);