Displaying 1 result from an estimated 1 matches for "modbraaksma".
2010 Jun 24
1
help, bifurcation diagram efficiency
...5815*sin((1/4)*pi*t) - 1.697435362*sin((1/2)*pi*t) - 1.570296922*sin
((3/4)*pi*t) + 0.3247901958*sin(pi*t) + 0.7962749105*sin((5/4)*pi*t) +
0.07812230515*sin((3/2)*pi*t) - 0.3424877143*sin((7/4)*pi*t) - 0.1148306748*
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 = 10...