How does one shade the basin(s) of attraction of a phase plot in *Mathematica*? I have been trying to do this using the system

$$\begin{align\*}
\dot x &= y\\
\dot y &= -9\sin(x) - 0.20y
\end{align\*}$$

but I have gotten nowhere.