function f = fcnCP5p19(x,gamma,delta)
f = [gamma*x(1)*x(2)-x(1)*(1+x(2)); -x(1)*x(2)+(delta-x(2))*(1+x(2))];
