function z=fnsys1(t,U);
z=[2*U(1)+9*U(2),U(1)+2*U(2)];
