function f = f14_6_8(t,x)
a = 1; b = 0.2; c = 0.2;
f = [-x(2)-x(3); 
     x(1) + b*x(2); 
     c + x(3)*(x(1)-a)];


