Simulation of the phase plane

The solution's form for the ODE $$\frac{d\vec{r}(t)}{dt\;\;\;\;} = k\;\vec{r}(t)$$ can be generalized like in this diagram:



Exist some program or some way of adjust some program of math for study the behavior of the solution through of the eigenvalues as parameter?