(a)
\begin{equation}
\left\{
\begin{array}{**lr**}
1-xy=0 & \\
x-y^{3}=0\\
\end{array}
\right.
\end{equation}
\begin{equation}
\left\{
\begin{array}{**lr**}
xy=1 & \\
x=y^{3}\\
\end{array}
\right.
\end{equation}
\begin{equation}
\left\{
\begin{array}{**lr**}
x=1 & \\
y=1\\
\end{array}
\right.
\end{equation} or \begin{equation}
\left\{
\begin{array}{**lr**}
x=-1 & \\
y=-1\\
\end{array}
\right.
\end{equation}
Therefore, the critical points are (1,1) and (-1,-1)
(b)
The Jacobian matrix of the vector field is:
\begin{align*}
J &= \begin{bmatrix}
-y & -x \\
1 & -3y^{2}
\end{bmatrix}\\
~\\
J(1,1) &= \begin{bmatrix}
-1 & -1 \\
1 & -3
\end{bmatrix}\\
~\\
J(-1,-1) &= \begin{bmatrix}
1 & 1 \\
1 & -3
\end{bmatrix}
\end{align*}
(c)
\begin{align*}
For (1,1), let A&= \begin{bmatrix}
-1 & -1 \\
1 & -3
\end{bmatrix}\\
~\\
A-\lambda I &= \begin{bmatrix}
-1-\lambda & -1 \\
1 & -3-\lambda
\end{bmatrix}\\
~\\
det(A-\lambda I) &=(\lambda+3)(\lambda+1)+1=0\\
~\\
\lambda_{1} &= \lambda_{2} = -2 \\
~\\
Then \ the \ system \ has \ a \ stable \ improper \ node \ at \ (1,1) \\
~\\
For (-1,-1), let A&= \begin{bmatrix}
1 & 1 \\
1 & -3
\end{bmatrix}\\
~\\
A-\lambda I &= \begin{bmatrix}
1-\lambda & 1 \\
1 & -3-\lambda
\end{bmatrix}\\
~\\
det(A-\lambda I) &=(\lambda+3)(\lambda-1)-1=0\\
~\\
\lambda = -1 \pm \sqrt{5} \\
~\\
Then \ the \ system \ has \ a \ unstable \ saddle \ point \ at \ (1,1) \\
\end{align*}
(d) In the attachment.