Date: 10 May 2013

a = 0.1
vnullrightmaxv: (1 + a + sqrt(1 - a + a^(2))) / 3
vnull(x) = (((-x) * (x - a)) * (x - 1))
epsilon = 0.01
gamm = 2
wnull(x) = x / gamm
vmax = 1.1
vmin = -0.3
vnullrightmax: vnull(vnullrightmaxv)
wmax: (vnullrightmax * 1.2)
wmin: ((-vnullrightmax) * 0.3)
DP2: (0.15 + a / 2, Min[wnull(0.15 + a / 2), wmax] / 2)
DP3: (0, wmax / 2)
DP4: ((2 * vmin / 3), 0)
DP1: (vnullrightmaxv, 0)
DPN1: ((0.3 * vnullrightmaxv) + 0.7, vnull((0.3 * vnullrightmaxv) + 0.7))
DPN2: (vmin / 2, vnull(vmin / 2))
wnullinv(x) = (x * gamm)
DPN3: (wnullinv((2 * wmax / 3)), (2 * wmax / 3))
DPN4: (wnullinv((2 * wmin / 3)), (2 * wmin / 3))
dvdt(v, w) = (((-v) * (v - a)) * (v - 1)) - w
dwdt(v, w) = (epsilon * (v - (gamm * w)))
vecfieldmag(v, w) = sqrt((dvdt(v, w) / vmax)^(2) + (dwdt(v, w) / wmax)^(2))
dvdtnorm(v, w) = dvdt(v, w) / vecfieldmag(v, w)
dwdtnorm(v, w) = dwdt(v, w) / vecfieldmag(v, w)
dirvect1 = Vector(DP1, DP1 + (dvdtnorm(DP1), dwdtnorm(DP1)) / 10)
dirvect2 = Vector(DP2, DP2 + (dvdtnorm(DP2), dwdtnorm(DP2)) / 10)
dirvect3 = Vector(DP3, DP3 + (dvdtnorm(DP3), dwdtnorm(DP3)) / 10)
dirvect4 = Vector(DP4, DP4 + (dvdtnorm(DP4), dwdtnorm(DP4)) / 10)
dirvectnull1 = Vector(DPN1, DPN1 + (dvdtnorm(DPN1), dwdtnorm(DPN1)) / 10)
dirvectnull2 = Vector(DPN2, DPN2 + (dvdtnorm(DPN2), dwdtnorm(DPN2)) / 10)
dirvectnull3 = Vector(DPN3, DPN3 + (dvdtnorm(DPN3), dwdtnorm(DPN3)) / 10)
dirvectnull4 = Vector(DPN4, DPN4 + (dvdtnorm(DPN4), dwdtnorm(DPN4)) / 10)
v01 = 0.23603603603603612
I1: (v01, 0)
v02 = 0.06261261261261268
I2: (v02, 0)
A = Intersect(vnull, wnull)