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


