clear all
clc
close all

syms s % cria variável simbólica

Ns = 5*s^2+13*s+10 
coefNum = [5 13 10]
Ds = (s+1)^2*(s+2) % Declara denominador
expand(Ds) % Expande denominador para imprimir os coeficientes na tela
coefDen = [1 4 5 2]

% A função residue aceita como entrada os coeficientes do numerador e
% do denominador de uma razão de polinômios,
% expande em frações parciais e devolve o coeficiente do numerador de cada
% parcela, os polos de cada parcela e o ganho DC caso o numerador tenha um
% grau maior ou igual ao do denominador.
% No caso de polos múltiplos, os coeficientes aparecem na ordem crescente
% de potência do denominador
[coef,polos,ganhoDC]=residue(coefNum,coefDen)

help residue