Differential equations in ecology (Levins function): how to create a FlowField diagram in R

Question

I found a nice looking function in R to return solutions to the Levins equation (useful in ecological modelling) for specific initial parameter values, a time series and an