%
% Heart Beat Dynamics Case Study and Figure 4.5
%
Td=0.4;
t=0:0.01:Td;
num=[1 12 -10]
den=[1 22 20 20]
[k,p]=residue(num,den)
x1=k(1)*exp(p(1)*t)+2*abs(k(2))...
*exp(real(p(2))*t).*cos(imag(p(2))*t+angle(k(2)));
% ... stands for the concatenation in MATLAB
plot(t,x1)
xlabel('time of the systolic period in [sec]')
ylabel('heart mulsle fiber normalized length')
grid
print -deps fig4_5.eps