Solution (source code)

= Solution

Set $y=\dot x$. Then $\dot x=y$, $\dot y=a+bx-y-x^2-y^2$.