function f = func4(x,a)
% Defines the function in dx/dt = func4(x,a), a : parameter
% b is a constant when computing the curve of equilibria x(a)
b = 0.89;

f = b + a + (a^2)*x + a*x.^2 - x.^3;