% Vypocet rozkladu nad R pomoci rozkladu nad C uzitim prikazu 'residue' v MATLABu
% Dosla, Kuben, Priklad 3.14b):

[R,P,K] = residue([1 1],[1 0 3 0 2 0])
R =
   0.2500 + 0.3536i
   0.2500 - 0.3536i
  -0.5000 - 0.5000i
  -0.5000 + 0.5000i
   0.5000          
P =
        0 + 1.4142i
        0 - 1.4142i
        0 + 1.0000i
        0 - 1.0000i
        0          
K =
     []

Q1=conv([1,-P(1)],[1,-P(2)])
Q1 =
    1.0000         0    2.0000
Q2=conv([1,-P(3)],[1,-P(4)])
Q2 =
    1.0000         0    1.0000
P1 = R(1)*[1,-P(2)]+R(2)*[1,-P(1)]
P1 =
    0.5000   -1.0000
P2 = R(3)*[1,-P(4)]+R(4)*[1,-P(3)]
P2 =
   -1.0000    1.0000
% Vysledny rozklad nad R:
% R(5)/(x-0) + P1(x)/Q1(x) + P2(x)/Q2(x) = 1/(2x)+ (x-2)/(2(x^2+2)) + (-x+1)/(x^2+1)
diary off