t=0:0.01:10;
y5=10.*sin(t).*exp(-t)-5.*sin(t).*exp(-0.5.*t);
plot(t,y5)
figure
y11=[zeros(1,250) -t(1:250) t(1:251)-2.5 zeros(1,250)];
plot(t,y11)


s=tf('s')
G1=1/(s+1)
G2=1/(s+2)
G3=2/s
G4=1/(s+2)
G5=zpk([],[1 -i i],-10)
G6=0.5
G7=0.1
G8=-1
G28=parallel(G2, G8)
G12=series(G1,G28)
G9=1
G39=parallel(G3,G9)
G34=series(G39,G4)
G34=minreal(G34)
G10=-1
G1234=G12-G34+G10
G12347=feedback(G1234,G7)
G123457=series(G12347,G5)
GF=feedback(G123457,G6)
M=minreal(GF)
c=zero(M)
p=pole(M)
ge=dcgain(M)


f=5
syms t s
y45=10*cos(2*pi*f*t+pi/2)*exp(-0.5*t)
figure
ezplot(y45)

y85=1-2*exp(-2*t)+exp(5*t)
tfy85=laplace(y85)
pretty(tfy85)

y85=(exp(-2*s))/(s^2+5*s+6)
itfy85=ilaplace(y85)