function f = f14_4_8(t,x)
f     = [10*(x(2)-x(1));  
	 28*x(1) - x(2) + x(1)*x(3);
	 -8*x(3)/3 + x(1)*x(2)];
