function z=sostiff1(t);
z=1+exp(-6*t);
