if( !require('phaseR') ){
    install.packages('phaseR')
}
library(phaseR)

my_yprime = function(t, y, parameters){
    y_prime = 5 - 3*sqrt(y)
    list ( y_prime )
}

#postscript("../../WriteUp/Graphics/Chapter2/chap_2_sect_7_prob_5_plot.eps", onefile=FALSE, horizontal=FALSE)

diff_eq_params = c()

logistic.flowField <- flowField(my_yprime, x.lim=c(0, 6), y.lim=c(0.5, 6), parameters=diff_eq_params, points=21, system="one.dim", add=FALSE)

#dev.off()