function dv = potential_eqn(t,v)
global a
dv = -v.*(v-a).*(v-1);