source("maysingle.R")
source("parabolasingle.R")
source("may.R")


r = 2.8
k = 100
n0 = 1
lag = 0
nostep = 100
noise =0.05

#simple dynamics
a=maysingle(r,k,noise,n0,nostep)
plot(a,pch=19,type="b")

#state space rendering of the simple dynamics
parabolasingle(r,k,noise,n0,nostep,4)

#go through different R values
may(1.1,4,0.1, k, noise,n0, nostep,2,0.5)