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