function f = func1(x,r)
% Defines the function in dx/dt = func1(x,r), r : parameter
f = -x + r*tanh(x);