//Example 3.14
//Partial fraction expansion for distinct real roots
clear;
clc;
//------------------------------------------------------------------
// Transfer function
s=%s;
num=2;
p1=(s+1);
p2=(s+2);
p3=(s+4);
sys=syslin('c',num/(p1*p2*p3))

//Partial fraction expansion is: sys= r1/p1 + r2/p2 + r3/p3
//residue calculation
r1=residu(num,p1,(p2*p3))
r2=residu(num,p2,(p1*p3))
r3=residu(num,p3,(p1*p2))

disp([r1 r2 r3]',"Residues of the poles p1, p2 and p3 are")
disp([roots(p1), roots(p2), roots(p3)]',"Poles p1, p2 and p3 are at")
disp('k=[]')
//------------------------------------------------------------------