Draw a phase portrait of a two-variable ODE system in Julia

Afterwards you can integrate a two-variable ODE system with DifferentialEquations.jl, find its fixed points, and draw its phase portrait over a direction field.