function f = func2(x,r)
% Defines the function in dx/dt = func2(x,r), r : parameter
f = r*x + x.^3 - x.^5;