 Solution

With $f(x,y)=-x{y}^2-y^5$, write $t=x_{n+1}=x_n+h$ and $z=y_{n+1}$. The Backward Euler method requiresor the monic quintic polynomial equationConstruct its $5\times 5$ companion matrixIts characteristic polynomial is $\det (zI-C_n)=p(z)$. Thus the polynomial roots needed for this update are exactly the eigenvalues of $C_n$.

At each step, assemble $C_n$, compute its eigenvalues, for example with the QR algorithm, and choose its unique positive real eigenvalue as $y_{n+1}$. To justify this selection, letIf $y_n>0$, then $F(0)=-y_n<0$, $F(y_n)=h{y}_n^5+ht{y}_n^2>0$, andThe intermediate value theorem and strict monotonicity prove there is exactly one positive polynomial root, located in $(0,y_n)$. Starting at $y_0=1$ therefore makes every update well defined, positive and decreasing. In numerical arithmetic, select the nearly real eigenvalue in this interval and check the polynomial residual $|F(z)|$; a final Newton root-finding iteration can refine it. This reduces each update to the requested eigenvalue problem.