id3 := .001;f := (x1,x2) ->((1+id3)*x1+(1-id3)*x2)*exp(-(x1+(1+id3)*x2)) -x1;f(x1,x2);dfdxA := diff(f(x1,x2),x1);dfdxA := (x1,x2) -> (1+id3)*exp(-x1-(1+id3)*x2)-((1+id3)*x1+(1-id3)*x2)*exp(-x1-(1+id3)*x2)-1;dfdxB := diff(f(x1,x2),x2);dfdxB := (x1,x2) -> (1-id3)*exp(-x1-(1+id3)*x2)-(1+id3)*((1+id3)*x1+(1-id3)*x2)*exp(-x1-(1+id3)*x2);dfdxB(x1,x2);