function z=sosys1(t);
z=[3*exp(5*t)-6*exp(-t),exp(5*t)+2*exp(-t)];
