num=[10 155 75] den=[1 161 6660 56500 50000] r_den=roots(den) C1=polyval(num,-1)/(-1+100)/(-1+10)/(-1+50) C2=polyval(num,-10)/(-10+100)/(-10+1)/(-10+50) C3=polyval(num,-100)/(-100+10)/(-100+1)/(-100+50) C4=polyval(num,-50)/(-50+10)/(-50+1)/(-50+100) [r,cc]=residue(num,den); r,cc