Displaying 1 result from an estimated 1 matches for "alphamin".
Did you mean:
alphai
2010 Jun 24
1
help, bifurcation diagram efficiency
...06748*
sin(2*pi*t) + 0.1063696962*sin((9/4)*pi*t) + 0.02812403009*sin((5/2)*pi*t)))
}
ModBraaksma = function(t, n, p)
{
dx.dt = (1/0.01)*(n[2]-((1/2)*n[1]^2+(1/3)*n[1]^3))
dy.dt = -(n[1]+p["alpha"]) + 0.032 * s_of_t(t)
list(c(dx.dt, dy.dt))
}
initial.values = c(0.1, -0.02)
alphamin = 0.01
alphamax = 0.02
alphas = seq(alphamin, alphamax, by = 0.00001)
TimeInterval = 100
times = seq(0.001, TimeInterval, by = 0.001)
plot(1, xlim = c(alphamin, alphamax), ylim = c(0,0.3), type = "n",xlab
= "Values
of alpha", ylab = "Approximate loop size for a lim...