function f = f14_4_9(t,x)
f     = [3*x(1) + x(3) - 3 - x(1)*x(3);  
	   x(1) - 1 + x(2)*x(3);
	   x(2) + x(3) - 2 + x(1)*x(2) - x(2)*x(3)];
