Dear all, I hope you are doing well.

Previously, I have tried to submit this question, however, I forgot to include my attempts which now I do so. Attached below, you may find my notebook with my doubts. It should be remarked that I have also read Mathematica's documentation. However, this latter did not help me with the Taylor and Laurent series expansion of multivariable functions.

With that said, suppose that we have the following multivariable function

$\displaystyle \psi(x,y)=\frac{-((2 (132 + 56 x^4 - 382 y + 394 y^2 - 171 y^3 + 26 y^4 + 6 x^3 (-47 + 31 y) + x^2 (608 - 782 y + 240 y^2) + x (-602 + 1130 y - 716 y^2 + 153 y^3)))}{((16 + 4 x^2 + 8 x (-2 + y) - 16 y + 5 y^2)^2 (5 x^2 + (3 - 2 y)^2 + x (-6 + 4 y))^2))}$

Based on the above, I ask the following questions:

1. How may I construct a 3D plot of $\psi(x,y)$ so that positive and negative values are represented by red and blue colors, respectively? As one can readily see from the notebook attached below, I have considered the ColorFunction command to represent the positive and negative values of $\psi(x,y)$. However, I am not quite sure whether this is correct. As we may observe, $\psi(x,y)$ is a negative function. Hence, should not we observe a predominantly blue color in the 3D plot of $\psi$?
2. Is there a way to determine whether or not $\psi(x,y)$ is symmetric using Mathematica?

3. I have read Mathematica's documentation, however, I did not find a proper approach to performing series expansion of multivariable functions like $\psi$. Consequently, How may one perform the Taylor or Laurent series of this function in the neighborhood of the point $\displaystyle (x,y)=(2,1)$ through Mathematica?

Thanks once more and I look forward to hearing from you.

Attachments:

MULTIVARIABLE_FUNCTION.nb