function z=fnsys2(t,U);
z=[U(2),6*cos(t)-U(1)];
