I want to calculate extremes of certain multivariable function $f(x,y)=(6−x−y)x^2y^3$. After solving system of derivatives $f_x=0$ and $f_y=0$ I got something like this:
$P_1=(x,0),x\in \mathbb R$
$P_2=(0,y),y\in \mathbb R$
$P_3=(2,3)$
First two conditions are satisfied with infinite number of $x$ and $y$. How am I supposed to act in such situation? Do I have to check the first two points in some way? If so, how should I do this?


Based on your comments I think that you have a problem with the second derivative test. Note that the Hessian of $f$ is
$$\begin{pmatrix} 12 y^3 - 6xy^3 - 2y^4 & 36 xy^2 - 9x^2 y^2 - 8xy^3 \\ 36 xy^2 - 9x^2y^2 - 8xy^3 & 36x^2 y - 6x^3 y - 12 x^2 y^2 \end{pmatrix}.$$
Now we want to check what happens at $P_1 = (x, 0)$ with $x \in \mathbb{R}$. Plugging this in yields the zero matrix
$$\begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix}.$$
Since this matrix has determinant $0$, the second derivative test fails. Similarly, for $P_2 = (0, y)$ with $y \in \mathbb{R}$, we get
$$\begin{pmatrix} 12y^3 - 2y^4 & 0 \\ 0 & 0 \end{pmatrix}$$
for every $y \in \mathbb{R}$. Note that this matrix also has determinant $0$ (for every $y$). So you have to use other methods for determining what happens there. Some ideas were already given in the answers, have a look also e. g. here.