After further coding I solved my question using this:
f=1000;
n=10;
T=1/f;
t=(0:T/100:n*T);
s = sin(2*pi*t*f);
e = (-expm1(-t*500));
e2 = (exp(-t*300));
tn = linspace(0,5*T,1001);
n = s .* e;
ts = linspace(5*T,8*T,1001);
tm = linspace(8*T,10*T,1001);
m = e2 .* s;
plot([tn ts tm],[n s m]);
However, if there is a way to clean this up I would be interested.