function vf = vfield2(t,vec)
% Defines the vector field in computing the bifurcation curves for f=func2
% vf(1) is df/dx, vf(2) is -df/dr
r = vec(1);
x = vec(2);

vf(1) = r + 3*x^2 - 5*x^4;
vf(2) = -x;
vf = vf';